AO-A041  864 


UNCLASSIFIED 


BOEING  commercial  AIRPLANE  CO  SEATTLE  WASH  F/G  9/5 

AIR  TRAFFIC  CONTROL  EXPERIMENTATION  AND  EVALUATION  WITH  THE  NAS— ETC (U) 
SEP  76  A D THOMPSON*  S C WILSON.  P F RIEDER  DOT-TSC-707-4 


06-44049 


FAA  RD-75-173-4 


I OF^ 

A0A04I864 


OA 04 1 864 


CO  a'r  traffic  control  experimentation  ano  evaluation 


i oo 

tH 

WITH  THE  NASA  ATS-6  SATELLITE 

Volume  IV:  Oata  Reduction  and  Analysis  Software 

o 

> 

A.D.  Thompson 

S.G.  Wilson 

i o 

P.F.  Rieder 

1 <2: 

. 

W.L.  Chu 
M.J.  Mardesich 
C.V.  Paulson 
P.  Alexander 

SEPTEMBER  1976 
FINAL  REPORT 


Document  is  available  to  the  U.S.  public  through  the 
National  Technical  Information  Service. 
Springfield,  Virginia  22161 


1 


^ Government  Accession  No 


pienT  s Catdioq  No 


'js-xiitiOi 


' Trrrran  <T  S • i nTTri"?  ■" 


MR IKAFFIC  CONTROL  hXl'fcRIMHNTATION  AND  F.VALUAT10N 
AITH  TKFIVASA  AIS-6  S^TELLITH'  5-  , 


AVITH  TKEJ^A  ArS-6  .satellite' 

Volume  ly.  Data  Reduction  and  Analysis  Software 


uOT-TSC-P.\-\-77-  22  .IV 


A.D./riiompson.  S.G./4'ilson 


/ilson,  P F-A* 
iilson,  P.  Alex 


IRieder,  W.L/hu/ 
lexander  ' ’ 


D6-  44049 


10  Work  Unit  No 


9 Perlorminq  Orqani/aiion  Name  and  Afldress 


Boeing  Commercial  Airplane  Company* 
P.O.  Box  .^707 
Seattle.  WA  98124 


FA711/R7106 


DOl-TSl-707-4 


12  Sponsonnq  Aqencv  Name  and  Address 

U.S.  Department  of  Transportation 
Federal  Aviation  Administration  ( J yy 

Systems  Research  and  Development  Service  ^ 

Washington.  DC  20.S91  j 

.5  Supoiemeniarv  Motel  (J.S.  Department  of  Transportation 
*Under  Transportation  Systems  Center 

contract  to:  Kendall  Square 

Cambridge.  MA  02142 


rl~'Tvtnro1  Ricmrl  enO  Pef.otiLu 

Final  Report,  , 
Sep.'tf7.?  to  Dec.  »75, 


t iy  t f ^ 


Software  used  for  the  reduction  and  analysis  of  the  multipath  prober,  modem  evaluation  (voice,  digital 
data,  and  ranging),  and  antenna  evaluation  data  acquired  during  the  ATS4i  field  test  program  is  described. 

Multipath  algorithms  include  reformatting  operations,  delay-spectra  time  histories.  delay-Doppler 
scatter  function' S(t,co)!  noise  determination  and  removal,  spread  calculations,  airborne  tape  analysis, 
and  uU)«r4c’tailed  processing  including  time-domain  analysis  and  various  integral  and  Fourier  operations 
oq  S(T,u) ).  Modem  and  antenna  evaluation  data  processing  software  includes  algorithms  for  the  deler- 
miitalion  of  ( 1 ) C/N^  and  multipath  interference  ratio,  S/I.  (2)  digital  data  bit-error  rates,  block  error 
statistics,  and  inter-error  spacing,  and  (3)  ranging  error  statistics  and  distribution.  Sample  outputs  are 
given.  Program  listings  and  other  information  arc  provided  in  an  auxiliary  software  data  package. 

The  report  consists  of  seven  volumes:  I Executive  Summary ; II  Demonstration  of  Satellite- 
Supported  Communications  and  Surveillance  for  Oceanic  Air  Traffic  Control;  III  Summary  of  IfS. 
Aeronautical  Technology  Test  Program;  IV  Data  Reduction  and  Analysis  .Software;  V Multipath 
Channel  Characterization  Test;  VI  Modem  Evaluation  Test;  VII  Aircraft  Antenna  Evaluation  Test. 


c/,v'  ^ 


1 7 Key  Words  (Suqqested  by  Authorlst  \ 

ATS-(i  Satellite.  Multipath.  Modem  Evaluation. 
Digital  Data.  Ranging,  Fast  Fourier  Transform. 
Data  Processing.  Algorithms.  Envelope  Detector 
Analysis.  Scatter  Function 


I 18  Oistrthutton  Statement 


Document  is  available  to  the  I'.S.  public  ihtougli  the 
National  Technical  Information  Service.  Spiingfield. 
Virginia  22l(il 


19  Secijfiiy  Cinss'i  lot  this  r^oorf* 

Unclassified 


20  Security  Classit  lot  this  paqet 

Unclassified 


21  No  ot  Paqas  1 22  Price* 


‘for  sale  by  the  National  Techn*cai  Information  Service  Sprmqfield  Virqima  221S1 


J 


AIR  TRAFFIC  CONTROL  FXPFRIMENTATION  AND 
LVALUATION  WITH  THE  NASA  ATS-6  SATELLITE 


FINAL  Rl  PORl 


I l^l^  report  eoMMNts  of  tne  lollowlne  \olumos. 


Volume  I 

F\eeutive  Siimm;ir\ 


X'oliime  II 

l)emoii>lr;ition  ol  S;itelliIe-Siipporteil  Comnuinieations 
aiiel  Surveillaiiee  lor  Oeeanie  Air  1 ralTie  C ontrol 


Volume  111 

Summary  of  L'.S,  Aeronautieal  Teelmolopy 
Test  Prosiram 


Volume  l\’ 

Data  Reeluetion  and  Analvsis  Software 


N'olume  V 

Multipath  ( haiinel  Charaeteri/ation  Test 


Volume  VI 

Modem  Evaluation  Test 


Volume  VII 

Aireraft  Antenna  Evaluation  Test 


II 


PREFACE 


I 


The  U S.  Uepartiiv.’iit  ol'  1 lansportation  ( DOT  I aeronautical  test  program  entitled  "Air  FralTie 
Control  Experimentation  and  Evaluation  with  the  NASA  ATS-(>  Satellite”  was  part  of  the  Integrated 
ATS-P  L-Band  Experiment.  The  overall  ,\ TS-(i  L-band  experiment  was  coordinated  by  the  NASA/ 
Cioddard  Space  Flight  Center  tCiSF'O  and  was  international  in  scope.  The  following  agencies  were 
participants  in  the  experiment:  NASA/Goddard  Space  Flight  Center;  DOT/Federal  Aviation  Admini- 
stration; DOT/Transportation  Systems  Center;  DOI'/U.S.  Coast  Guard;  DOC/Maritime  Administration; 
European  Space  Agency  (ESA);  and  the  Canadian  Ministry  of  I ransport  and  Department  of  Communi- 
cations. Each  participant  performed  tests  in  one  or  more  of  three  categories:  aeronautical,  maritime 
safety,  and  maritime  lleet  operations.  All  tests  were  conducted  in  accordance  with  an  overall  integrated 
test  plan  coordinated  by  NASA/GSFC. 

The  U.S.  DOT  Aeronautical  test  program  was  undei  the  direction  and  sponsorship  of  the 
Federal  Aviation  Administration.  Systems  Research  and  Development  Service  (SRDS),  Satellite 
Branch,  with  the  DOT'TSC  conducting  the  technology  tests  and  the  FAA/NAFEC  conducting 
the  ATC  demonstration  tests.  The  technology  tests  included  multipath  channel  characterization, 
modem  evaluation,  and  aircraft  antenna  evaluation.  Results  of  these  tests  are  presented  in  volumes 
III  through  VII,  and  the  results  of  the  ATC  demonstration  tests  are  presented  in  volume  II.  The 
DOT/TSC  test  program  was  supported  by  the  Boeing  Commercial  Airplane  Company  under  contract 
DOT-TSC-707.  Mr.  R.  G.  Bland  was  the  TSC  Project  Engineer  and  Contract  lechnical  Monitor. 

This  volume  describes  the  software  used  for  reduction  and  analysis  of  the  technology  test  data. 
All  work  described  was  performed  under  contract  DOT-TSC-707.  Sections  4.21  and  appendix  B of  this 
volume  were  prepared  by  CNR,  Inc.,  under  subcontract  to  the  Boeing  Commercial  Airplane  C'ompany. 
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1.  INTRODUCTION 


Work  performed  under  eontruct  DOT-TSC-707  includes  all  software  development  (for  both 
quick-look  and  detailed  data  reduction  and  analysis),  data  processing,  analysis,  and  reporting  of  results 
for  the  following  U.S.  DOT  aeronautical  technology  tests; 

Multipath  channel  characterization 

Modem  evaluation 
Voice 
Data 
Ranging 

Antenna  evaluation 

This  report  gives  a detailed  description  of  the  software  and  associated  hardware  used  for  the 
reduction  and  analysis  ot  data  acquired  during  the  ATS-6  program  for  the  above  tests.  Sections  „ 
through  5 describe  software  and  hardware  for  reduction  and  analysis  of  multipath  prober  data;  sec- 
tions 6 through  1 1 provide  equivalent  information  for  the  modem  and  antenna  evaluation  tests. 


Program  listings  and  other  information  relevant  to  the  software  described  are  provided  in  an 
auxiliary  software  data  package.  The  software  data  package  was  delivered  to  DOT/TSC  during  April 
1975.  Software  developed  after  April  1975  will  be  delivered  upon  completion  of  the  DRandA  task. 


2.  MULTIPATH  DATA  PROCESSING  FUNCTIONAL  SEQUENCE 


2.1  SYSTEM  BLOC  K DIAGRAM 

Figure  2-1  depicts  the  processing  steps  involved  in  tlie  reformatting,  reduction,  and 
analysis  of  the  recorded  multipath  SACP  signal  arrays  and  the  aircraft-transmitter  parameter  tape. 

Received  SACP  signals*  are  direct  recorded  in  a standard  telemetry  analog  format  (serial 
PCM,  NRZ-L).  On  the  analog  source  tape,  two  tracks  are  in  general  occupied  by  the  multipath 
data  while  another  of  the  seven  available  tracks  contains  an  IRIG-A  modulated  time  code.  At  the 
Boeing  ground  station  facility,  the  data  and  time  tracks  are  initially  processed  by  telemetry  front 
end  (TFE)  equipment,  which  for  this  particular  application  routes  the  played-back  serial  signals 
through  its  PCM  subsystem  to  the  programmable  data  distributor  (PDD),  The  PDD  merges  time 
words  with  the  data  and  distributes  the  information  to  one  or  both  of  the  PDP  1 1/45  computer 
I/O  buses.  The  dual  PDP  computer  system  is  called  upon  to  perform  three  basic  functions: 

a.  Conversion  of  the  analog  recorded  data  tapes  into  digital  format  computer-compatible 
tapes.  The  digital  tapes  are  1/2-in.  width  with  nine  parallel  tracks  written  in  1600-bpi 
density. 

b.  Quick-look  processing  of  the  multipath  data.  This  consists  of  a playback  real-time 
analysis  of  the  recorded  tap  outputs  and  produces  individual  tap  Irequency  spread, 
delay  spectra,  and  receiver  parameter  dumps  for  all  pertinent  test  data. 

c.  Calculation  of  the  time-ordered  delay-spectra  arrays  that  are  used  to  generate  the 
time  history  of  the  multipath  channel’s  delay  spectra.  The  normalization  and  three- 
dimensional  plottingof  this  data  are  performed  in  the  ( DC  b(i00. 

The  prober  data  is  analyzed  in  detail  by  the  CDC  6600  computer.  This  analysis  provides 
a comprehensive  characterization  of  the  multipath  channel  for  horizontal  and  vertical  polariza- 
tion, gathered  over  a down-looking  antenna.  The  primary  output  of  this  routine  is  the  delay- 
Doppler  scatter  function  of  the  channel.  Also  included  are  the  channel's  time-frequency  auto- 
correlation function . total  scattered  intensity ; delay  spectrum ; Doppler  spectrum : frequency 
autocorrelation  function;  time  autocorrelation  function;  parameter  spread  values;  tap  1,  Q,  and 
phase  distributions;  goodness-of-fit  test  on  the  tap  distributions;  tap  cross-correlation;  tap  I and  Q 
dependency  tests;  and  the  tap-gain  autocorrelation  function. 

In  addition  to  the  SACP  receiver  data  reduction,  this  report  also  includes  details  pertain- 
ing to  software  designed  to  reduce  and  analyze  the  data  tape  containing  airplane  and  transmitter 

'complex  I.Q  outputs  from  two  delay  tap  correlator  banks,  each  consisting  of  1 1 2 elements. 
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parameters.  The  airplane  system  ami  transmitter  parameters  are  recorded  in  a time  multiplex-mode  of 
seven-track  (six  data,  one  parity  I tape  that  is  directly  compatible  with  the  Boeing  CDC  6600.  Data 
extracted  from  the  airborne  system  parameters  tape  will  serve  as  a data  collection  integrity  measure 
and  will  be  used  to  augment  the  logged  llight  test  data  for  normalization  purposes. 


SYSThM  DE  SCRIPTION 

In  this  section,  details  pertaining  to  the  computer  systems  used  to  analyze  the  multipath  data 
are  provided. 


-.2. 1 Boeing  Flight  Test  PDP  1 1/45  Telemetry  Ground  Station 

Tlie  Boeing  Flight  Test  Center  ground  station  is  used  to  convert  analog  PCM  telemetry  data  to 
formatted  digital  information,  to  perform  quick-look  reduction  including  delay-spectra  and  tap 
frequency-spread  calculations,  and  to  set  up  the  data  base  for  the  delay-spectra  time  history  calculation. 

The  hardware  configuration  of  the  ground  station  is  depicted  in  figure  2-2.  Vertical  and  hori- 
zontal data  and  time  code  are  read  simultaneously  on  the  analog  tape  transport.  These  data  are  routed 
to  various  components  of  the  telemetry  front  end  (TFE)  by  the  EMR  2793  data  matrix,  which  is  man- 
ually set  up  at  the  start  of  a test  series. 

Multipath  data  is  routed  to  the  EMR  720  bit  synchronizer  where  the  PCM  bit  stream  is  dig- 
itized and  synchronized  to  a clock  signal.  The  EMR  720  output  is  input  to  the  frame  synchronizer, 
which  searches  for  the  frame  sync  bit  pattern  that  heads  each  frame  of  data.  This  action  allows  for  the 
correlation  of  specific  stream  bit  segments  to  real  data  elements  that  can  be  processed  by  computer. 

Frame  synchronizer  output  is  processed  by  the  programmable  data  distributor  (PDD),  which 
removes  blank  words  and  merges  time  code  information  into  the  data  stream.  The  system  setup  con- 
troller determines  the  source  of  the  merged  time  code  for  the  PDD  and  interrupts  for  the  processor. 

Test  condition  liming  control  and  time  code  data  are  generated  by  the  Datum  tape  search  unit 
and  time  code  translator.  A start  and  stop  time  is  set  in  this  unit,  which  is  used  to  delimit  the  test 
condition  interval.  When  the  operator  has  the  TFE  set  up.  system  control  is  transferred  to  this  unit. 

Most  of  the  TFE  components  are  set  up  by  control  signals  from  the  computers  that  interface 
to  the  TFE  through  an  array  of  2763  interfaces.  These  signals  arc  generated  through  interpretation  of 
commands  either  made  on  the  DEC  writer  terminal  by  the  operator  or  read  in  from  punched  cards  by 
the  card  reader. 


CONTIOL  SlfiNM. 


Igure  2-2.  Boeing  Flight  Test  Data  Processing  EMR  Ground  Station 


Formatted  telemetry  data  from  the  FDD  is  input  to  the  PDF  I 1/45  memory  via  a 2763  DMA 
channel.  The  data  is  processed  and  written  to  the  TMA-2. 

Following  the  reduction  of  a test  condition,  the  digital  tape  files  containing  level  quick-look 
data  are  read  back,  with  graphic  output  produced  on  the  Tektronix  4012  CRT  display  and  printed 
summaries  generated  on  the  Gould  4800  printer. 

The  RK05  and  RP02  are  used  to  store  programs  for  the  data  processing  and  system  operation. 
2.2.2  CPC  6600  Computer  System 

A CDC  6600  computer,  with  KRONOS  2.1  operating  system,  is  used  to  process  the  multipath 
data  tapes. 

Figure  2-3  gives  the  data  How  for  the  airborne  tape  processing.  Program  inputs  consist  of  a 
seven-track  tape  and  control  data  cards.  The  program,  which  is  run  in  batch  mode  on  the  6600  com- 
» puter.  produces  numerical  output  and  a plot  tape  that  is  processed  off-line  by  the  CalComp  plotter. 

This  plotter  was  chosen  for  the  quality  and  size  of  its  plot.  Tlie  processing  program  is  written  in 
FORTRAN  IV.  To  format  the  plot  tapes.  CalComp  plot  routines  from  the  system  library  are  used. 

* Figure  2-4  gives  the  data  How  for  the  delay-spectra  time  history  and  the  detailed  analysis  proc- 

essing. Program  inputs  consist  of  the  reformatted  tape  and  control  data  cards.  The  program  is  run  in 
batch  mode.  Temporary  storage  on  disk  is  used  to  record  the  input  data.  Program  outputs  include  a 
listing,  a save  tape  (for  the  scatter  function),  and  a plot  tape  that  is  processed  off-line  on  a SC  4020 
plotter.  The  graph  produced  by  the  SC  4020  has  resolution  and  size  adequate  for  multipath  output. 
The  processing  programs  are  written  in  FORTRAN  IV  and  CDC  6600  assembly  language  (27?  of  the 
coding).  The  subroutines  used  from  the  system  library  are  relative  to  random  access  (R/A)  input/ 
output  operations  and  plot  tape  formatting. 
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.i.  MULTIPATH  ALGORITHM  LXECUTION  SEQUENCE 


I muros  -G1  ;iiul  depict  the  aljtoriihm  execution  xctiuence  for  processing  the  multipath 
channel  data  and  the  recener  s>  stent  paranieters.  Detailed  descriptions  of  the  indicated  block  opera- 
tions are  gisen  in  section  4,  Processing  of  the  airborne  system  paranieters  tape  is  carried  out  independ- 
ently from  this  seciiience  and  is  discussed  in  section  4.23. 


3.1  AUiORl  niM  I XLCUTION  Sl.QUl.NCL 

The  operations  described  in  figures  3-1  and  3-2  are  loosely  categorized  as  A/D  reformatting, 
quick-look,  delay-spectra  time  history,  extraction  of  delay-Doppler-related  channel  parameters,  and 
tap  statistical  analyses.  Common  to  the  four  latter  categories  is  in  the  processing  of  the  SACP  receiver 
and  analog  tape  status  parameters.  The  tape  parameters  (parity  error  count,  frame  lock  loss  count, 

» data  track  being  processed,  polarization  of  data  track)  are  obtained  from  the  PDP  1 1/45  software 

interaction  with  the  FH.  equipment,  whereas  the  SAC  P receiver  parameters  are  embedded  in  the  data 
frame  that  contains  the  correlator  outputs. 

• Through  use  of  the  dual-system  capabilities  of  the  PDP  1 1 '45  ground  station,  the  algorithm 

execution  sequencing  described  below  is  followed: 

a.  The  dual  system  is  used  to  simultaneously  quick-look  process  the  H and  V data  tracks; 
this  includes  generation  of  tape  D3. 

b.  From  the  quick-look  outputs,  periods  of  data  hiterest  and  receiver  configuration  station- 
arity  are  identified. 

c.  The  dual  system  is  used  to  simultaneously  reformat  the  H and  V data;  each  reformatted 
data  stream  is  recorded  on  its  own  separate  digital  tape.  Dl . 

d.  Tape  D3  is  hatch  processed  for  the  three-dimensional  delay-spectra  time  history. 

e.  The  reformatted  tapes  (1)1  ) are  batch  processed  for  the  detailed  data  analysis. 

Although  this  execution  sequence  is  considered  to  be  optimum  for  processing  the  simultane- 
ously received  horizontal  and  vertical  polarization  data,  it  does  not  have  to  be  strictly  followed.  For 
example,  one  may  simultaneously  (luick-look  process  and  reformat  a single  polarization  data  stream  or 
reformat  one  track  while  i|uick-look  processing  the  other.  Similarly,  any  of  the  processing  How  seg- 
ments could  be  eliminated  without  affecting  the  other  operations. 
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ANALMIC  AL  OBJl  C TIVhSOI  ALC'-ORUUM  PROCtSSING  BLOCKS 


A brief  ileseriplion  of  the  ;malytieal  ohjeetive  of  eaeli  of  tlie  major  processing  algorithms  is 
given  below. 

-C2.1  Quick-Look  Real-Time  Playback  Data  Analysis 

This  output  is  directly  available  from  the  PDP  1 1 '45  system  and  provides  both  oscilloscope 
display  plots  and  hard-copy  numerical  output.  Trom  this  the  operator  may  investigate  receiver  param- 
eter configuration,  tap  frequency  spread,  delay  power  spectral  density,  and  analog  magnetic  tape 
status  descriptors.  Primary  use  of  these  data  is  to  identify  time  intervals  that  provide  both  steady-state 
receiver  conditions  and  multipath  scatter  phenomena  of  particular  importance. 

Reformatted  SACP  Digital  Tapes 

Source  analog  tapes  are  processed  to  provide  computer-compatible  digital  tapes  that  represent 
the  comple.x  tap  voltage  time-domain  data  of  the  SACP  receiver.  Typically,  those  periods  of  data 
identified  by  the  quick-look  analysis  as  being  of  particular  interest  were  reformatted  for  the  CONUS 
tests.  For  the  oceanic  tests,  the  bulk  of  the  valid  data  was  converted  to  digital  format. 

3. 2. .5  Delay-Spectra  Time  History 

For  all  periods  of  valid  data  collection,  the  scatter  channel's  delay  power  spectral  density  (psd) 
is  determined  in  a time-running  nonoverlapping  manner  with  psd  estimates  being  calculated  over  a 
2-sec  interval.  The  outputs  that  occur  once  every  2 sec  are  given  in  both  numerical  and  three- 
dimensional  plotted  formats.  Respectively,  these  data  provide  both  a quantitative  and  a comprehensive 
overview  description  of  the  channel’s  time-variant  delay-spectra  characteristics.  This  analysis  is  of 
particular  importance  for  the  CONUS  scatter,  where  terrain  roughness  and  electrical  characteristics 
vary  rapidly  with  distance.  One  use  of  the  3-D  overview  plot  is  for  isolation  of  time  and  tap  bank 
intervals  that  possess  data  for  which  it  is  desirable  to  either  ( I ) reference  the  numerical  output  to 
obtain  quantitative  delay-spectra  information  or  (2)  subject  the  data  input  string  to  detailed  ilelay- 
Doppler  psd  computer  processing. 

3.2.4  Delay-Doppler  Scatter  Power  Spectral  Density.  S(r.co) 

This  function  represents  the  distribution  of  diffusely  scattered  power  arriving  at  the  receiver 
with  Doppler  frequency  w and  time  delay  t.  For  the  zero-mean  complex  Ciaussian  random-scatter 
proce.sS,  Str.w)  completely  characterizes  the  channel  -.tatistics.  For  each  test  condition,  including 
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proho>  >.oiuliult\l  (uor  iIk‘  multip.ilh  ;iiul  opcnitiotial  antenna  systems,  at  least  one  detailed  Str.cjl 
Uinet\on  venerated  Hie  output  is  iiiven  in  numerical  and  three-dimensional  plotted  form.  Further- 
more. since  this  parameter  is  ot  such  fundamental  importance  for  the  interpretation,  modelinj:.  and 
application  ol  the  scatter  channel  phenomena,  S(r.w)  is  also  preserved  on  mattnetic  tape.  This 
provides  a convenient  and  compact  basis  Ironi  which  lulure  analysis  may  be  conducted  without  going 
through  the  tmie-coiisummg  and  expensive  computer  processing  steps  required  to  obtain  the  delay- 
Doppler  psd  from  the  retdrniatted  S.\('P  iligital  tapes. 

-TT.v  Integral  and  Fourier  Operations  on  Sir.u;) 

The  scatter  function  Sir.uil  contains  all  the  ingredients  needed  for  derivation  ofequiva- 
, lent  and  lower  order  channel  parameters.  Software  modules  are  contained  within  the  detailed 

‘ ^ analysis  routines  to  derive  the  following  parameters. 

a.  Joint  tmie-lrequencv  autocorrelation  function. 

RtJ.m  = // Slr.cc)e‘'“’"^^‘^’drda;  (3-1) 

b.  Time  autocorrelation  function, 

R(f.O)  = //S(r.u;)dTe'^‘^du)  (3-2) 

c.  Frequency  autocorrelation  function. 

RtO.m  = //str.widcce'^^^dr  (3-3) 


d.  Doppler  spectrum. 
D(  w I = j StT.cc  I dr 


(3-4) 


e.  Delay  spectrum. 

0(rl  = ^ S(T.cc)  du) 


(3-.S) 


f.  Total  rnis  scattered  energy. 


I! 


Sir.cj)  drdto  . 


(3-6) 


The  time-frequency  autocorrelation  function  and  its  respective  axial  cuts  measure  the 
degree  of  correlation  between  two  signals  delayed  in  lime  by  ^ and  offset  in  frequency  by  . 
These  functions  are  complex  and  therefore  are  formatted  in  terms  of  their  amplitudes  and  phases 
in  the  software  output. 
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3.2.6  Channel  Spread  Parameters 


f-rom  the  D(cj),  Ott-),  K(J.O),  and  R(0.S2)  distributions,  lower  echelon  I'irst-order 
channel  parameters  such  as  the  Doppler  spread,  delay  spread,  decorrelation  time,  and  coherence 
bandwidth  ol'  the  scatter  channel  are  easily  estimateil. 


3.2.7  Noise  Determination  and  Removal  (NDatuIR)  i 

The  outputs  of  the  SACP  multipath  correlator  contain  desired  signal  data,  spurious  signal  i 

terms,  low-pass  additive  thermal  noise,  and  receiver  arithmetic  noise.  Under  normal  SACP  oper-  ^ 

ating  conditions,  the  NDandR  algorithm  statistically  eliminates  these  noise  terms  from  the  data 
on  a tap-by-tap  basis.  The  noise-free  estimate  of  the  delay-Doppler  function  may  then  be  nor- 

mali?ed  for  the  low-pass  filter  attenuation  and  subjected  to  the  integral.  Fourier  transform,  and  i 

spread  determination  steps  as  previously  outlined. 

3.2.8  Antenna  Pattern  liffects  Removal 
, ] 

Design  criteria  for  the  forward  multipath  antenna  were  established  to  provide  a radiation 
pattern  with  nearly  uniform  coverage  over  the  effective  scatter  region.  However,  for  certain 
night  direction  headings,  the  fidelity  of  the  channel  measurement  may  be  enhanced  by  applying 
the  antenna  effects  removal  algorithm.  This  routine  operates  on  the  noise-free  Str.w)  estimate 
to  provide  an  equivalent  scatter  function  that  would  be  measured  with  a uniform  gain  antenna. 

The  multipath  process  performs  a 2-into-l  mapping  operation  wherein  two  surface  returns  are 
mapped  into  one  delay-Doppler  point.  The  application  of  this  algorithm  is  thus  restricted  to  the 
in-plane  flight  goemetry  cases  where  the  cosymmetry  of  the  delay  contours,  Doppler  contours, 

and  scatter  cross-section  (isotropic  surface)  is  exploited.  i 

3.2.*)  Tap  Amplitude  and  Phase  Distributions 

Determination  of  the  fundamental  statistical  properties  of  the  scatter  process  requires 
that  the  received  signal’s  time-vlomain  lluctuations  be  characterized.  The  probability  distribution 
and  the  associated  means  and  variances  ol  these  lluctuations  are  derived  for  each  tap's  I compo- 
nent, O component,  an.'  phase  angle  Hie  processes'  composite  signal,  which  is  constructed  by 
vectorially  summing  the  delay  tap  outputs,  is  also  subjected  to  these  operations.  The  experimen- 
tal distributions  are  compared  with  prediction  for  the  complex  (iaussian  channel.  These  corrcia- 
tions  are  implemented  via  the  x"  goodness-of-fit  test. 
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.V2.10  Tup  Process  Bank  Cross-Correlations 


The  degree  of  coherency  between  any  two  taps  in  a particular  tap  bunk  or  in  cross- 
polari/ed  banks  (i.e.,  one  is  horizontal,  the  other  vertical)  is  measured  through  use  of  the  normal- 
ized cross-correlation  function: 


I 

I 


R 


<X(t)Y*(t-£)> 

•v/<XmX*(t)xY(t)Y*(t)> 


(3-7) 


where: 

X = complex  tap  process  in  either  bank 
Y = complex  tap  process  in  either  bank 
J = time-lag  variable. 

3.2.11  Tap  I and  Q Dependency 


For  random  rough  surface  scattering  where  the  electromagnetic  wave  undergoes  deep 
phase  modulation  at  the  multipath  interface,  we  expect  that  the  I and  Q components  of  the 
received  signals  are  statistically  independent.  This  condition  is  explored  by  determing  the  zero- 
lag  normalized  correlation  coefficient  between  a tap's  orthogonal  components,  i.e.: 


Rl  q(0) 


<10> 


(3-R) 


3.2.12  Tap-Gain  Autocorrelation  Function 

An  estimate'  of  the  channel's  tap-gain  autocorrelation  function.  U(t.J),  is  derived  through 
application  of  the  following  operation: 

U(r.^)  = <X^(t)X};(t-i)>  . (3-^)) 

where: 

r = delay  tup  value 

X^  = complex  process  of  tap  t in  bank  X (i.e..  horizontal  or  vertical) 

I = time-lag  variable  . 

I'his  function  measures  the  autocorrelation  function  of  the  multipath  process  on  a tap-by-tap 
basis  and  is  available  (magnitude)  as  a three-dimensional  output  plot  from  the  software  package. 

'An  alternate  derivation  of  this  function  may  be  obtained  by  inverse  Fourier  transformation  of 
the  S(r.co)  function  with  respect  to  the  co  variable. 
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Two  progiam  options  t‘\ist  in  this  aluoritlim:  ( 1 ) oross-polarized  estimate  of  UtT.JI  and 
(2)  noise  estimate  iemo\al. 

Hie  eross-polan/ed  estimate  of  U(7-,$)  is  ealeiilated  by  replacing  the  conjugated  variable  in 
equation  with  the  appropriate  tap  of  the  bank  containing  the  orthogonal  polarized  return. 

.Assuming  the  independence  of  multipath  and  noise,  we  arrive  at  a noise-free  estimate  of  the 
tap-gain  autocorrelation  function 


U (r.^)  = <X  (t)\*(t-?  )'>-!£  <.\.(T)Xf(r-{) 

N i = ■ 


(3-10) 


1=1 


where  .Xj  is  the  complex  output  of  a multipath-free  tap  (usually  taken  from  a region  of  the  bank 
preceding  the  specular  point  return). 


3.2.13  System  Calibration  Parameter  Data 


Magnetically  recorded  data  pertaining  to  receiver  system  operation  (i.e.,  direct  and  multipath 
channel  gains,  etc.),  transmitter  power  amplifier  outputs,  and  aircraft  (light  parameter  descriptors 
are  computer  reduced  to  aid  in  the  normalization  of  the  scatter  channel  power  returns.  These 
data  also  serve  as  a data  collection  integrity  measure  and  are  used  primarily  to  augment  the  logged 
(light  test  data. 


3.3  MULTIPATH  DATA  BASE  COMPONENTS 

The  data  base  delivered  to  DOT/TSC  relative  to  the  multipath  channel  characterization  test 
consists  of  the  following  items; 

a.  Library  of  analog  SACP  receiver  source  data  tapes,  l liese  data  were  recorded  at  the 
NAS  ground  station  on  dbOO-  or  7200-ft  reels  of  I '2-m.  magnetic  instrumentation  tape 
using  .Ampex  LR-ldOO  and/or  ER-200n  recorders.  Ihe  recorded  data  corresponds  to  the 
unprocessed  data  obtained  from  the  outputs  of  the  SA(  P receiver  during  test. 

b.  Library  of  reformatted  digital  data  (apes  (Dl  ).  Lhesc  loOO-bpi.  nine-track  digital  tapes 
contain  the  SACP  receiver  output  data  in  computer-compatible  digital  formal  Ilie  DI 
tapes  are  generated  by  the  Boeing  PDP  I I M.S  facility  and  have  a formal  as  described  in 
section  4.22. 1 . 
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c.  Library  of  delay-spectra  time  history  tapes  (D3).  These  IbOO-bpi.  nine-track  digital  tapes 
contain  the  delay-spectra  time  history  arrays  computed  by  the  dual  I’HF*  I I '45  computer 
system.  In  most  cases,  the  files  containing  these  data  are  physically  located  on  the  same 
tape  as  the  corresponding  reformatted  digital  data  described  as  1)1  (above),  lape  format 
is  described  in  section  4.22.3. 

d.  Library  ot  Sj^tr.cc)  save  tapes  (1)2).  Tape  1)2  contains  the  computed  noise-present  esti- 
mate of  the  channel’s  delay-Doppler  scatter  function.  The  tape  fonnat  is  described  in 
section  4.22.2. 

e.  Library  of  airborne  system  parameters  source  tapes.  I'hese  8t)t)-bpi.  seven-track  computer- 
compatible  digital  tapes  were  recorded  onboard  the  KC-135  aircralt  terminal  using  a 
Kennedy  8707  recorder.  The  tapes  contain  power  level  calibration  data  and  other  infor- 
mation relevant  to  the  test  parameters  during  test.  Infomiation  content  and  tape  format 
are  described  in  section  4.23. 

I.  Computer  programs  (card  decks)  and  listings  corresponding  to  the  programs  used  for  de- 
tailed analysis  of  the  reformatted  multipath  data  (Dl  ) tapes.  Description  of  algorithms 
are  given  in  sections  3 and  4. 

g.  C omputer  programs  (card  decks)  and  listings  corresponding  to  the  physical  optics  vector 
scatter  model  described  in  section  5. 

h.  Punch  paper  tapes  providing  radiation  distribution  plots  of  antenna  range  data  for  the 
front  multipath  antenna  for  various  polarizations  and  pointing  angles. 

The  following  additional  data  base  resulted  from  a subcontract  to  CNR,  Inc. 

a.  Retormatted  800-bpi,  seven-track  multipath  computer-compatible  digital  data  tapes 
generated  from  selected  Boeing  PDP  I 1/45.  I(i00-bpi.  nine-track  1)1  tapes. 

b.  Transposed  multipath  digital  data  tapes  listed  in  volume  V',  section  5.4.(i. 

c.  Sea-state  buoy  800-bpi,  seven-track  digital  data  tapes  containing  measured  buoy  variables 
stripped  from  the  analog  buoy  data  source  tapes.  ( Lapes  correspond  to  data  aaiuired  on 
January  2‘)  and  30.  I‘)75.  plus  March  25  through  .April  2.  1475.)  Lhese  PCM  data  strip- 
ping operations  employed  a capability  existing  at  DOT  I SC. 

d.  Backup  program  tape  containing  the  programs  used  for  sea-stale  buoy  data  reduction. 
tl-V  correlation  and  statistical  analysis,  and  transpose  and  scattering  function  DRandA 
These  programs  are  discussed  in  sections  4.2  I ami  appendix  B.  Programs  can  be  run  on 
the  DOr/TSC  PDP  10  in  time-share  or  batch  processing  modes. 
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4 Ml  LTIPATH  AKXJRITHM  DESCRIPTION 


Ihix  M-aion  p^null•^  .1  more  lU-lailod  lU'scnptioii  ol  the  multipath  algorithms  presented  in  fig- 
ures .'-I  anil  .'-2 


4 1 II  I 1 Ml  IK't  ini  A IM’L'l 

Data  proeessmg  test  eonditions  on  the  PDP  1 1 '45  are  detuu'd  hy  ( I ) program  type  and  option, 
( 2>  data  traek  or  traeks  to  he  analyzed,  and  (3)  the  start  and  stop  time  of  the  interval  to  be  processed. 
Telemetry  data  input  to  the  PDP  II  45  is  preceded  hy  transferring  the  EMR  telemetry  front  end  con- 
trol to  the  tape  search  unit,  which  causes  the  analog  tape's  time  code  track  to  be  searched  for  the 
specified  interval  start  lime.  When  the  tape  search  unit  enters  the  interval,  the  digitized  PCM  data  is 
routed  to  the  PDP  I 1 45  memory  in  blocks  of  four  frames  (each  frame  word  is  input  as  a 16-bit  entry: 
the  most  significant  seven  bits  are  always  0 while  the  least  significant  ^ hits  contain  the  data). 

The  Input  data  is  double  buffered.  When  the  input  of  a block  of  data  is  completed,  an  inter- 
rupt occurs  and  causes  two  events:  ( 1 ) new  data  starts  entering  the  second  buffer  and  (2)  the  process 
software  begins  to  operate  on  the  buffer  that  has  just  been  filled.  If  a TF'E  frame  lock  loss  event 
occurs,  the  processor  is  notitled  hy  another  interrupt  that  enables  software  to  maintain  status  logic. 
For  both  the  reformatting  and  nuick-look  analysis'delay-spectra  time  history  program  operations,  data 
input  is  typically  run  at  a rate  equal  to  roughly  80,000  words, 'sec.  ,At  a sampling  rate  of  300  frames/ 
sec.  this  corresponds  to  a 2-to-l  processing  slowdown  factor. 


4.2  DATA  DhCODINX..  I RROR  RECOVERY,  AND  REFORMATTING 

To  return  the  data  to  the  form  as  output  from  the  SACP  receiver,  it  must  be  parity  error 
checked,  stripped  of  its  parity  hit.  and  odd-bit  complemented.  These  three  steps  are  accomplished  in 
one  operation,  which  consists  of  using  the  data  as  an  index  to  a table  of  reformatted  data  values  and 
parity  error  tlags.  Parity  errors  are  corrected  by  data  replacement  from  the  corresponding  data  in  the 
preceding  frame.  An  accumulative  count  of  detected  parity  errors  is  maintained. 

When  the  reformatted  data  is  written  on  magivtie  tape,  it  is  double-buffer  stored  in  vectors  of 
length  equal  to  I 2 fr.imes.  Tliese  output  buffers  also  contain  data  pertaining  to  analog  magnetic  tape 
traek.  Irame  lock  loss  count,  parity  error  count  and  directory,  and  polarization.  The  digital  tape  for- 
mat IS  described  in  section  4.22. 
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4 3 QUICK-l.CX)K  DATA  ANALYSIS 

Prior  to  quick-look  analysis  of  the  SAC'P  data  tapes,  we  must  select  the  time  interval  to  be 
processed  and  establish  the  subinterval  time  length  over  which  sample  space  construction  is  periormed. 
lypically.  the  processing  interval  consists  ot  the  entire  llight  test  condition  time  segment  over  which 
log-look  records  indicate  valid  data  was  obtained.  VVe  also  typically  select  a subinterval  length  of  2 sec 
because  this  appears  to  give  adequate  resolution  lor  the  CONUS  data  runs,  which  encounter  surlaces 
whose  electrical  and  physical  characteristics  change  rapidly. 

bor  each  subinterval  sample  space,  the  quick-look  process  algorithm  perfomis  the  tbllowing 
operations:  ( I ) extracts  the  receiver  and  analog  tape  status  parameters  associated  with  the  first  and 
last  frame  of  the  subinterval.  (2)  estimates  the  multipath  tap's  frequency  spread,  (3)  estimates  the 
direct  and  multipath  tap's  delay  spectra,  and  (4)  writes  the  above  as  one  record  on  digital  tape.  The 
quick-look  analysis  results  are  not  intended  to  provide  definitive  quantitative  channel  parameter 
results.  Their  primary  use  is  to  direct  attention  toward  test  intervals  that  warrant  further  detailed 
analysis. 

An  algorithmic  flow  diagram  for  the  above  operations  is  shown  in  figure  4-1 . By  intercepting 
the  receiver  and  analog  tape  status  parameters  at  the  beginning  and  end  of  each  subinterval,  a measure 
of  receiver  parameter  stationarity  is  provided.  These  data  occupy  the  first  24  bytes  of  each  data  set 
(a  data  set  consists  of  the  SACP  data  frame  preceded  by  the  analog  tape  system  parameters),  as 
described  in  section  4.22. 

The  indirect  tap  bank's  frequency-spread  estimate  is  derived  by  counting  the  number  of  times 
a particular  tap's  I component  crosses  the  zero  voltage  level.  This  is  implemented  in  the  PDP  I I 45 
by  maintaining  a counter  array  and  a previous  sign  array  in  memory.  Each  element  of  the  two  arrays 
corresponds  to  a particular  SAC'P  tap  output.  The  previous  sign  array  contains  sign  information 
(positive  or  negative)  for  the  last  frame  processed.  This  information  is  compared  with  the  sign  of  the 
current  array  element;  if  they  differ,  the  corresponding  element  of  the  counter  array  is  incremented 
and  the  sign  of  the  element  in  the  previous  sign  array  is  complemented.  The.se  operations  take  place 
on  each  frame  that  is  input  to  the  quick-look  program.  .At  the  end  of  a subinterval  sample  space,  the 
counter  array  is  written  on  magnetic  tape  without  normalization  and  is  subsequently  zeroed  for  reuse. 

For  the  case  of  a signal  possessing  complex  Gaussian  amplitude  statistics,  the  zero-crossing 
technique  yields  a good  approximation  to  its  frequency  spread  as  defined  in  reference  4-1 . Because 
the  L-band  multipath  signal  is  known  to  have  Rayleigh  amplitmle  statistics,  the  frequency-spread  data 
is  used  as  a measure  of  each  tap's  Doppler  spread.  T his  is  valid  for  the  taps  that  possess  high  signal-ti)- 
noise  ratio  but  is  invalid  for  poor  signal-to-noise  ratios  since  the  zero-crossing  counter  is  activated 
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pnin.inly  hy  the  noise  lluelii.ilion.  I or  e\;implc.  the  I eomponent  ol  .i  parlieular  delay  tap  may  be 
presented  as 


DItt » = n|(t  t + mp|(l»  . (4-1  ) 

where. 

1)1(1 1 = I eomponent  ol  a partienlar  ilela\  tap 

iijiti  = I component  ol  the  noise  content  ol  the  tap 

mp|ili  = I component  ol  the  imillip.ilh  content  of  the  tap 

hot  the  maiorit\  ol  dela\  laps,  n|ll»  Is  predominantly  thermal  and  wonkl  he  expected  to  possess 
the  complex  (iaiissiaii  sialelKs  and  to  have  a /ero-erossing  rate  that  is  roughly  erpiivalent  to  the  band- 
width ot  the  low -pass  niter  opi  ration  The  combination  of  the  two  processes  thus  always  gives  an 

upper  limit  measure  to  the  lap  s l)oppler  spread.  In  the  region  ot  the  postspeeular  point  return,  the 
.ilgorilhm  output  vields  a relatively  good  estimate  ot  the  true  Doppler  spread.  Kor  prespecular  point 
c.ipiure  and  captuie  tar  delayed  Irom  the  specular  return,  the  freiiuency-spread  measure  approaches 
the  noise  haitdwulth  These  char.icicrisiics  are  illustrated  in  the  output  data  described  in  section  4.24. 

rile  ipiKk-hiols  delay -spectra  calculations  are  b.ised  on  a sample  space  constructed  Irom  the 
lust  I'r.ime  ol  each  butler  input  li  e.,  one  out  ot  every  lour  Irames).  This  was  implemented  to  speed 
up  the  algorithm  oper.ilion  .md  at  the  same  lime  improve  the  independency  of  the  sample  space  ele- 
ments ( ly  pically  the  multipath  dal.i  is  overs.impled ).  f or  each  lap.  including  both  the  direct  and 
multipath  banks,  the  process  algonlhin  proceeiK  as  follows: 

a.  Hie  squares  of  the  I and  Q components  for  the  particular  tap  are  determined  by  using  the 
v.ihie  as  an  index  into  a table  of  squares.  This  operation  is  considerably  faster  than 
multiply  ing  the  number  by  itself. 

b.  Hie  squares  of  the  I and  Q components  are  added  together. 

c.  Ihe  1"  + Q‘  sum  IS  accumulated  from  frame  to  Iraine  by  adding  it  to  a memory  location 
corresponding  to  that  tap.  A ,^2-bit  result  is  maintained. 

d \t  the  end  of  a submterval  condition.  Ihe  delay-spectra  array  is  written  on  magnetic  tape 

I ike  the  frequency-spread  estimate.  Ihe  delay-spectra  me.isure  is  contaminated  by  receiver 
noise.  In  this  case,  however,  the  elTect  is  not  nearly  so  significant.  I he  receiver  noise  may  be  esiimatcvl 
by  examining  Ihe  process  output  for  laps  prior  to  the  specular  point  return  to  derive  a mean  estimate  of 
the  noise  content  of  a multipath-free  l.ip  This  value  can  be  subtracted  from  the  entire  induced  tap 
bank.  It  should  be  noted,  however,  that  on  occasion  certain  laps  are  possessed  by  an  anomalous  noise 
content  (e.g..  taps  m the  vicinity  of  1,17,  2(>.  HI . aiul  1 1 I have  been  observed  to  tall  into  this  category  I 
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4.4  ni  L.A^-SPK  1 R.A  1 IMP  HISTORN 

Hie  ilela\ -spectra  time  history  algorithm  operates  on  the  data  records  created  during  the 
quick-look  analysis  and  preserved  on  the  D3  storage  tape.  A quick-look  analysis  data  record  is  gener- 
ated once  every  AS  sec  and  contains  ini'ormation  on  ( I ) the  receiver  and  analog  tape  status  param- 
eters at  the  heginning  and  end  ot  the  subinterval.  (2  I unnornialized  frequency-spread  estimates  for  the 
indirect  tap  bank,  and  t.H  imnormalized  delay-spectra  estimates  for  both  the  direct  and  indirect  tap 
banks. 


[-igure  4-2  illustrates  the  processing  steps  associated  with  this  algorithm.  Analysis  of  the 
receiver  and  analog  tape  status  parameters  is  covered  separately  in  section  4.10.  The  frequency-spread 
and  delay-spectra  estimates  are  normalized  as  follows: 


FSliiT-.tjl 


(4-21 


Q(nr.tj)  = (4'M)Oj(nr) 


(4-.^) 


where: 


b'S(  iiT.tj 

FS|(  iiTt 


OjOlTl 

n 

7 

'i 

AS 

'start 

,\1 

SR 


I = normalized  frequency-spread  estimate  of  the  n'*’  tap  for  time  interval  tj 


= unnormalized  frequency-spread  estimate  of  the  n'”  tap  as  contained  on  the  i**'  record 

th 


Qtnr.tj)  = normalized  mean  square  energy  of  the  n"'  tap  for  time  interval  tj 


imnormalized  mean  square  energy  of  the  n"'  tap  as  contained  on  the  i"’  record 
tap  number 

chip  width  of  PN  sequence 
'start  ^ 

time  length  of  subinterval  (typically  2 sec) 
start  time  of  the  first  record  processed 
SR  .\  AS 

frame  rate  (sec”'  ) 


The  FS(nr.tj)  array  is  presented  numerically  for  all  I 12  taps  Similarly,  the  numerical  array 
for  0(  nr.tj ) is  given  for  all  six  direct  taps  and  all  1 I 2 indirect  taps.  Any  contiguous  set  of  taps 
including  both  direct  and  indirect  may  be  selected  for  the  three-dimensional  plotted  delay-spectra  time 
history  array.  This  plotted  output  format  gives  a concise  visual  representation  of  the  channel’s  time- 
variant  delay  spectra. 
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Figure  4-2.  Delay-Spectra  Time  History  Algorithm 


Incoherent  power  siinimations  of  the  direct  and  multipath  total  energies  are  also  cojitained  in 
this  algorithm  and  are  calculated  as  follows: 


6 

PD(tj)  = 4/M  Y,  Oj(njT)  , (44) 


Pl(tj)  = 4/M  2 Oj(n|T)  , (4-5) 

"1=1 

where: 

PD(tj)  = incoherent  power  summation  estimate  of  the  direct  tap  energy 
Pl(tj)  = incoherent  power  summation  estimate  of  the  indirect  tap  energy 

Note  that  this  routine  operates  on  the  noise-contaminated  frequency-spread  and  delay-spectra 
estimates.  Thus,  the  comments  in  section  4.3  relative  to  the  effects  of  receiver  noise  also  apply  to  this 
algorithm. 
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4.5  TAP-TIME  MATRIX  TRANSPOSITION 

The  tap-time  matrix  transposition  operation  consists  of  reordering  the  input  data  sequence. 
Although  this  operation  does  not  affect  the  data  analysis  computations,  it  is  mentioned  here  for  two 
reasons:  ( I ) to  clarify  the  presentation  of  the  algorithm  implementation  and  ( 2)  to  indicate  the  level 
of  coding  and  computer  time  that  it  requires.  For  the  software  discussed  in  this  document,  the  output 
of  the  matrix  transposition  is  used  primarily  to  compute  the  delay-Doppler  scatter  tunction.  However, 
it  was  implemented  as  an  independent  computer  program;  thus  it  may  he  readily  used  for  additional 
data  reduction  routines  such  as  tap-amplitude  statistics,  etc. 

The  transposition  operation  consists  of  reordering  the  data  from  a time  sequence  to  a tap 
sequence.  The  input  data  may  be  viewed  as  a two-dimensional  matrix.  A column  of  the  matrix  repre- 
sents the  time  sequence  of  one  component  of  one  tap.  A row  of  the  matrix  corresponds  to  a simulta- 
neous sample  of  the  real  or  imaginary  parts  of  a particular  tap  in  the  hank.  The  data  is  input  to  the 
computer  on  a row-by-row  basis  (i.e„  time-ordered).  However,  for  efficient  data  analysis,  the  matrix 
must  be  accessed  by  column  or  tap.  The  matrix  si/e  presents  a problem:  e.g..  for  a 5 I 2-poinl  complex 
FFT  over  the  total  direct  and  indirect  tap  banks,  we  must  work  with  a data  array  on  the  onlerof 
I 20.000  values  of  the  SACP  receiver’s  eight-bit  words. 
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Memory  and  disk  space  with  random  access  are  used  to  perform  the  transposition  operation. 
For  practical  reasons,  two  other  operations  are  performed  at  the  same  time  the  input  matrix  is  being 
transposed;  1 1 ) data  validity  checks  and  (2»  tap  reordering  (from  the  recorded  interleaved  order  to  a 
monotonically  increasing  order). 

The  program  high-level  algonthm  (including  the  coincident  operations)  is  given  in  figure  4-3. 
Program  input  consists  of  reformatted  data  tape,  number  of  frames  per  segment,  number  of  segment 
tape  selected,  and  20  characters  for  case  identification.  Program  outputs  consist  of  the  file  with  tap- 
ordered  data  and  a history  and  summary  of  validity  check  messages. 

The  transposition  operation,  shown  in  figure  4-4.  consists  ol  the  tollowing  basic  steps; 

a.  Data  is  input  to  the  CDC  6600  core  memory  frame  by  frame. 

b.  For  each  frame,  the  tap  data  words  (1.0  Point  pairs)  are  repositioned  in  core  memory  to 
represent  a subset  of  the  overall  transposed  matrix  as  diagrammed  in  figure  4-5.  In  this 
operation  the  number  of  frames  processed  is  sized  to  occupy  a core  space  of  50K  octal. 

c.  W hen  the  allotted  core  space  is  filled  or  when  the  last  frame  of  a data  interval  segment  has 
been  processed,  the  transposed  subsegment  is  transferred  from  core  to  the  random  access 
disk  file,  tact)  tap  segment  is  a.ssigned  a unique  record  number. 

d.  After  completing  the  frame  loop,  the  random  access  file  is  acces.sed  and  the  tap  records 
are  rearranged  to  yield  a sequential  Tde  representation  of  the  tap-time  matrix. 

4.6  DATA  VALIDITN  ( HtCK 

The  data  validity  check  algorithm  provides  information  pertaining  to  time  (merged  IRIG-A) 
and  subframe  index  counter'  between  adjacent  frames  that  arc  input  to  the  detailed  data  analysis 
routine.  For  proper  operation  one  expects  that;  ( 1 ) between  adjacent  Irames  the  time  interval  will 
be  constant  and  equal  to  the  reciprocal  of  the  sampling  rate  and  (2)  the  subframe  counter  advances  by 
one  count.  These  criteria  are  checked  for  each  frame  of  data  input  to  the  program.  If  an  error  occurs, 
(lags  are  set  without  interrupting  the  program  execution.  .At  the  end  ot  a data  segment  input,  an  error 
summary  is  provided  to  delineate; 

Total  number  of  sidsframe  error  counts 
Total  number  ol  lime  sequence  errors 

Mean  and  standard  deviation  of  lime  intervals  between  frames 

' The  SAC  Pdata  frame  contains  a suhfrai  ie  counter  that  runs  f rom  0 to  5 m a continuous  cyclic 
manner  by  advancing  one  count  per  frame 
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SUBSEGMENT  INPUT  STRING  OF  K FRAMES 


Figure  4-5.  Subinterval  Transposition  Structure 
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T.\P  ( ONTR.UTION 


Tap  contraction  is  a software  option,  liesigned  primarily  to  save  computer  time,  which  is  avail- 
able to  reduce  the  effective  size  of  the  indirect  tap  bank  array.  This  may  take  on  any  or  all  of  the 
followiiii:  operations: 

a.  Eliminate  taps  from  the  beginning  of  the  bank. 

b.  I liininate  taps  from  the  end  of  the  bank. 

c.  Coherently  sum  any  number  of  adjacent  taps  to  yielil  a new  array  that  has  a decreased 
tap  resolution  width  (i.e..  one  may  choose  to  sum  over  taps  to  procure  an  array  of  16 
taps  as  opposed  to  the  original  set  of  112). 

The  first  two  operations  are  implemented  in  the  matrix  transposition  algorithm  to  eliminate 
taps  that  are  known  to  contain  no  significant  amount  of  multipath  signal.  The  coherent  tap  summa- 
tion option,  if  implemented,  is  performed  immediately  prior  to  inputting  the  data  to  the  scatter 
function  algorithm.  Since  the  resultant  decrease  in  resolution  is  undesirable,  this  latter  option, 
although  available,  is  normally  not  used. 


4.8  DlRhCT  TAP  PROCESSING 

In  the  interpretation  of  the  multipath  signal  strength  data,  it  is  useful  to  reference  the  scat- 
tered energy  with  respect  to  the  energy  arriving  over  the  direct  line-of-sight  path.  The  total  energy 
contained  in  the  direct  signal  is  estimated  by  coherently  summing  the  direct  tap  complex  voltage  out- 
puts to  form  a composite  direct  signal.  The  mean  square  value  of  this  signal  is  then  determined  over 
the  time  interval  of  interest,  i.e.. 


where: 

<101- '> 
Olj(i) 
l)Oj(i) 
N 


mean  square  direct  composite  signal 
in-phase  component  of  i**’  sample  of  the  j’*'  direct  tap 
quad-phase  component  of  i'*’  sample  of  the  j**’  direct  tap 
the  number  of  points  in  the  sample  space  . 


(4-(i) 


In  general,  the  direct  tap  processing  algorithm  acts  on  the  same  data  space  that  is  processed  for 
the  S(r.u))  determination  (.see  fig.  4-6).  The  effects  of  receiver  noise  on  the  <II)1‘-'  estimates  are 
assumed  to  be  negligibly  small  since  the  direct  link  typically  operates  at  a G-'Nj,  density  in  excess  of 
60  dB-llz  and  the  processing  bandwidth  is  on  the  order  of  7.5  llz. 


Figure  4 6 Delay  Doppler  Scatter  Function  Algorithm 


4.0  DELAY  SPECTRA  (TIME-DOMAIN  ANALYSIS) 

To  serve  primarily  as  a cheek  on  the  integrity  of  the  operations  used  to  derive  S(r,w)  we 
perform  a time-domain  analysis  on  the  complex  tap  outputs  to  derive  an  estimate  of  the  channel’s 
delay  spectra.  Tliis  estimate  may  then  be  compared  with  the  spectra  obtained  by  integrating  S(t.co) 
(noise  not  removed)  over  the  Doppler  variable. 

Inputs  to  this  routine  consist  of  the  identical  complex  tap  voltage  points  utilized  in  the 
ensemble  array  over  which  S(t,cj)  is  calculated  (see  fig.  4-6  for  processing  sequence).  These  points 
are  operated  on  as  follows  to  derive  the  delay-spectra  estimate: 


Otda<"’->  = I/nE  + , 


Otda*’’’’^  = mean  square  energy  arriving  at  receiver  with  delay  nr 

n = tap  number 

T = chip  width  of  transmitted  PN  sequence 

l^^(i)  = in-phase  component  of  n^’’  tap’s  i**’  ensemble  array  value 

C?„^(i)  = quad-phase  component  of  n**’  tap’s  i‘*^  ensemble  array  value 

N = total  number  of  points  in  the  array  ensemble  . 


4.10  RECEIVER  AND  ANALOG  TAPE  STATUS  PARAMETERS 


The  receiver  and  analog  tape  status  reduction  algorithm  operates  on  the  first  I hits  of  a data 
record  contained  either  on  the  D1  or  D2  digital  tape.  With  the  exception  of  one  minor  mathematical 
operation,  the  algorithm  process  consists  of  decoding,  labeling,  and  printing  the  status  parameter  bits. 
The  mathematical  operation  entails  determining  the  ratio  of  the  SACP's  direct  to  indirect  receiver  tap 
bank  gains.  Prior  to  proce.ssing.  the  bit  pattern  has  the  format  described  in  section  4.22.  Decoding  of 
the  frame  counter,  lock  loss  frequency  counter,  and  parity  error  counter  is  straightforward  and  con- 
sists of  converting  the  16-bit  binary  words  into  their  decimal  equivalent.  The  merged  time  code  words 
are  converted  from  this  BCD  format  into  the  hours,  minutes,  seconds,  and  decimal  fractions  of 
seconds  distribution.  Bytes  7 and  H contain  data  track  and  polarization  status;  the  polarization  byte 
is  not  operated  on  whereas  the  data  track  byte  is  converted  to  its  decimal  equivalent.  SACP  receiver 
parameters  occupy  bytes  (eight-bit ) I .f  through  22  plus  byte  24  of  the  data  record.  Lhey  are  decoded 
according  to  tables  4-1  and  4-2. 
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TABLE  4-1.  TAPE  FORMAT  FOR  SACP  RECEIVER  AND  TAPE  STATUS  PARAMETERS 


i 

1 


] 


! 

J 


■i 

Bit 

7 

6 

5 

4 

3 

2 

1 

0 

(LSB) 

B 

CR2 

CRl 

CR0 

wn 

A’ 10 

SCI 

SC0 

TLL 

OMl 

OM0 

PT3 

FT  2 

FTl 

F'7 

■9 

DS4 

DS3 

DS2 

DSl 

DS0 

3W2 

B'Nl 

BW0 

15 

MS4 

MSB 

MS2 

'ISl 

•1S0 

ifl 

— 

13 

TV7 

TVC 

TVS 

TV4 

TV3 

tv: 

TVl 

17 

FV15 

FV14 

FV13 

FV12 

FVll 

FV10 

FV9 

FV3 

20 

0E7 

DE6 

DEB 

DE4 

DEB 

DE2 

NEl 

DE0 

19 

DE9 

0E3 

HR23 

HR10 

HR3 

HR4 

HR2 

MO  j 

22 

EL 

MI40 

•’120 

MI10 

MIS 

MI4 

■112 

■'ll 

21 

— 

SE40 

SE20 

SE10 

SE3 

SE4 

SE2 

SEl 

24 

— 

— 

FV7 

FV6 

FV5 

FV4 

FV3 

f\/2 

.'Jote:  Tiie  entries  in  the  bit  ■natrix  are  defined  in  Table  4-2. 


TABLE  4-2.  LEGEND  FOR  SACP  RECEIVER  AND  TAPE  STATUS  PARAMETERS 


TSP  0=T,  l=T/2 

CR  CHIP  RATE  0 = irHZ,  1 = 5,  2=2.5,  3=1.25,  4=1.0,  5=0.5 

AM  AMTP\:iA  M02E  0=::'JX,  1=NQ>’,M(V),  2 = ‘iRRM(H) 

SC  SYSTEM  C0'Fiai’’ATIO';  0=A/C  X'TT,  1=A/C  RCV  A,  2=A/C  RCV  5 

FLL  FREQUENCY  LOOP  .'ODE,  0=AUTO,  1=MANUAL 

TLL  TIMING  LOOP  'iODE  . 0=/U'TO,  1 = ‘1ANUAL 

on  OPERATING  MODE,  0--TEST.  1=SEARCH,  2 = RUN 

FT  FREQUENCY  LOCP  THRESHOLD,  THR£SH0LD=2**(N+5) , 

IF  N=1001  = 9,  TilEN  THP.ESHOL 0=2**14=16334 
OS  DIRECT  TAP  SCALE  FACTOR,  THIS  IS  A 5 BIT  TWOES  COMPLI'IENT 

INTEGER  'IHICH  CAN  BE  C'^NSIDERED  TO  BE  AN  EXPONENT  FOR  THE 
DIRECT  PATH  TAPS  V.'HE'i  THEY  ARE  TREATED  AS  BLOCK  COMPLEX  NUMBER 
IF  = -3  then  the  SCALE  FACTOR  IS  .125. 

OW  MULTIPATH  TAP  LGX  PASS  FILTER  BANDWIDTH,  0=G00HZ,  1=300HZ. 

2=150HZ,  3=75HZ.  4=37. 5HZ 

MS  '1ULTIPATH  TAP  SCALE  FACTOR,  SA'DIS  INTERPERTATION  AS  DS 

IFL  INTERFACE  CLOCK  PHASE  L''CKED  LOOP,  0=NOT  LOCKED.  1=L0CKE0 
03Q  0=FREQ  LOOP  BELOW  THRESHOLD,  1=A30VE  TH.dEPHOLD  LEVEL 

TV  TIMING  LOOP  D-TC-A  CONVERTER  VOLTAGE 

FV  frequency  loop  D-TO-A  :'''’''ERTl=  V'ILTAGE  ;153ITS  data; 

DE  DIRECT  PAT^'  TO  MULTIPAT-  T'^  SPACING  IN  CHIPS  - DELTA  L - 

(10  BITS) 

EL  EARLY -LATE  SWITCH,  0=T  . I ” 

HR  TDE  CODE,  HOURS 

•II  TIME  CODE.  MINUTES 

S£  TIME  CODE,  SECONDS  , 
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4. 1 1 DELAY  -DOPPLFR  SCATTER  FUNCTION 


The  multipath  channel’s  delay-Doppler  scatter  function  estimate.  S(T.tu).  is  derived  by  calcu- 
lating the  discrete  power  spectral  density  of  the  recorded  SACP  time-domain  tap  outputs.  Figures  4-6 
and  4-7  illustrate  the  fundamental  processing  steps  associated  with  this  algorithm.  Five  operations  are 
involved  in  each  PSD  calculation:  ( 1 ) application  of  a data  window  to  the  complex  time-domain 
input.  C)  FFT  calculation  of  the  complex  input  signal's  discrete  PSD  periodogram.  (3)  spectral 
averaging  of  each  periodogram  output.  (4)  ensemble  averaging  of  a set  of  periodograms.  and  (5) 
normalization. 


4,11.1  Input  Data  Window 

A generalized  raised  cosine-arch  function  is  used  to  weight  the  gated  complex  fast  Fourier 
transform  (FFT)  data  input.  The  purpose  of  this  operation  is  primarily  to  reduce  side  lobes  in  the 
periodogram’s  output  (i.e.,  it  is  assumed  the  FFT  is  operating  on  one  cycle  of  a periodic  time-domain 
process,  and  thus  any  difference  between  the  values  of  the  first  and  last  data  inputs  appears  as  a sharp 
signal  discontinuity;  the  data  window  overcomes  this  by  forcing  the  leading  and  trailing  values  to  be 
nearly  identical).  The  weighting  function.  W(i).  may  be  expressed  as: 


W(i) 


= (l/2)  I - cos  ((i/aN)w) 


i < aN 

aN  < I < N(  1-a)  ' 
i > N(  1-a) 


(4-8) 


where: 

i = the  input  array  point  number 

a - the  traction  ot  points  over  which  the  leading  and  trailing  edges  of  the  window  are  applied 
N = total  number  of  points  in  the  array  . 


Parameter  a is  a program  input:  by  setting  it  to  0.5  we  obtain  the  Hanning  window. 


4.1  I .:  Discrete  PSD 

The  fast  Fourier  transform  algorithm  is  used  to  estimate  the  time-domain's  power  spectral 
density.  The  FFT  is  merely  an  efficient  algorithm  for  performing  the  well-understood  discrete  Fourier 
translorm.  Inputs  to  the  program  are  in  the  torm  ot  complex  tap  voltage  point  pairs.  The  array  size 
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INPUT  : A(jAt)  jtP  SAMPLE  OF  A COMPLEX  TAP  OUTPUT.  GIVEN  EVERY 
At  SECONDS 


TIME-DOMAIN  WINDOW 


DISCRETE  FOURIER  TRANSFORM 


POWER  SPECTRAL  DENSITY 


SPECTRAL  AVERAGING 


A( jAt) 


E(Wi  ) 


OUTPUT  : UNNORMALIZED  PSD  FOR  A COMPLEX  TAP  INPUT 


Figure  4-7.  Tap  Power  Spectral  Density 


7 


may  take  on  any  reasonable  value;  zero-filling  is  used  to  extend  the  array  to  the  nearest  2'^  value.  Out- 
puts from  the  FFT  periodogram  occur  every  I /AT  Hz  (AT  = input  gate  width  in  seconds)  and  are  in 
the  form  of  complex  numbers;  i.e.. 


C(kAf)  = E BUAt)e"-’^j*^^'‘'^‘';  k=l.N  ’ 


i=0 


where: 

C(kAf)  = complex  periodogram  output  at  frequency  kAf 
Af  = 1 /AT 

AT  = NAt 

At  = (sampling  rate)”* 

B(iAt)  = i**’  complex  input  point  as  modified  by  W(i) 

j =v^. 

The  magnitude  squared  of  the  periodogram  output  represents  the  PSD;  i.e.. 


D(kAf)  = C(kAf)C*(kAf)  , 


where  D(kAf)  is  the  power  at  frequency  kAf. 


4. 1 1 .3  Spectral  Averaging 


It  is  well  known  that  the  basic  periodogram's  prediction  uncertainty  does  not  approach  zero  as 
the  number  of  input  points  becomes  large.  One  technique  to  overcome  this  estimation  instability 
involves  convolution  processing  of  the  PSD  with  a spectral  gate  of  specified  bandwidth.  Referred  to 
herein  as  “spectral  averaging,”  this  procedure  is  represented  as; 


F(kAf)  = ^ SG(i)D((k-ti)Af)  . 
i=-M 


where: 

F(f)  = smoothed  PSD  estimate  at  frequency  f 

S(i(i)  = spectral  gate  function 

D(f)  = unsmoothed  PSD  estimate  at  frequency  f 


Program  options  allow  for  specification  of  the  gale's  bandwidth  (M)  and  for  use  of  a rectan- 
gular function. 
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4.11.4  Knsemble  Averaging 

Under  conditions  of  input  data  stationarity,  FFT  output  stability  may  be  enhanced  by  aver- 
aging the  PSD  outputs  of  separate  time  segmentt;.  This  procedure  consists  of  sectioning  a time  interval 
of  interest,  applying  a separate  data  window  to  each  section,  transforming  each  section,  spectral 
averaging  each  section's  PSD,  and  forming  the  ensemble  average  of  the  smoothed  periodograms 

4.1  1 .5  Normalization 

The  output  of  the  ensemble  average  algorithm  is  normalized  to  account  for  the  effects  of  the 
input  data  window  and  the  spectral  averaging  process.  To  account  for  these  two  effects,  the  peri- 
odogram's  PSD  is  multiplied  by  the  ratio  of  A|/A-,  and  divided  by  the  quantity  2M+1,  where 

N 

=^>  r i 

i=i 
N 

(t)dt=C>^  W'(i)  (4-12) 


4.12  INTEGRAL  AND  FOURIER  OPERATIONS  ON  S(t.cj) 

Both  the  noise-present  and  noise-free  estimates  of  the  scatter  channel's  delay-Doppler  function 
may  be  processed  by  an  algorithm  that  performs  integral  and  Fourier  operations  to  yield  delay  spectrum. 
Doppler  spectrum,  joint  time-frequency  correlation  function,  frequency  autocorrelation  function,  time 
autocorrelation  function,  and  mean  square  total  multipath  power.  These  operations  are  described 
below. 


4. 1 2. 1 Delay  Spectrum 

For  the  case  where  the  Sir. to)  function  is  continuous  over  both  the  r ami  to  variables,  the 
channel's  delay  spectrum  may  be  derived  through  the  following  integral: 


In  our  case.  Str.ui)  is  discrete  in  both  dimensions;  thus,  the  delay-spectra  estimate  is  given  by 


0(rj)  = • (4-13) 

where  Q(rd  is  the  power  spectral  density  ot'energ>’  arriving  with  delay  . 

As  previously  described  in  section  4.‘>.  a time-domain  analysis  on  each  tap  output  may  be  used 
to  derive  an  estimate  of  QItj).  That  analysis,  however,  is  limited  to  the  noise-contamined  signal 
structure  whereas  the  integral  operation  on  S(t,ui)  allows  one  to  operate  on  noise-removed  data  arrays. 


4. 1 2.2  Doppler  Spectrum 


The  estimate  of  the  multipath's  Doppler  spectrum  is  derived  by  summing  the  delay-Doppler 
function  over  the  delay  variable: 


D(u)j)  = Z)S(Tj.cjj).  (4-14) 

i 

where  D(Wj)  is  the  energy  arriving  at  the  receiver  with  Doppler  frequency  shift  coj  . 

In  general,  the  Doppler  spectrum  is  highly  asymmetrical  and  thus  both  the  negative  and  posi- 
tive frequency  components  of  D(cj)  are  maintained.  The  discrete  outputs  are  given  at  the  following 
spectral  locations: 


cjj  = (2wAf)j.  radians/sec 
SR 

Af  = VI  hertz  . 

N 

where: 

SR  = sampling  rate  of  S.ACP  receiver  laps 
N = number  of  points  in  I TT  (sec.  4,1  I ) ’ 

M = window  width  of  spectral  smoothing  (sec.  4.1  I ) . 


(4-15) 


These  outputs  fall  between  the  limits  of  iSR,  2. 
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Joint  Time-FrequcnLy  Autocorrelation  Function 

In  tliat  the  channel's  delay-Doppler  scatter  function  and  joint  time-frequency  autocorrelation 
function  are  Fourier  transform  pairs,  we  derive  an  estimate  ot  RtJ.fl)  through  application  of  a 
double  inverse  FFT  to  the  S(r.cj)  array: 


R(|, 


S' 


M-1 

= £ 

m=0 


N-l 

£ 

n=0 


S(u.,^.Tn)exp 


( xm 

yn^ 

•.71  I 

^ M 

Ny 

x = 0.1, 
y = 0.1, 


. M- 
. N- 
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where: 

time-frequency  autocorrelation  function  at  time  frequency  0,^ 

M = number  of  elements  in  Doppler  array  extended  to  a value  that  is  an  integer  power  of 

2 (zero-fill  if  necessary) 

N = number  of  elements  in  delay  coordinate  array  extended  to  an  integer  powv.i  ui  ; 

(zero-fill  if  necessary) 

. = 


Parameters  and  Hy  take  on  the  following  discrete  values: 


t = — — : X = 0.±1 , 

AfM 


M-1 


y N-l 

^T- 


where: 

Af  = frequency  separation  between  adjacent  Doppler  spectra  outputs 
r = the  effective  tap  width  of  the  correlator  output . 


Rt^.ni  is  output  (magnitude  and  phase)  in  both  numerical  and  three-dimensional  plotted 
form.  Since  Rt^.Sl)  = R*(-t  >•  only  half  of  the  periodogram  output  is  retained. 


I 

4.12.4  Frequency  Autocorrelation 

An  estimate  of  the  channel’s  frequency  autocorrelation  function  is  derived  by  applying  the 
inverse  FFT  operation  to  the  delay  spectrum 


R(0.S2:)  = j = 0, 

k=0 


,N-1 


(4-17) 


where: 

R(0,f2j)  = frequency  autocorrelation  function  output  at  frequency  f2j 

N = number  of  delay  taps;  zero  order  extended  to  be  equal  to  the  nearest  value  that  is  an 

integer  power  of  2 

i =V^ 

Since  Q(T|^)  is  real  valued  for  all  T|^,  the  R(0,f2)  distribution  is  symmetric  with  respect  to 
n = 0.  Thus,  the  algorithm  outputs  only  the  positive  components  of  the  distribution,  which  are  given 
at  the  following  frequency  locations: 

flj  = j(l/rN)  j = 0.1 N/2  , 

where  r is  the  effective  tap  width  of  the  correlator  output.  Outputs  are  given  in  both  numerical  and 
plotted  formats. 


4.12.5  Time  Autocorrelation 

Application  of  the  inverse  FFT  algorithm  to  the  multipath  channel  Doppler  spectrum  provides 
an  estimate  of  the  time  autocorrelation  function.  Tliis  is  represented  by 

R(£i.O)  = S j = o.l N-1  . (4-18) 

^ k=0 


“As  a program  option  R({.f2)  may  be  bypassed;  thus  we  derive  the  zero-axes  cuts;  i.e.,  R(J,0)  and  R(0.S2). 
via  a method  that  does  not  depend  on  the  availability  of  R(J.f2). 
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where : 

R(Jj.O)  = time  autocorrelation  function  at  time 

M — number  ot  D(w)  array  input  points  extended  to  a value  that  is  an  integer  power  of  2 
(zero-fill  if  necessary) 
i = v^. 

The  R(Jj,0)  function  is  symmetric  with  respect  to  J = 0 since  D(cO|.)  is  real  valued.  Thus, 
only  the  positive  components  of  R(^0)  are  output.  These  are  given  at  the  following  discrete  time 
intervals: 


?j  = j/AfM  ; j = 0.1 M-1, 

where  Af  is  the  frequency  separation  between  adjacent  Doppler  spectra  outputs.  The  output  formal 
of  the  complex  R(J,J2)  array  is  identical  to  that  used  to  represent  RlJ.H). 


4.12.6  Mean  Square  Total  Multipath  Power 


The  multipath  channel’s  mean  square  power  is  derived  by  integrating  the  delay  and  Dopper 
spectra  over  their  respective  variables;  i.e.. 


N 

S D(oj.) 
i=l 

' M 

H 0(Ti) 
. i=l 


(4-1  dt 


where: 

N = number  of  periodogram  outputs 
M = number  of  delay  taps  . 

Both  methods  are  used  to  give  a partial  validity  check  on  the  spectral  transformation  operations  that 
have  occurred  up  to  this  point  in  the  overall  computer  program  (i.e..  if  the  alternate  procedures  yield 
different  results,  then  a program  error  is  indicated). 

It  is  noted  that  the  <|l|‘>  estimation  technique  is  in  essence  equivalent  to  an  incoherent 
power  summation  ot  the  delay  t:ip  outputs.  This  does  not  give  a t rue  estimate  since . ( I ) adiacent  t.ips 
(for  the  usual  one-chip  tap  spacing)  have  a .SO'r  commonality  in  the  delay  area  to  which  lhe\  are 


partiall>  correlated  and  (2)  the  autocorrelation  function  of  the  PN  m-sequence  is  not  uniform  over  a 
chip  width  hut  has  trianjiular  distribution  with  peak  equal  to  the  sequence  length  and  base  equal  to  -1 
at  locations  ±1  chip  from  the  peak.  Assuming  the  multipath's  delay  spectra  is  piecewise  constant  over 
a chip  width,  the  above  effects  are  easily  calculated  to  yield  an  <|llS  estimate'^  that  is  on  the  order 
of  1 .8  dB  lower  than  would  be  the  case  for  a coherent  tap  summation.  A coherent  summation,  in 
effect,  provides  a single  tap  experiment  with  uniform  peak  correlation  over  the  entire  tap  bank  range. 
The  tap  lumping  (contraction)  feature  of  the  program  allows  this  to  be  performed;  however,  it  cannot 
be  employed  simultaneously  with  the  detailed  single  tap  resolution  SIt.cjI  calculation. 


4. 1 3 spre;ad  ( alculations 

Spread  calculations  are  performed  on  the  multipath  channel's  delay  spectra.  Doppler  spectra, 
frequency  autocorrelation  function,  and  time  autocorrelation  function.  The  algorithm  operates  on 
these  distributions  to  determine  the: 

a.  .Maximum.  3-dB.  and  I /e  levels 

b.  Spread  between  the  3-dB  levels  j 

c.  Spread  between  the  I /e  levels  J 

d.  First  moment  of  the  distribution  1 

e.  Second  moment  of  the  distribution  with  respect  to  the  mean  level  j 

f.  Second  moment  of  the  distribution  with  respect  to  the  maximum  level  ; 

1 

1 

With  the  exception  of  the  second  moment  parameters,  these  quantities  are  diagrammed  in  figure  4-8 
for  an  assumed  Doppler  spectrum. 

The  spectra  and  autocorrelation  functions  upon  which  the  spread  calculation  algorithm  oper- 
ates are  in  the  form  of  discrete  sample  arrays.  Assuming  an  array  distribution  of  the  form  yj  = f(Xj). 

the  processing  begins  with  a search  for  the  Xj  value  that  produces  the  maximum  y'j  . These  values  , 

are  designated  by  and  respectively.  Two  searches  are  then  originated  from  the 

location  (one  with  i decreasing,  the  other  with  i increasing)  to  determine  the  Xj  values  where  the 

distribution  first  crosses  the  yj.,,„v  '2  and  yi.„,  ^/e  levels.  The  corresponding  differences  between 

*mjx  iTiux 

the  Xj  '2  and  xj  e levels  are  then  computed  to  yield  the  3-dB  and  I A' delay  spread  values, 
niiix  niux 


-^Assuming  an  infinile  sampling  rate  in  the  digital  implementation  of  (he  SACP  receiver. 
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Figure  4-8.  Spread  Parameter  Illustration 

In  deterniinini!  the  distribution  mean  value,  the  first  step  is  to  calculate  the  total  value  of  the 
array  sum : 


N 

^=E 

1=1 


(4-:o) 


where: 

S = array  sum 

N = number  ol  points  m the  arra\. 

The  mean  value  of  the  array  's  \-coonlinale  is  then  calculated  as 


P 


E 


(4-: 


( « 
M 


t 


Two  I'orms  of  second  moment  calculations  are  performed,  one  with  respect  to  the  array's 
mean  value  and  the  other  with  respect  to  the  quantity: 


(4-22t 


(4-23) 


Note  that  the  algorithm  operates  on  both  sides  of  the  delay  and  Doppler  spectra  but  considers 
only  the  positive  x-values  for  the  autocorrelation  function. 


4.14  NOISE  DETERMINATION  AND  REMOVAL  (NDandR) 

The  outputs  of  the  SACP  multipath  correlator  contain  desired  signal  data,  low-pass  additive 
thermal  noise,  spurious  signal  terms,  and  receiver  arithmetic  noise.  Determination  of  the  bulk  of  the 
tap  noise  is  aided  by  two  factors:  ( I ) for  a given  time  delay,  there  are  upper  and  lower  Doppler  fre- 
quency limits  beyond  which  it  is  physically  impossible  for  multipath  returns  to  exist  and  (2)  .SACP 
receiver  operator  procedures  entail  replica  sequence  delay  settings  such  that  the  specular  pa'li  return  is 
captured  somewhat  after  the  lO’^  tap  (oceanic)  or  30**’  tap  (overland  CONUS)  in  the  ban)..  See 
figure  4-9  for  a plan  view  of  these  realms  as  superimposed  upon  the  delay-Doppler  plane  of  the  SACP 
receiver.  Also  included  are  regions  X and  Y,  which  are  referred  to  in  the  following  discussion. 

The  algorithm  for  noise  determination  and  removal  is  illustrated  in  figure  4-10.  Delineation  ol 
the  major  features  of  each  processing  step  is  given  below. 

4.14.1  Eilter  Erequency  Response  Determination 

The  SACP  complex  correlators  use  digital  realizations  of  Butterworth-Thompson  low-pass 
filters  at  their  outputs.  Their  Irequency  response  is  determined  in  the  NDandR  process  by  emulating 
the  recursive  filter  algorithm  for  the  case  of  interest  and  then  subjecting  it  to  a series  of  unit  magni- 
tude sinusoidal  inputs  that  occur  at  the  frequencies  of  the  IT'1  periodogram  for  the  S(r,to(  function 
The  amplitude  square  of  the  filter  output  is  used  to  represent  the  magnitude  of  the  discrete  tfequcncN 
response  of  the  low-pass  filter.  This  is  liesignated  ERt  f ). 
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Figure  4-9.  Delay  Doppler  Realm  of  SACP  Receiver 


4. 14.2  Filter  Response  Including  Aliased  Components 


\ 


When  dealing  with  a tap’s  low-pass  thermal  noise  content,  we  consider  not  only  the  basic 
trequency  response  of  the  filter  but  also  the  aliased  components  cau.sed  by  the  implicit  periodicity 
assumption  of  the  FFT  algorithm.  We  include  the  first  four  periods  of  the  aliased  low-pass  filter  noise, 
which  is  rather  significant  due  to  the  effective  tap  sampling  rate  at  only  twice  the  filter's  -I  dB  fre- 
quency. Assuming  a white  thermal  noise  content  in  the  SACP  receiver,  the  filter  response  to  thermal 
noise  including  the  bulk  of  the  aliased  components  may  be  expressed  as 


iM(f)  = 53  FR(f+iAf),  (4-24) 

i=-2 

where: 

M(t)  = amplitude  ot  filter  response  to  white  noise  including  first  four  aliased  components 
Af  = basic  tap  sampling  rate  . 


4. 1 4.3  Multipath-Free  Total  Noise  Spectra 

Once  a multipath-free  delay  region  has  been  determined  (i.e.,  region  X of  fig.  4-d).  its  total 
noise  spectra  is  determined  by  averaging  the  S^(r.cj)  array  over  the  delay  variable: 


N(f) 


, XN 

= - E 


N i=Xl 


Sn'^i 


(4-25) 


where: 

N(f)  = multipath-free  region’s  total  noise  spectra 
XI  = first  tap  of  multipath-free  region 
XN  = last  tap  of  multipath-free  region  . 


An  example  of  a multipath-free  region’s  total  noise  spectra  is  given  in  figure  4-1 1.  For  this 
case,  the  low-pass  filter  bandwidth  was  on  the  order  of  300  Hz  and  we  note  the  effect  of  the  , lased 
thermal  noise  (i.e..  the  increase  in  the  density  at  the  tails  of  the  distribution).  Also  in  evidence  is  a 
large  spurious  dc  noise  component,  which  is  caused  primarily  by  the  SACP  receiver’s  two’s-comple- 
ment  arithmetic. 
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Figure  4- 1 1.  Multipath-Free  Noise  Spectrum  Showing  Aliasing  and  Arithmetic  Noise 
4.14.4  Average  Thermal  Noise  Level  of  Nd) 

A tap’s  thermal  noise  psd  may  be  represented  by  AgM(0,  where  is  tap  dependent.  To 
determine  A^,  we  analyze  the  tap’s  psd  in  a frequency  region  devoid  of  multipath  (i.e.,  region  of 
tig.  4-9).  This  region  is  also  well  removed  from  the  receiver  arithmetic  noise  and  the  PN  code’s  low- 
pass  spurious  response  (which  has  bandwidth  characteristics  of  the  multipath  signal).  Thus,  we  may 
assume  it  to  be  contaminated  by  thermal  noise  alone.  For  N(f),  the  A^  estimate  is  obtained  from 
the  following  operation: 

YN 

E N(fj) 
i=Yl 


where: 

Y1 

YN 


= first  spectral  component  in  region  Y 
= last  spectral  component  in  region  Y 
YN 

= E M(f.). 

i=YI 


4.14.5  Residual  Noise  Determination 

By  subtracting  Aj,M(f)  from  N(f),  we  obtain  what  is  referred  to  as  the  tap  bank  residual 
noise  estimate: 


Rtf)  = Ntf)- A^M(f)  . 
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This  parameter  thus  represents  the  nonthermal  noise  content  of  the  psd.  It  is  assumed  to  be 
identical  for  eacli  tap  in  the  bank.  The  procedures  for  calculating  the  residual  noise  includes  the  PN 
correlator’s  low-pass  spurious  response  to  each  and  every  signal  that  arrives  at  the  receiver.  When  the 
residual  noise  is  subtracted  from  a tap  containing  multipath  return,  we  are  thus  removing  a portion  of 
low-pass  spurious  noise  that  is  not  present  in  that  particular  tap  initially.  For  example,  if  the  tap  in 
question  intercepts  delay  energy  E(Tj),  then  we  remove  a component  corresponding  to  the  noise 
determination  interval’s  uncorrelated  PN  response  to  Elrj).  For  the  sequence  length  typically  used  in 
the  multipath  test  (5  1 2),  this  error  is  very  small  and  amounts  to  roughly  four  parts  in  10*’  (power)  of 
the  correlated  signal. 


4.14.6  Residual  Noise  Removal 


Once  the  residual  noise  estimate  is  determined,  the  NDandR  algorithm  then  enters  into  a tap 
loop  and  proceeds  to  derive  a noise-free  estimate  of  the  channel's  delay-Doppler  function.  The  first 
step  involves  removing  R(f)  from  the  particular  tap  in  the  loop: 


Sj(ri,f)  = S,^(rj,f)-R(f), 


(4-28) 


where: 

Sj(Tj,f)  = the  residual  noise-free  estimate  of  the  i^*’ tap’s  psd 
Snl^i-f)  ==  the  noise-contaminated  estimate  of  the  i^*i  tap’s  psd  . 


4.14.7  Average  Thermal  Noise  Level  of  Sj(f|,f) 


After  the  residual  noise  content  has  been  removed  from  a particular  tap’s  psd,  its  thermal  noise 
level  AqItj)  is  derived  in  a manner  identical  to  that  discussed  for  the  N(t)  ease;i.e.. 


V^i'  = 


YN 

j=YI  ■' 

SM 


(4-2‘)) 


where  Aj,(Tj)  is  the  thermal  noise  level  of  the  i**’  delay  tap’s  psd. 
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4.14.8  Thermal  Noise  Removal  and  Filter  Attenuation  Compensation 


rite  final  steps  in  NDandR  involve  removing  .•\y(Tj)M(f)  from  Sj(Tj.f)  and  dividing  the  result 
by  the  amplitude  response  of  the  low-pass  digital  filter  ( whieh  compensates  for  the  filter  attenuation 
on  the  multipath  signal  component).  Tliis  is  represented  by: 

S-p(7j.f)-M(f)  xA^(t.) 

S(T;.f)  = — • (4-40) 

‘ N(f) 

where  StTj.f)  is  the  noise-free  filter  compensated  estimate  of  the  multipath  psd  for  the  i**’  receiver  tap. 

The  effectiveness  of  the  NDandR  algorithm  is  illustrated  by  the  sample  data  presented  in 
section  4.44  (figs.  4-20  through  4-24).  In  these  figures,  the  delay-Doppler  scatter  functions,  delay 
spectra,  and  Doppler  spectra  for  both  noise-present  and  noise-removed  conditions  are  illustrated. 
Visually,  very  little  difference  appears  to  exist  between  the  noise-present  and  noise-removed  estimates 
of  S(r.co’)  (although  for  lower  signal-to-noise  cases,  the  dc  offset  bias  for  the  noise-removed  case 
exhibits  the  expected  random  uncertainty  tluctualion  about  the  zero-power-level  x-y  axis).  On  the 
other  hand,  there  is  a pronounced  difference  between  noise-present  and  noise-free  Doppler  and  delay 
spectra.  The  noise-free  delay  spectra  have  essentially  no  energy  prior  to  the  specular  return,  approach 
zero-energy  asymptotically  as  the  delay  variable  becomes  much  greater  than  the  specular  return,  and 
do  not  possess  the  extraneous  noise  in  tap  81.  By  comparison,  the  noise-free  Doppler  spectra  estimate 
IS  virtually  signal  free  for  frequencies  greater  than  roughly  I 50  Hz  (the  physical  upper  limit  for 
Doppler  return  over  the  tap  bank  extent  for  the  test  condition)  and  the  dc  spurious  noise  contribution 
has  been  completely  eliminated. 

In  general,  the  NDandR  algorithm  will  handle  most  noise-contamination  cases  encountered  in 
the  multipath  tests.  It  will  handle  tap-dependent  thermal  noise  conditions  (e.g..  tap  81  in  the  above 
discussion  I but  will  not  properly  account  for  conditions  where  the  statistics  of  the  residual  noise  content 
vary  from  tap  to  tap.  The  latter  cases  have  been  observed  to  occur  during  conditions  of  receiver 
malfunction  and  when  unwanted  signals  have  entered  the  tap  bank. 


4. 1 5 AIRBORNl  PROBFR  ANTFNNA  IT  FI  CTS  RFMOVAL 

I'o  remove  the  prober  aircraft  antenna  pattern  effects  from  the  channel  measurements,  it  is 
first  necessary  to  make  some  assumptions  regarding  the  mechanism  b>'  which  the  electromagnetic 
waves  are  scattered  from  the  earth  surface.  In  line  witli  the  most  witlely  accepted  theoretical  treat- 
ments and  the  good  agreement  between  experiment  and  theory  as  presented  in  volume  V.  we  base  our 
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analysis  upon  the  vector  formulation  of  the  physical  optics  very  rough  surface  scatter  rr.odel.  Under 
this  assumption  the  received  power  from  an  increment  of  surface  area  is  expressed  as 


dPrcvd“'Tcf'‘T?f)Pss'%-^e' 


(4-31) 


where: 


dPrcvj  “ power  received  from  the  incremental  area 


^ss'^rj'^e*  “ probability  of  incremental  area  having  slope  6^  and  6^  in  the  rj  and  e direc- 


tions. 


Tcf 


IS.  respectively  and  6^  provide  favorable  reflections  from  the  surface  to 
satellite) 

= complex  coefficient  that  accounts  for  coupling  between  incident  polarization 
vector,  tilted  surface  facet,  and  receiver  polarization  vector. 


The  aircraft  antenna  influence  is  embodied  only  in  parameter  which  may  be  represented 

as  (see  sec.  5.2): 


= R||  (a)  1 

i .A  cos  8j  + B sin  6j  j 

j C cos  6^  + D sin  6^  j 

-Rp (a) 

j B cos  5j  - A sin  5j  j 

j-D  cos  5,  + C sin  6..  ] 

(4-32 

where : 

Rii  (al 


Rj^  (a) 


«s 


A.BC.D  = 


parallel  Fresnel  reflection  coefficient 
perpendicular  Fresnel  reflection  coefficient 

angle  between  incident  wave  and  normal  of  properly  tilted  (to  produce  reflection 
into  receiver)  surface  facet 

angle  between  theta  component  of  incident  wave  and  the  incident  parallel  unit 

vector  (i.e..  vertical  polarization  with  respect  to  tilted  surface) 

angle  between  theta  component  of  scattered  wave  and  the  reflected  parallel  unit 

vector 

transmitter  and  receiver  complex  antenna  polarization  vector  coefficients  as 
described  below  . 


Polarization  vectors  tor  the  transmit  (Pjl  and  receive  (Pj^)  antenna  system  are  given  by : 


P p A i-j-p  + B ipp 


Pr  - C i p^  + I)  ip^ 
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(4-33) 


1 


i 

. 

' 

i 


where: 

^TT'^PT  ~ 'f*  ’I’*-'  Pi'i  tlirections.  respectively,  with  respect  to  coordi- 

nates centered  on  the  transmitter 

^TR'^PR  ~ theta  and  phi  directions,  respectively,  with  respect  to  coordi- 

nates centered  on  the  receiver. 

In  our  case  transmission  originates  from  the  kC-135  and  is  received  at  the  ATS-t>  satellite.  The 
satellite  antenna  is  assumed  to  be  right-hand  circular  (RHC)  with  0-dB  ellipticity : thus: 

C = 

D = jvcy: , 

where  ('i^  is  the  gain  magnitude  of  the  satellite  antenna. 

Parameters  .-\  and  B for  the  aircraft  antenna  pattern  are  derived  from  scale-model  antenna 
range  measurements  (see  vol.  \'ll  of  this  report).  We  apply  the  antenna  effects  removal  algorithm  to 
the  specific  case  where  the  prober  is  operating  in  the  simultaneous  horizontal 'vertical  polarization 
mode.  These  probe  polarization  descriptions  are  nominal  in  that  for  aircraft  antenna  systems  one 
usually  expects  the  presence  of  cross-polarization  tenns.  Thus  for  the  nominal  horizontal  polarization 
probe,  B represents  the  principal  polarization  term  and  A is  a measure  of  the  cross-polarization  tenn 
(vice  versa  for  nominal  vertical  polarization).  From  the  antenna  range  radiation  distribution  plots,  it 
is  observed  that  for  the  various  prober  antenna  pointing  directions  the  ratio  of  principal-to-cross- 
polarization  terms  is  typically  on  the  order  of  20  dB  over  the  effective  "glistening"  area  of  the  surface. 
Hence,  to  good  approximation  we  consider  only  the  principal  tenns  of  the  polarization  vectors  and 
write  T^,|-  for  the  horizontally  polarized  probe  as: 


R|l  (a)  sin  5j  (cos  + j sin  6^)  - Rj^  (a)  cos  5j  (sin  5^,  - j cos  6^)  . 


(4-.U) 


and.  for  the  vertically  polarized  probe. 


I 


cf 


(a)  cos  6j  (cos  8^  -)-  j sin  8^)  + R ^ (n)  sin  8^  (sin  8^ 


I cos  8^^) 


(4-.’'.si 


where: 

G|l  = gain  of  horizontally  polarized  probe 
Ciy  = gain  of  vertically  polarized  probe. 


4-,’^4 


Parameters  Gjj  and  Gy  do  not  have  constant  gain  over  the  scattering  surface.  To  remove 
this  spatial  filtering  effect,  we  operate  in  the  delay-Doppler  coordinate  system  of  the  receiver.  The 
scatter  function.  S(t,w),  is  related  to  dP|.(-v(j  through  the  following  formulation: 


S(r,w)  = J|(T,w)dP^^.yj  , + 2 ’ <4-361 

where  jfr.w)  is  the  Jacobian  of  the  transformation  from  surface  coordinates  to  (t,w)  receiver 
coordinates. 

In  the  above  equation,  subscripts  I and  2 denote  the  fact  that  the  scattering  operation  is  a 
2-into-l  mapping  process  (that  is.  two  surface  points  are  mapped  into  one  delay-Doppler  location).  In 
general,  even  for  isotropic  antennas.  dP^^,yj  , does  not  equal  thus  without  making  some 

a.ssumptions  with  regard  to  the  distribution  it  becomes  impossible  to  properly  adjust 

S(T.o))  for  the  antenna  filtering  effects.  To  circumvent  this  ambiguity,  we  restrict  our  attention 
entirely  to  the  data  gathered  during  the  in-plane  legs  of  the  test  scenario.  The  symmetry  properties  of 
the  delay  contours,  Doppler  contours,  and  function^  are  thereby  exploited  and  we  may  write 

S(r,cj)  as: 


S(T.w)  = J|<^-UJ)dP;.^yj  I (G,  +G^)  . (4, 

where  tlPJ-j-vd,  ~ ‘^*’rcvd'^^*^*'  ' c.,  the  received  engery  from  incremental  area  (j3)  given  that  the 
antenna  illumination  is  isotropic. 

Now  under  the  same  geometry  and  surface  assumption  we  may  express  S(r,w)  for  the  iso- 
tropic antenna  g.iin  case  as: 


S,(r.w)  - Ji<^.w)dP;^,^j  , +J,(r.w)dP;^,^j  , 


= 2J,(r.w)dP'^,^j^|. 


4 • 

Kxamination  of  equation  (4-31 ) reveals  that  for  in-planc  Hight,  dP^^^^^j  will  be  cross-plane  symmetri- 
cal (apart  from  antenna  asymmetry  effects)  if  ( I ) the  sea-surface  slope  distribution  is  isotropic,  a valid 
assumption  per  data  from  Cox  and  Munk  (ref.  4-2)  and  (2)  the  polarization  planes  of  the  transmit 
and  receive  antennas  are  linear  and  lie  either  in  or  perpendicular  to  the  plane  containing  the  KC-l3,‘i. 
specular  point,  and  ATS-(i  satellite.  The  latter  restriction  is  necessary  to  ensure  that  T ,j  j = T . 

In  our  case  this  condition  is  not  met  since  the  satellite  antenna  polarization  is  circular.  ^^owev^,^.Tl 
can  be  shown  that  for  grazing  angles  greater  than  a few  degrees  the  induced  error  is  relatively  small 
when  compared  to  other  factors  such  as  the  antenna  gain  characterization. 
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Thus  for  in-plane  llight  the  isotropic  antenna  scatter  function.  S|(r.aj).  is  related  to  the  antenna-gain- 
contaminated  version  S(r.co)  by  the  simple  relationship: 


S|(r.u)) 


SiT.UJl 

(ilT.Wl 


where; 


Cl  I G 1 

Ci(T.U))  = ; 


(4-39) 


(440) 


In  the  commonly  used  steepest  descent  integration  analysis  of  forward  scatter  multipath 
(mathematically  appropriate  under  the  constraint  that  the  grazing  angle  is  much  greater  than  the  sur- 
face rms  slope),  one  usually  assumes  that  the  antenna  coverage  is  adequately  described  by  its  specular 
point  gain.  For  this  reason  we  have  chosen  to  relate  our  antenna-corrected  version  of  the  scatter  func- 
tion to  the  gain  directed  toward  the  specular  point;  i.e.. 


Syy(^(r.tJ)  = S|(T.cj)  G(O.O) 


Slr.tc) 

G(r.cJ) 


G(O.O)  , 


(441  ) 


where: 

S^(^  (r.cj)  = antenna-corrected  version  of  the  scatter  function 
G(O.O)  = the  specular  point  antenna  gain  (i.e..  r = CO  = 0) . 

In  addition  we  also  calculate  a difference  function.  SSiT.co).  between  the  corrected  and  iin- 
corrected  versions  of  the  scattering  lunction: 


S(r.cj)  = S(T.w)  - S^^  tT.'uj)  . 


(442) 


I'his  function  is  available  as  a three-dimensional  plot  and  provides  a compact  ilescription  ol  the 
prober's  spatial  filtering  characteristics  as  translated  'o  the  delay-|)oppler  coordinates  ol  the  receiver 

['he  algorithm  for  delemiining  G(t.co)  has  the  following  parameters  as  its  basic  input  ( I ) air- 
craft antenna  pattern  in  the  classical  theta,  phi  coordinate  system.  (2)  aircraft  altitude,  velocity,  and 
height.  (3  ) elevation  angle  to  the  satellite,  and  (4  ) tap  number  associated  with  specular  point  return 


I 


4 
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Using  this  information  the  program  determines  the  realm  of  physically  possible  (r.w)  point 
pairs  into  which  the  scatter  surtace  can  be  mapped.  Each  S(r.w)  density  within  this  realm  is  then 
compared  to  a noise  threshold  level  (to  bypass  this  test,  the  threshold  may  be  set  to  zero).  If  the 
threshold  test  is  not  passed.  S^(;-(r.a;l  is  set  to  zero:  otherwise  standard  geometric  techniques  are 
used  to  map  (t.oj)  into  the  two  sets  ot  surtace  coordinates  that  are  responsible  for  the  scatter.  These 
locations  are  in  turn  translated  into  antenna  pattern  coordinates  from  which  a table  look-up  routine  is 
used  to  derive  the  gains  Gj  and  G^.  which  are  then  used  in  equation  (4-40)  to  form  the  G(t.w) 
function.  As  a secondary  output  from  the  algorithm,  one  may  obtain  GIt.cj)  in  both  numerical  and 
three-dimensional  plotted  formats. 


4.16  TAP  I AND  0 .AMPLITUDE  DISTRIBUTION 

Amplitude  probability  distribution  in  the  form  of  histograms  is  derived  for  the  I and  0 com- 
ponents ot  a given  delay  tap  through  the  application  of  standard  numerical  grouping  techniques.  A 
block  diagram  illustrating  the  fundamental  components  of  the  APD  algorithm  is  given  in  figure  4-1  2. 
Prior  to  histogram  sorting,  three  operations  are  performed  on  the  input  data: 

a.  Tlie  data  may  be  optionally  desampled  to  improve  the  independency  of  the  sample  space 
array . 

b.  A dither  signal  is  added  to  each  of  the  desampled  data  points.  This  is  necessary  since  the 
quantizing  levels  of  the  digital  SACP  receiver  are  not  significantly  smaller  than  the  bin 
widths  of  the  histograms.  In  accord  with  the  two's-complement  arithmetic  of  the  receiver, 
the  dither  signal  is  chosen  from  an  ensemble  whose  values  are  uniformally  distributed  be- 
tween 0 and  1 . A new  dither  signal  is  used  for  each  signal  input  value. 

c.  The  mean  and  standard  deviation  of  the  desampled.  dithered  array  is  determined: 


(4-4.4) 


a = 


(4-44) 


4-.47 


Tap  I or  Q 
component 


Numerical  output:  histogram 
value  and  X probability 


Figure  4- 12.  Tap  / and  Q Amplitude  Distribution  Algorithm 
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After  the  sample’s  mean  and  standard  deviation  have  been  determined,  the  algorithm  then 
sorts  the  input  sample  array  into  three  histogram  distributions.  All  three  histograms  are  calculated  for 
each  tap  component.  The  sorting  techniques  are  described  as  follows: 

a.  Histogram  1 - consists  of  20  bins  ranging  between  -a>  and  +=o  , with  the  18  interior  bins 
being  of  equal  width  and  spanning  the  range  [-3a,  +3a[ . 

b.  Histogram  2 — similar  to  No.  1 except  bins  are  combined  with  their  nearest  neighbors  if 
the  observed  frequency  is  less  than  a specified  (input)  percentage. 

c.  Histogram  3 — bin  widths  are  selected  so  that  the  expected  frequency  (see  below  for  dis- 
cussion pertaining  to  expectation  model)  for  each  bin  is  S7f.  Thus,  there  are  20  bins  total. 

The  three  histogram  distributions  are  then  tested  for  normality.  This  is  carried  out  via  the 
X“  goodness-of-fit  test.  The  x"  statistic  is  calculated  as: 


^ (EFj-OFj)- 


i=l  EFi 


(445) 


where: 

N 


FFi 

OF: 


number  of  bins  in  the  histogram 

the  expected  number  of  data  points  falling  into  the  i*E  i^jj, 
the  observed  number  of  data  points  falling  into  the  i*E  bin 


Parameter  FFj  is  based  on  the  assumption  of  a Gaussian  distribution  whose  mean  and 
variance  are  derived  from  the  sample  space  array ; i.e.. 


cn  - ' f*' 

FFj  - I expl )dx  , 

V,  \2o^  / 


(44(1 ) 


where  tj  is  the  upper  limit  of  the  voltage  range  associated  with  the  i**’  bin. 

Knowing  the  x~  statistic  and  the  number  of  degrees  of  freedom  (N-3)  = K associated  with 
each  particular  histogram,  we  are  then  able  to  detennine  the  following  cumulative  probability: 


F(x-)  =J 


2'^nK/2)  o 


^dy 


(447) 
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This  represents  the  probability  that  a sample  space  selected  frotn  a Gaussian  distribution  with  the  same 
mean  and  variance  as  the  test  case  will  produce  a x"  value  less  than  that  of  the  observed  sample  (i.e., 
if  F(x‘i  i'  consistently  greater  than  0.9,  we  would  have  reason  to  doubt  the  nonnality  of  the  data). 


! i 


4. 1 7 TAP  PHASt;  DISTRIBITION 


We  define  the  phase,  4>.  of  the  complex  tap  received  signal  as 


tan 


Oi/Ii 


(448) 


The  algorithm  for  determining  the  probability  distribution  histogram  of  4>  operates  on  the 
same  desampled  dithered  1 and  0 data  values  that  are  input  to  the  amplitude  statistics  routine. 
Twenty  equal-length  bins  are  utilized  for  the  histogram  structure,  which  has  a total  extent  from  -180° 
to +180°.  Via  the  x"  goodness-of-fit  test  as  described  in  section  4. lb,  the  observed  bin  densities  are 
compared  to  a hypothesis  resulting  from  the  two-dimensional  Gaussian  distribution  whose  quadrature 
components  are  independent  with  nonzero  means  and  nonidentical  variances: 


W(1,0) 


:trO|OQ 


exp 


(i-/i,)-  (Q-tiq)- 


(449) 


where; 

W(1,0)  = joint  probability  density  of  the  tap  output  taking  on  the  specific  I and  Q values 
°r  °0  ~ standard  deviations  of  the  I and  Q processes 
/i|,jiQ  = means  of  the  1 and  Q processes. 


The  phase  probability  density  may  be  derived  by  transferring  the  above  relationship  to  polar 
coordinates  and  integrating  over  the  total  extent  of  the  radial  dimension;  i.e., 


P(0)  = 


:rra|OQ 


-/”fdre> 


xp  - 


(f  cos  0 - /aj)-  (fsin0-fjy)- 


Jo 


0 


(4..S0) 


Knowing  P(0).  the  expected  bin  frequencies  are  then  calculated  by  integrating  over  the  (0) 
range  and  using  steps  similar  to  those  described  in  the  previous  section  to  yield  the  x“  statistic  and 
cumulative  probability. 
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4.18  TAP  PROCESS  CROSS-CORRKLATION 
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Given  any  two  complex  tap  processes,  this  algorithm  computes  the  normalized  cross-correla- 
tion function  as: 


RxY<k) 


£ ^i^ilk 

1=1 


N N 


'A 


(4-51) 


where: 

*^XY**^*  = normalized  cross-correlation  function  between  tap  processes  X and  Y for  lag 
offset  k 

Xj  = i*^  sample  of  the  X tap  process  (as  a program  option  the  mean  may  he  removed 

prior  to  processing)  . 


The  program  is  written  so  that  .X  and  may  he  taken  from  either  the  same  tap  bank  or 
from  cross-polarized  tap  banks,  A high-level  How  diagram  for  this  algorithm  is  shown  in  figure  4-13. 
As  indicated,  there  are  two  major  loops  in  the  routine.  The  case  loop  is  used  to  delineate  the  tap 
correlation  pairs  and  the  number  of  points  to  be  considered  in  each  calculation.  For  example,  one 
may  request  that  all  110  taps  of  the  horizontal  bank  be  correlated  with  their  seven  closest  neighbors 
of  the  vertical  bank.  Because  the  program  perfonns  all  correlation  simultaneously  (i.e.,  770  in  the 
above  example),  the  software  has  been  configured  to  segment  the  input  sequence  length  into  blocks 
that  make  machine  core  usage  realistic.  We  therefore  require  a segment  loop  that  operates  by  auto- 
matically inputting  tape  frames  to  fill  the  core  allotment.  Temporary  results  are  then  calculated  and 
the  process  repeats  itself  until  the  total  requested  sample  space  has  been  analyzed. 


4. 1 4 TAP  1 AND  0 DEPENDENCY 

Tlie  degree  of  coherence  between  a tap's  I and  Q components  is  calculated  for  the  zero-lag 
case  through  application  of  an  .ilgorithm  similar  to  that  described  in  section  4.18  Since  I and  O are 
both  real  functions,  their  normalized  cross-correlation  zero-lag  value  is  given  by 


N 
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Figure  4- 13.  Cross-Correlation  Program  High-Level  Algorithm 


This  parameter  is  given  as  a numerical  output  together  with  the  following  secondary  results: 


'rms 

^rms 


N 

(1/N)^^  I.- 

(1/N)^  0j“ 

i=l 


A 


rms 
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4,20  TAP-GAIN  AUTOCORRELATION  FUNCTION 

Prior  to  application  of  the  noise-removal  option,  the  tap-gain  autocorrelation  function,  U(t,J). 
is  calculated  as  a subset  of  the  algorithm  described  in  section  4.18.  That  is,  for  the  zero  tap  offset  the 
unnormalized  tap  correlation  function  represents  the  noise-present  component  of  the  tap-gain  auto- 
correlation function: 


U(T,f) 


I N 

= —T. 


^.i  Ki-i 


(4-54) 


where: 

T = delay  value  associated  with  the  particular  tap 
f = time-lag  variable 

X = tap  complex  process  from  either  the  horizontal  or  vertical  polarization  bank. 

As  the  program  option,  the  above  process  may  be  applied  to  the  cross-polarized  tap  banks 
(X  and  Y);  i.e.. 


UxY<i‘  i>  - TriZ 


i=l 


^.i 


(4-55 » 
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A second  program  option  allows  one  to  estimate  and  remove  the  noise  contribution  to  the 
above  tap-gain  autocorrelation  function.  The  noise  estimate  is  made  by  applying  equations  derived  in 
section  4. 14  to  a multipath-free  region  of  the  tap  bank;  i.e.. 


^ I 
S "n 


Un(J) 


i=l 


ln-'  + 


where: 

U^(^)  = estimate  of  average  complex  noise  contribution  to  the  U(t.J)  function 
Xj  = tap  process  in  multipath-free  region  of  the  correlator  bank  . 


(4-56) 


By  subtracting  Up(J)  from  U(r.J)  we  then  obtain  an  estimate  of  the  noise-free  tap-gain  auto- 
correlation function 


Unf^r.J)  = U(r.J)-U^(f)  . 


(4-57) 


4.2 1 C'OVARIANCh  AND  CROSS-SPKCTRAL  DENSITY  ESTIMATES  OF  HORIZONTAL- 

VERTICAL  CHANNEL  MULTIPATH  DATA 

The  software  described  below  assumes  that  the  tap  transposition  programs  (sec.  4.5)  have  been 
executed  for  both  the  horizontal  and  vertical  channel  multipath  tapes  and  that  the  processed  data  has 
been  placed  in  a tile  on  either  magnetic  tape  or  disk.  In  addition,  the  output  record  length  of  the  tap 
transposition  must  be  256  words  (512  packed  complex  samples).  The  number  of  records  per  segment 
must  be  less  than  or  equal  to  four. 

4.21.1  Data  Files 

The  programs  generate  eight  temporary  data  files  during  execution.  Their  structures  are 
detailed  in  figure  4-14.  HORIZ.DAT  and  VliR  [ .DAT  are  horizontal  and  vertical  channel  buffers. 

Ttie  variable  IREC  has  a value  I greater  than  the  number  of  records  per  segment  in  the  tap  transposi- 
tion output.  N I .APf  represents  the  /nimber  of  ru/»s  to  be  used  for  cross-channel  eslimates.  ( The  first 
tap  to  be  used  is.  of  course,  an  input  variable.)  Hie  last  column  in  the  file  structure  is  a pseudo-tap 
that  has  been  labeled  "carrier."  It  represents  the  response  of  the  system  to  a cw  tone  at  the  carrier 
frequency.  It  is  obtained  by  summing  a set  of  user-selected  taps  that  ostensibly  contain  all  significant 
energy.  All  tap  processes  have  gone  through  a mean-removal  operation  before  placement  in  the  file 
data  buffers. 


L 


4-44 


HORIZ.DAT.VERT.DAT  HTRAN . DAT , VTRAN . DAT 


Figure  4- 14.  H and  V Data  File  Structure 


H I RAN, DAI  and  VTRAN.DAT  LOiitain  one-record  delayed  transforms  of  the  data  in  HORIZ. 
DAI  and  VRRI'.DAT,  Kach  transfomt  is  two  transposition  records  in  length.  Intermediate  cross- 
channel covariance  estimates  are  stored  in  CRCOV.DA  f.  Intermediate  cross-spectral  density  estimates 
are  stored  in  CROSF.DAT.  All  phase  information  is  retained  to  allow  calculation  of  both  cross- 
power and  cross-phase  estimates. 

As  a result  of  the  relatively  small  amount  of  additional  work  required,  the  program  also  calcu- 
lates intertap  covariances  for  both  the  horizontal  and  vertical  channels.  These  estimates  are  located  in 
file  INCOV.DAT.  The  number  of  records  in  this  tile  is  2 x NTAPA,  where  NTAPA  represents  the 
/(umber  of  ra/zs  to  be  used  for  intra-channel  estimates.  A user-selected  pivot  tap  is  correlated  with 
each  of  the  NT  APA  taps.  Horizontal  channel  intertap  covariances  are  stored  in  the  first  set  of  NT.AP.A 
records;  vertical  channel  intertap  covariances  are  stored  in  the  second  set  of  NTAPA  records. 

DtNS.D.AT  is  temporary  storage  for  a joint  density  description  of  the  horizontal  and  vertical 
channel  pivot  tap  processes. 


4.21.2  Processing  Hxample 

All  covariance  and  spectral  estimates  are  derived  via  the  FFT.  These  techniques  are  well 
known  and  need  not  be  described  here.  However,  we  will  brietly  describe  a novel  procedure  that 
allows  both  estimation  algorithms  to  be  executed  at  the  same  time  while  maximum  use  is  made  of 
common  infonnation  requirements. 

Figure  4-1  5 contains  a record  description  of  two  complex  channel  processes.  For  simplicity, 
suppose  a statistical  description  of  only  segment  D of  the  vertical  channel  data  is  required.  .Assume 
that  lag  results  are  required  from  -5  I 2 to  -t-5  I I . Let  9 represent  the  forward  FFT  and  let  JT'*  repre- 
sent the  inverse  operation.  Defining  the  following  quantities: 

FA  = J"(A)  FD  = J^(D) 

FB  = ^(B)  FVIP  = 7(V;i 

FC  =7(0  FVIIP  =7(V;+|)  . 

then, 

NFGATIVF  I AGS=  FVIP*)J;;il  +^)L-^-'(FB  ■FVllP*)J5n  (4-,ss» 

POSITIVE  LAGS  = „L^-'(FB  -FVIP*)  , + „L^-'(FC  1 VllP’iJ^,,  (4-.S'M 
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where  Ihe  [ [ iiKlieates  that  the  first  half  of  tlie  inverse  transform  is  to  he  used.  I he  last  half 

of  the  inverse  transform  eontams  aliased  correlation  information  and,  therefore,  is  useless.  Cross- 
speetral  density  information  is  easily  derived  hy  noting  that 


H)(k)  = rVIPtkl  + (-1  i*^  I VI  IPikl, 


(4-()0l 


where  explicit  indexing  notation  has  been  introduced  lor  clarity.  I hus.  Hi)  is  easily  derived  by  a set 
of  adds  and  subtracts 


Unsmoothed  spectral  density  = FB  • FD*  . 


(4-bl  » 


.A  data  window  has  not  been  used,  to  allow  the  use  ot  common  FFT  operations  lor  both  esti- 
mation procedures.  As  a result,  it  is  necessary  to  use  a spectral  window  for  spectral  density  side-lobe 
suppression.  Since  cross-spectral  estimates  are  often  characteri/ed  by  a sinusoidal  structure,  it  is 
necessary  to  choose  a spectral  window  free  of  negative  e.xcursions  to  avoid  resonance.  For  this  reason 
a P.ARZFN  window  is  employed. 

It  is  clear  that  if  this  procedure  is  followed  for  NT.AK'  sets  of  horizontal-vertical  channel  tap 
processes,  all  the  infomiation  required  for  intertap  covariance  estimates  is  present  in  transform 
representation.  Only  one  additional  FF  T per  channel  per  record  is  required  for  each  intratap  covari- 
ance desired. 


4.:i..^  Flowchart 


.A  How  chart  for  the  basic  program  processing  loop  is  detailed  in  figure  4-lo.  I he  terms 
W,  .\.  V.  Z represent  data  storage  arrays  of  length  equivalent  to  two  tap  transposition  records.  I he 
notation  L • J denotes  that  the  first  half  of  the  data  is  to  be  extracted.  The  notation  P • “1  means 
that  the  first  half  of  the  d.ita  is  to  be  extracted  and  then  shifted  in  index  by  one  record  length.  The 
eciuation 


W = L,\J  + PN'n  . 


(4-«:) 


indicates  that  the  first  half  of  V\  is  the  same  as  the  first  half  of  \ while  the  second  half  of  VV  is  the  same 
as  Ihe  first  half  (d  A 
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Figure  4- 16.  (Continued! 
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Figure  4-16.  ( Concluded) 
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Smoothing  operations  are  not  explicitly  indicated  on  the  Ilow  chart  because  they  are  done 
outside  the  loop.  I he  tlow  chart  assumes  that  processing  of  record  i is  occurring  and  that  HTRAN. 
DAT  and  VTRAN.DAT  contain  transform  information  corresponding  to  the  record  i-1.  In  addition, 
the  file  transfer  indication 

W- RfcT|(-(HORlZ;i;i+l).  (4-63) 

means  that  records  i and  i+1  for  record  1C  in  file  HORlZ.DAT  are  to  be  transferred  to  W.  This  nota- 
tion is  adopted  for  convenience  only  since  the  data  buffering  required  make  the  fetch  operations  much 
more  involved.  The  notation 

W^REC|(^-(HORIZ:i,  0).  (4-64) 

means  that  the  second  record  transferred  to  W is  identically  zero.  Tap  1C  is  an  index  of  the  pivot  tap 
specified  by  the  user. 


4.22  DIGITAL  TAPE  EORMATS 

Three  digital  magnetic  tapes  are  created  in  the  overall  multipath  algorithm  execution  sequence. 
Referred  to  as  Dl.  D2.  and  D3  in  figure  3-1,  these  tapes  are  described  below. 

4.22.1  Dl . The  Reformatted  Digital  Tape 

This  tape  contains  the  reformatted  SACP  frame  words  together  with  merged  time  code,  frame 
count,  frame  lock  loss  count,  parity  error  count,  and  a parity  error  directory.  Alt  physical  records 
written  on  the  tape  have  a length  equal  to  I 2 of  Ihe  above  data  sets.  A data  set  is  comprised  of  2‘>2 
eight-bit  bytes  as  described  in  table  4-3. 

The  frame  counter  starts  at  0 at  the  first  frame  of  a file  and  is  increased  by  1 for  every  input 
frame  of  the  file.  A file  is  used  to  contain  Ihe  data  for  one  test  condition. 


The  lock  loss  counter  starts  at  0 and  is  incremented  by  1 for  each  frame  lock  loss  encountered 
in  tile  TEE  equipment  when  data  is  input  from  Ihe  analog  tape. 


TABLE  4 3.  FORMAT  OF  REFORMATTED  DIGITAL  TAPE,  D1 


Byte 

number 

Number  of 
words 

Bits  per 
word 

Word  description 

1-2 

1 

16 

Frame  counter 

3-4 

1 

16 

Lock  loss  count 

5-6 

1 

16 

Parity  error  count 

7-8 

1 

16 

Data  track  and  polarization 

9-12 

3 

16 

Merged  time  code 

13-22,  24 

11 

8 

Receiver  parameters 

23 

1 

8 

Parity  error  tipoff 

25-36 

12 

8 

Direct  tap  samples 

37-260 

224 

8 

Indirect  tap  samples 

261-262 

2 

8 

Subframe  data  and  index 

263-292 

31 

8 

Parity  error  directory 

The  parity  error  count  gives  the  number  of  parity  errors  (from  the  analog  input)  detected  in 
the  frame. 

The  correspondence  between  data  set  byte  number  and  SACP  receiver  parameter  word 
number  is: 


Data  set  byte  14  13  l(i  15  18  17  20  10  22  21  24 

SACP  receiver  word  1 2 3 4 5 b 7 8 0 10  11. 

Tile  SACP  receiver  words  are  de.scribed  in  section  4. 10. 

The  data  track  and  polarization  descriptors  are  encoded  as  numeric  and  ASCII  words;  they 
occupy  the  first  four  bits  of  the  7*^  and  8**'  data  set  bytes,  respectively. 

The  merged-time  code  words  are  encoded  as  shown  in  table  4-4. 

The  parity  error  tipoff  byte  is  encoded  as  0 if  no  parity  errors  were  detected  in  the  frame,  it  is 
nonzero  otherwise.  The  parity  error  directory  contains  encoded  addresses  where  parity  errors  were 
detected  one  per  Ki-bit  word.  The  addresses  are  decoded  as  tollows;  receiver  parameters  are  num- 
bered I3()|()  to  I2()|()  (ilescending);  the  taps  are  numbered  ID'io  to  I where  eacli  value  corre- 
sponds to  a Q,1  point  pair.  The  last  entry  In  the  table  is  followed  by  a wonl  conlaming  - I. 

I he  direct  tap  samples  are  in  the  order 


yi.  II,  04,  14.  02.  12.  05.  15.  03.  13.  0('.  K' 
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TABLE  4.4.  MERGED  TIME  CODE  FORMA  T 


Bytes 

15  14  13  12 

00 

CO 

O) 

o 

t 

7 6 5 4 

3 2 10 

9, 10 

Units 

BCD  seconds 

Hundreds 

Tens 

3CD  milliseconds 

Units 

Bytes 

15  14 

O 

CM 

CO 

Bit 

9 8 7 

6 5 4 3 

2 1 0 

B 

Tens 

BCD 

Units 

minutes 

Tens 

BCD  seconds 

The  indirect  tap  samples  are  in  the  order: 


01.  II.  057,  157,  02,  12 056,  156,  0112,  1112. 


Tile  subframe  data  word  is  blank  whereas  the  index  counts  from  0 to  5 and  then  recycles. 


4.22.2  D2:  The  Sj^lr.co)  Save  Tape 

Tape  D2  contains  the  noise-present  estimate  of  the  scatter  channel’s  delay-Doppler  function. 
Also  present  on  the  tape  is  a directory  that  lists  the  identity  of  the  data  on  the  files.  The  directory 
file  delineates  the  identification,  creation  date,  and  time  associated  with  all  the  tape  files.  Each 
Sjj(T,w)  array  occupies  its  own  tape  file.  The  scatter  function  array  is  accompanied  by  data  and 
polarization  test  descriptors  plus  the  receiver  and  analog  tape  status  parameters  that  occur  at  the  begin- 
ning and  end  of  the  sample  space  interval.  Under  the  CDC  6600  KRONOS  2. 1 operating  environment, 
the  following  FORTRAN  binary  write  statements  fully  describe  the  data  arrangement  of  each  Sj^tT.cj) 
file; 

WRITE  (lUHBUF(l).  1=1,20) 

WRITE  (lUMBUFd),  1=1,20).  (TlT(i),  1=1.4) 

WRITE  (lUlNFREO.  NTAP 
WRITE  (IU)(FR(I).  1=1,  N FREQ) 

WRITE  (lU)  (TNB(I).  1=1,  NTAP) 

DO  10  1=1,  NTAP 

10  WRITE  (lU)  (SE(J.I).  J=l,  NFREQ) 
lU  is  the  FOR  TRAN  output  unit  (2). 
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BUF  is  an  array  of  20  words.  These  contain  the  receiver  and  analog  tape  status  parameters 
associated  with  the  first  frame  of  the  ensemble  used  to  calculate  the  scatter  function. 

BUFl  is  an  array  of  20  words  corresponding  to  the  receiver  and  analog  tape  status  parameters 
associated  with  the  end  of  the  sample  space. 

TIT  is  a four-word  array  identifying  date  and  polarization  of  the  test  condition. 

NFREQ  is  the  number  of  frequency  points  saved  from  the  periodogram  output. 


NTAP  is  the  number  of  taps  processed. 

FR  is  the  frequency  array  in  hertz.. 

TNB  is  the  tap  number  array. 

SF  is  the  scatter  function  array  with  SF(  1.1)  corresponding  to  FR(  I ),  TNB(  I ) 

SF(2.1)  FR(2).  TNB(l) 

SF(1,2)  RF(l),  TNB(2). 


4.22.3  D3:  The  Delay-Spectra  Time  History  Storage  Tape 

Tape  D3  is  used  primarily  to  contain  the  record  elements  of  the  delay-spectra  time  history 
array.  It  is  generated  by  the  quick-look  process  algorithm,  with  each  subinterval  condition's  output 
being  written  as  a separated  physical  record.  The  record  structure  is  as  shown  in  table  4-5.  A tape 
file  is  used  to  contain  all  the  records  for  a particular  test  condition. 


4.23  AIRBORNE  SYSTEM  PARAMETERS 


In  its  normal  mode  of  operation,  the  multipath  prober  test  utilizes  wideband  PN-PSK  sequence 
transmissions  from  the  KC-135  aircraft.  Airborne  system  parameters  necessary  for  proper  interpreta- 
tion and  normalization  of  multipath  scatter  return  include  ( 1 1 modulator  mode  characteristics.  I 2)  air- 
craft pitch,  roll,  heading,  groundspeed.  and  altitude.  (3)  aircraft  location.  (4)  antenna  mode  and 
polarization,  and  ( 5)  RF  power  input  levels  at  the  antennas.  With  the  exception  of  aircraft  location, 
which  is  obtained  from  a manual  log  of  INS  waypoints,  these  data  are  recorded  over  all  Ihght  test 
intervals  during  which  scatter  channel  probing  takes  place.  In  this  section  the  recorded  parameters, 
the  recording  format,  and  the  data  conversion  are  described. 
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TABLE  4-5.  FORMAT  OF  DELAY-SPECTRA  TIME  HISTORY  TAPE.  D3 


Byte  number 

Number  of  words 

Bits  per  word 

Description 

t-24 

24 

8 

Receiver  and  analog  tape  status  data  at  start  of 
test  subinterval  (see  bytes  1-24  of  D1 1 

25-48 

24 

8 

Receiver  and  analog  tape  status  data  at  end  of 
test  subinterval 

49-272 

112 

16 

Indirect  tap  frequency-spread  in  the  order: 

Tap  1,  tap  57,  tap  2,  tap  58, ... , tap  56, 
tap  112 

273-296 

6 

32 

Direct  tap  delay  spectra  in  the  order: 
Tap  1 , tap  4,  tap  2,  tap  5,  tap  3,  tap  6 

297-742 

112 

32 

Indirect  tap  delay  spectra  in  the  order: 
Tap  1,  tap  57,  tap  2 tap  112 

4.23.1  RecorcleJ  Parameters  und  Recording  F-ormat 

.Airborne  system  parameters  are  recorded  on  computer-compatible  800-bpi.  seven-track  tape 
(six  bits  plus  parity  I using  a Kennedy  8707  recorder.  Byte  size  is  six  bits.  Data  is  segregated  into 
records  containing  484  bytes.  The  record  format  is  given  in  table  4-b. 

.Apart  from  the  time  code  words,  the  sixth  bit  of  each  data  byte  contains  no  infomiation. 
Thus,  the  power  monitor  data  and  aircraft  parameters  are  digitized  at  a resolution  level  of  10  bits. 


4.23.2  Data  Conversion  and  Processing 

Figure  4-1  7 illustrates  the  basic  CDC  (>600  softwaie  algorithm  processing  steps  associated  with 
reduction  and  analysis  of  the  airborne  system  parameters  tape.  I'he  first  step  involves  unpacking, 
decommutation,  and  processing  of  the  time-multiplexed  calibration  signals,  which  in  general  will 
precede  the  airborne  system  data.  The  calibration  file  is  analyzed  to  yield  upper  and  lower  calibration 
levels  for  the  appropriate  channels.  I'he  signal  data  is  subsequently  processed  according  to  the  follow- 
ing paragraphs.  (The  output  description  is  given  in  section  4.24.) 
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TABLE  4-6.  AIRCRAFT  PARAMETERS  RECORD  FORMA  T 


Bytes 

Description 

0-3 

IRIG  B time  code 

4 

SACP  modulator  mode 

5 

RF  subsystem  mode 

6,7 

Spare  (no  data) 

8,9 

Internal  power  monitor,  direct  channel  PA 

10,11 

Power  meter,  direct  channel  (HP  435A) 

12,13 

Internal  power  monitor,  V channel  PA 

14,15 

Power  meter,  V channel  (HP  435A) 

16,17 

DC  power  supply,  1(X)-W  PA 

18,19 

Internal  power  monitor,  H channel  PA 

20,21 

Power  meter,  H channel  (HP  435A) 

22,23 

Aircraft  altitude  (feet) 

24,25 

Aircraft  coarse  heading  (degrees) 

26,27 

Aircraft  fine  heading  (degrees) 

28,29 

Aircraft  roll  angle  (degrees) 

30,31 

Aircraft  pitch  angle  (degrees) 

32,33 

Aircraft  ground  speed  (knots) 

34,35 

Spare  (no  data) 

36-483 

1 5 additional  scans  of  parameters  contained  in  bytes  6 to  35  with  the 
last  scan  truncated  at  aircraft  groundspeed 

4.23.3  SACP  Modulator  Mode  Code 

The  five  data  bits  contained  in  the  fourth  byte  of  each  record  are  used  to  designate  the  SACP 
modulator  chip  rate  (bits  2,3.4)  and  code  length  (bits  O.l  ).  Switches  that  set  these  parameters  provide 
coded  bits  to  the  recording  system  according  to  table  4-7. 

4.23.4  RF  Subsystem  Mode 

The  aircraft  antenna  and  the  nominal  polarization  over  which  SA('P  transmissions  take  place 
are  delineated  by  the  four  least  sigi.ificant  bits  of  the  fifth  record  byte.  Coding  is  shown  in  table  4-S. 


4-.SS 


RUN  INIT. 


CALIBRATION 

COEFFICIENTS 

FILE 


JNPUT  CASE  DATA 


CASE  LOOP  {TIME  SEGMENT) 


TERMINATE 

RUN 


Figure  4-17.  Airborne  Tape  Processing  - High-Level  Program  Algorithm 


TABLE  4-7.  SACP  MODULATOR  MODE  CODE 


Bit  pattern  and 
decimal  no.  equiv 

Code  length 

001  = 1 

255 

01  = 1 

010=  2 

511 

10  = 2 

on  = 3 

1023 

11=3 

100  = 4 

101  = 5 

110  = 6 

TABLE  4-8.  ANTENNA  TYPE  AND  POLARIZATION  CODE 


Antenna  type  (bits  1,2) 

Bit  pattern  and  decimal  no.  equivalent 

Operational  (slot-dipole) 

00  = 0 

Side  multipath  (SMP) 

01  = 1 

Front  multipath  (FMP) 

10  = 2 

Polarization  (bits  3,4) 

Bit  pattern  and  decimal  no.  equivalent 

Vertical/horizontal 

00  = 0 

Left-hand  circular 

01  = 1 

Right-hand  circular 

10  = 2 

I 

I 

4.23.5  Power  Monitors 

With  the  exception  of  the  three  power  monitors  internul  to  the  RF  umplifiers  used  lor  the 
direct,  hori^ontul.  and  vertical  channels,  the  recorded  power  monitor  values  from  tlie  external  RF' 
power  meters  are  converted  directly  to  a voltage  level  and  output  without  further  reduction,  i.e.,  out- 
put voltage  = (word  value,' 2“  ) x 5.0  V.  The  internal  power  monitor  data  are  converted  to  their  respec- 
tive channel  power  amplifier  readings  through  the  following  linear  relationship: 

PA  output  = (\/U('l.A.  walls, 

where: 

X = data  word  value 

UC  = upper  calibration  signal  (recorded  on  the  first  tape  file) 

A = power  amplifier  output  (as  measured  when  calibration  tile  is  created,  input  to  routine 

via  data  card) 
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4.23.6  .Aircraft  Paraiticiers 


.Aircraft  groundspeed  is  derived  through  the  following  relationship; 


Groundspeed  = (X/5  1 2)  (500).  knots. 


where  x is  the  data  word  value. 


All  the  other  aircraft  parameters  are  obtained  through  use  of  the  following  relationship: 


A/C  parameter 


X-  LC 
UC  - LC  ' 


(Ai  — At)  + A- 


where: 

X 


X = data  word  value 

UC  = upper  calibration  voltage  level 

LC  = lower  calibration  voltage  level 

A I = parameter  value  corresponding  to  UC 

At  = parameter  value  corresponding  to  LC  . 


The  UC  and  LC  values  are  obtained  from  the  tape  file  containing  calibration  readings.  The  Aj  and 
At  values  are  set  internally  in  the  computer  program  and  have  been  derived  from  information  sup- 
plied by  NAFEC. 

As  previoulsy  noted,  tine  and  coarse  headings  are  components  of  the  aircraft  system  parameters 
data.  Fine  heading  varies  from  0°  to  1 0°,  whereas  coarse  heading  is  given  over  the  full  360°  range.  In 
general,  the  aircraft  actual  heading  is  obtained  by  rounding  down  the  coarse  reading  to  the  nearest  10° 
and  then  adding  the  fine  heading.  Due  to  calibration  uncertainties,  this  procedure  may  cause  large 
errors  if  the  actual  heading  is  close  to  a mod- 10°  value  (e.g..  if  coarse  heading  reads  321  and  fine  head- 
ing reads  ‘f,  then  the  actual  heading  is  undoubtably  closer  to  3 I6°  than  329°).  To  circumvent  such 
potential  problems,  the  following  logic  may  be  applied  to  the  output  to  derive  the  aircraft  heading: 


= CH^j-  10  + FTI 

for  A+2  < FM 

= ( "ru  ^ m 

forA-2  > FH 

= ^ ”rd  ^ " 

(7lherwisc . 
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where: 


CH 

CHr 

FH 

A 


rd 


= rounded-down  (nearest  10^)  coarse  heading 
= rounded-up  (nearest  lO*^)  coarse  heading 
= fine  heading 


= CH-CH 


rd 


4,24  SAMPLE  MULTIPATH  ANALYSIS  OUTPUT 

In  this  section,  sample  numerical  and  plotted  outputs  are  illustrated  and  described  for  each  of 
the  previously  presented  algorithms. 


4.24.1  Quick-Look  Analysis 

Immediately  following  the  processing  of  a quick-look  test  condition,  the  results  are  displayed 
in  graphic  and  tabular  form  on  the  Tektroni.x  401  2 display  and  Gould  4800  printer,  respectively.  An 
example  of  the  numerical  hard-copy  output  is  given  in  table  4-4,  The  top  line  of  the  output  delineates 
the  analog  tape  description:  name  (e.g.,  AT  1 27).  file,  and  record  number  from  which  the  data  was 
dumped.  We  then  have  two  rows  of  data  words  corresponding  to  the  analog  tape  and  SACP  receiver 
parameters  (see  sec.  4.10  for  the  format  description)  as  they  existed  at  the  beginning  (start)  and  end 
(stop)  of  the  subinterval  sample  space.  Data  arrays  for  the  indirect  tap  frequency  spread  (titled 
DOPPLER),  the  direct  delay  psd  (titled  DIRECT  DELAY),  and  the  indirect  tap  bank  delay  psd 
(titled  DELAY)  follow  the  start/stop  receiver  parameters.  The  direct  delay  taps  are  interleaved  and 
thus  the  row  array  represents  tap  1,4.2.  5.  3.  b.  The  indirect  tap  Doppler  and  delay  data  are  arranged 
in  order,  reading  from  left  to  right  on  a row-by-row  basis  (each  row  is  preceded  by  the  tap  number 
associated  with  the  adjacent  data  value).  As  indicated  in  section  4.1 . the  quick-look  results  are 
unnormalized  for  the  size  of  the  sample  space.  For  the  example  output  presented,  a 2-sec  sample 
space  was  processed  and  the  sampling  rate  was  approximately  610  sec'* . Thus,  for  nonnalization.  the 
frequency-spread  data  is  divided  by  2 and  the  delay  spectra  data  by  30.S  (only  every  fourth  point  is 
processed  for  delay-spectra  calculation). 

The  401  2 quick-look  display  presents  a graph  of  the  numerical  data  arrays  described  above. 

The  fields  (e.g..  DOPPLER,  DIREC1  DELAY,  and  DEL.A^’)  are  separately  and  automatically  scaled  to 
improve  detail  and  resolution.  Up  to  five  subinterval  segments  may  be  displayed  at  one  time.  The 
data  field  arrangement  is  shown  in  figure  4-1 K.  For  this  particular  example,  the  indirect  taps  are 
occupied  not  only  by  multipath  return  but  also  by  the  direct  signal  from  the  multipath  proberantenna. 
located  roughly  at  tap  10. 
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TABLE  4-9.  QUICK-LOOK  NUMERICAL  OUTPUT 


tape;  at  127  DUMP  CF  FILE  Oi  RECCRD  340 
START 


JS644  0 

STOP 

0 

44903 

10T>64 

£1112 

20020  -40 -4003 

4144 

20111 

304 

55150  0 

0 

44003 

30561 

21112 

20020  -4057  -4003 

4144 

21111 

304 

DOPPLER 

1 594  £.11  £.?4  604  646  596  646 


0 < 

626 

650 

645 

652 

652 

617 

623 

15 

656 

645 

637 

651 

666 

633 

650 

22 

634 

628 

636 

622 

625 

649 

634 

29 

642 

633 

664 

644 

645 

612 

642 

36 

646 

613 

316 

3.90 

403 

453 

400 

43 

545 

554 

571 

593 

646 

622 

610 

50 

598 

607 

631 

620 

624 

643 

596 

57 

613 

644 

603 

592 

676 

647 

622 

64 

602 

641 

633 

609 

626 

647 

625 

71 

639 

622 

629 

593 

644 

610 

633 

78 

642 

646 

621 

665 

633 

633 

661 

05 

655 

631 

663 

637 

629 

615 

657 

92 

620 

624 

623 

609 

622 

621 

667 

99 

642 

643 

600 

622 

633 

657 

634 

106 

647 

660 

6o7 

613 

6.36 

655 

617 

DIRECT  DELAY 

£344  4904  1£11969  £025  1294765  1076 


DELAY 


1 

37153 

35532 

37351 

34181 

35575 

32679 

38228 

8 

34824 

36381 

32  309 

34641 

36957 

33283 

35553 

15 

37304 

35220 

35.200 

33244 

33315 

34777 

36907 

22 

37791 

38298 

36363 

37679 

3477£ 

36362 

39355 

29 

35011 

35059 

32871 

37131 

39159 

41732 

33432 

34091 

152195 

14i303 

15  S5253  77AcO  •&c'04 

43 

52515 

4£029 

459.59 

55-M6 

4.1.137 

390.^6 

50 

57 

41254 

39261 

35420 

35475 

3*9624 

402*92 

39543 

57 

39154 

35862 

33564 

36346 

40120 

38077 

38485 

64 

35543 

36028 

37567 

36430 

355*55 

33199 

•53191 

71 

35234 

4o;..j0 

0.i,:ig3 

3374? 

35^»?  0 ? 

30427 

?tl06 

78 

35164 

33092 

3I.-09 

33341 

38544 

36661 

34155 

85 

35720 

35806 

35238 

35693 

53033 

39445 

36643 

52 

32413 

f 

33?.?' 6 

■35263 

37750 

i 

34342 

99 

3."’445 

St-OcO 

32386 

32321 

35047 

33755 

36339 

106  38203  36014 

. ::S563 

3598/ 

’ 34623  203159  35307 
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4.24.2  Delay-Spectra  Time  History 

Three-dimensional  plots  with  hidden  lines  removed  and  numerical  output  representations  are 
provided  for  the  delay-spectra  time  history  distribution.  Figure  4-19  represents  a sample  plotted  out- 
put for  a overland  CONUS  data  run.  The  plot  title  contains  information  pertaining  to  test  data, 
polarization,  tap  resolution  or  width,  and  the  time  interval  over  which  the  time  history  is  plotted. 
Typically,  the  subintervals  over  which  a delay  spectrum  is  calculated  are  on  the  order  of  2 sec.  and  60 
distributions  per  graph  are  usually  plotted.  Thus,  the  time  interval  subtended  spans  roughly  2 min. 
The  grid  is  labeled  according  to  the  (x.y.z)  coordinate  system,  with  axis  descriptors  given  at  the 
bottom  of  the  output  page.  This  form  of  axis  description  is  used  to  facilitate  the  labeling  software 
required  to  support  the  option  available  in  the  3-D  program  that  allows  for  the  arbitrary  choice  of  any 
desired  point  of  view. 

Corresponding  to  each  subinterval  delay  spectrum  is  a set  ol  numerical  output  containing  the 
following  information: 

a.  Receiver  and  analog  tape  status  parameters  corresponding  to  conditions  existing  at  the 
beginning  and  end  of  the  subinterval  over  which  the  following  parameters  are  determined 
(see  following  paragraph  for  description ) 

b.  The  frequency-spread  estimate  associated  with  each  of  the  indirect  channel  taps 

c.  The  direct-tap  bank  delay  spectra 

d.  The  indirect-tap  delay  spectra 

e.  The  total  rms  energies  in  the  direct  and  multipath  tap  banks  (signal  plus  noise)  . 

An  example  of  a numerical  output  is  presented  in  table  4-10. 


4.24.3  Receiver  and  Analog  Tape  Parameters 

A typical  output  from  the  algorithm  that  analyzes  the  SACP  receiver  and  analog  tape  param- 
eters is  given  in  table  4-10.  The  output  begins  with  a descriptor  indicating  the  track  designation  from 
which  the  analog  data  was  extracted  Following  are  two  rows  of  parameters  delineating  the  receiver 
status  descriptors,  the  time  merged  from  the  IRKi-A  time  code  track,  the  frame  counter,  the  frame 
lock  loss  counter,  and  the  parity  error  counter.  The  counters  are  referenced  to  the  first  frame  pro- 
cessed for  a given  test  condition  on  the  PDP  I I 4.3. 
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4.24.4  Mean  Square  PireLt  t'nergy 

I'lie  mean  square  direet  energ>  Is  ealeulated  willim  the  program  loop  that  performs  the  Sit.coI 
ealeulation.  A numerieal  output  from  this  algorithm  is  provided  at  the  end  of  the  total  ensemble  eal- 
eulation.  .An  example  of  such  an  output  is  included  m section  4.24.7. 


4.24.5  Delay-Spectra  Time-Domain  Analysis 


As  mentioned  in  section  4.4,  the  delay  spectra  as  determined  from  a time-domain  analysis  of 
the  complex  tap  outputs  are  calculated  to  serve  as  an  integrity  check  on  the  delay-Doppler  function 
and  associated  integral  operations.  The  output  from  this  algorithm  is  given  numerically  and  is  treated 
further  in  section  4. 24.'. 


4.24.6  Delay-Doppler  Scatter  Function 


Identical  output  formats  are  used  for  the  noise-present  and  noise-free  estimates  of  the  scatter 
channel's  delay-Doppler  function.  Numerical  and  three-dimensional  plotted  outputs  (with  or  without 
the  removal  of  hidden  lines)  are  used  to  represent  this  function.  1-xamples  of  numerical  output 
segments  are  given  in  table  4-1 1 (noise  present)  and  table  4-12  (noise  removed).  The  page  begins  with 
a heading  descriptor  delineating  test  date,  polarization,  and  whether  the  estimate  is  noise  free.  F-ollow- 
ing  is  the  S(r.  cc)  array  presented  in  matrix  form.  The  left-hand  column  corresponds  to  (he  frequency 
in  hertz  of  the  smoothed  periodogram  output  I'his  is  given  once  every  Wi,,Af  hertz,  where  W,/.  is 
the  half  width  of  the  spectra  smoothing  window  and  Af  is  the  unsmoothed  periodogram  output 
sampling  rate  that  is  equal  to  the  frame  rate  divided  by  the  number  of  complex  points  used  for  the 
FFT  input.  With  the  exception  of  the  frequency  delineator,  the  top  row  of  the  matrix  corresponds  to 
the  delay  tap  number  associated  witii  the  PSD  (linear  relative  watts)  given  in  the  below  column.  From 
these  figures  it  is  seen  that  the  form  of  the  numerical  output  gives  an  easily  recognized  display  of  the 
S(T,cc)  energy  concentration. 

Figures  4-20  and  4-21  represent  examples  of  the  hidden-line-removed,  three-dimensionally 
plotted  outputs  of  S(T.cc)  for  the  noise-present  and  noise-removed  cased,  respectively.  The  hidden- 
line-present  plot  for  the  noise-removed  case  is  given  in  figure  4-22.  Fach  plot  contains  a title  describ- 
ing the  test  date,  polarization,  tap  width  or  resolution,  time  interval  over  which  the  scatter  function 
was  calculated,  and  whether  the  plot  corresponds  to  a noise-removed  or  noise-present  estimate.  The 
grid  is  labeled  according  to  the  (x.y.z)  coordinate  systems,  with  axis  labels  given  at  the  bottom  of  the 
output  page  This  form  of  axis  description  is  used  to  facilitate  the  I ’heling  software  required  to  sup- 
port the  option  available  in  the  ,TD  program  (hat  allows  for  the  arbitrary'  choice  of  any  desired  point 
of  view. 
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TABLE  4-11.  DELA  Y-DOPPLER  SCATTER  FUNCTION  NUMERICAL  OUTPUT 
(NOISE  PRESENT) 


TABLE  4-12.  DELA  Y DOPPLER  SCATTER  FUNCTION  NUMERICAL  OUTPUT 
(NOISE  REMOVED) 

ct  Zf  5 .-1-  ‘OlSt 


on  Av  I . 4.  H(f  •user  . ^ 
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300.00  300.00  100. OC 


Figure  4-20.  Three-Dimensional  Plot  of  Delay-Doppler  Function  (Noise  Present) 


The  intended  purpose  ol  plotting  the  Str.cu)  (unctions  is  twofold:  ( 1 ) to  give  the  observer  a 
quick  and  concise  description  ofthe  energy  distribution  characteristics  of  the  Str.cu)  spectra  and  (2) 
to  provide  a visual  directory  from  which  one  may  access  the  numerical  results  to  key  in  on  areas  ol 
interest  or  anomolous  behavior  of  the  scattering  function. 

4.24.7  Integral  and  Fourier  Operations  on  Str.co) 

Both  the  noise-present  and  noise-free  estimates  of  the  S(t,w)  function  are  subjected  to  opera- 
tions that  yield  the  scatter  channel’s  delay  spectra,  Doppler  spectra,  frequency  autocorrelation  func- 
tion, time  autocorrelation  function,  and  total  rms  energy.  For  the  delay  and  Doppler  spectra,  numeri- 
cal and  plotted  output  are  given,  as  shown  in  figure  4-23  and  4-24  and  table  4-13.  As  previously  men- 
tioned in  section  4.14.  plotted  spectra  for  both  the  noise-present  and  noise-removed  algorithm  outputs 
are  presented.  This  is  done  to  give  an  indication  of  the  effectiveness  of  the  ND  and  R operations. 

.Also  included  in  the  numerical  outputs  of  table  4-13  are:  1 1 ) the  delay  spectra  trom  input  (or  time- 
domain  analysis),  (2)  the  mean  square  total  indirect  energy  (integration  over  delay  variable  = T.  inte- 
gration over  Doppler  variable  = F).  (3)  the  mean  square  total  indirect  energy  for  the  noise-present  case, 
and  (4)  the  composite  mean  square  energy  of  the  direct  signal.  The  T and  F integrations  are  identical 
as  would  be  expected  and  they  are  substantially  less  than  the  signal-plus-noise  integration.  The  direct 
and  multipath  mean  square  energies  are  not  normalized  to  remove  differences  in  antenna  ERP.  receiver 
channel  gains,  etc.  This  step  is  performed  in  the  data  interpretation  phase  of  the  analysis,  which  is  not 
a DRandA  software  component. 

Numerical  and  plotted  output  formats  are  used  to  represent  the  multipath  R(^f2)  function. 

An  example  illustrating  the  distribution  of  |R(J,f2)l  is  given  in  the  three-dimensional  plot  of  figure 
4-25.  To  enhance  the  fidelity  of  the  output,  the  plot  routine  has  in  this  case  been  operated  in  a mode 
whereby  both  the  constant  { and  constant  contours,  which  occur  at  the  discrete  periodogram 
output  values,  are  drawn  with  hidden  lines  removed. 

The  frequency  and  time  autocorrelation  functions  are  given  in  the  form  of  numerical  and 
plotted  output.  For  both  numerical  outputs,  the  amplitude  and  phase  are  tabulated  relative  to  the  fre- 
quency (hertz)  or  time  (seconds)  variable,  respectively.  Numerical  examples  of  the  RlO.H)  and 
R(^,0)  distributions  are  given  in  tables  4-14  and  4-15.  Sample  plotted  results  of  the  autocorrelation 
function  magnitudes  arc  presented  in  figure  4-26. 

rile  frequency  and  time  autocorrelation  functions  are  given  only  in  the  form  of  numerical  out- 
put. For  both  distributions,  the  voltage  amplitude  and  phase  in  degrees  are  tabulated  relative  to  the 
frequency  (hertz)  or  time  (seconds)  variable.  Examples  of  the  algorithm  outputs  for  the  frequency 
and  time  autocorrelation  functions  are  given  in  figures  4-25  and  4-26, 


4-74 


I 

i 


dclay  and  dopplcr  spectra 

OATC  2/i7/n  U M C POLARI2ATION 

TAP  RCSOLUTICN  .2  MICROSCC 

IIHC  interval  IO/44/S8.9GS  TO  10/49/  4.838 
NOISE  PRESENT 


Figure  4-23.  Delay  and  Doppler  Spectra  Plot  (Noise  Present) 
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Figure  4-24.  Delay  and  Doppler  Spectra  Plot  (Noise  Removed) 


n 

n 


'■dJiW>*i^^Uj  iJ  ’./  iaJ  m 
N.'T<NJK> 

^0  •*•£>■»  9 •>»  0 9 0 9^ 


■ikj  uj  ^ <ij  uj  ^ 
^t\lt^J>yM\f\4  9 


^ Ui  W 1^  ^ 

M »vj»^-»>w%  n,T*«»*  ««^.nM* 

(\JfM**tr>J«  9 J -I  ■<> 


•«J  U>  MJ  W <4i  '«l  '■><  ‘y  ^ UJ  W ■ J 
am 

■N  0 O*^9»4^t0  T rs  f\t  ^ 


■a 

ykAa^ift^  9K. 
<M  «4  I * ^ ^ 'VJ 

I I 


UJ  Ui  'U  y Ui  'aj  aJ  O 
H^«^lfX*44.9i0 


>al  oJ  '«i  Ui  U1  UJ  U ‘ij  '«J  Ll  U -U 
«4^.0^  a«4  t 0«^mK.  «« 
INJ 

•*yi  n r 9 ♦-<•»»  > 


t K.  ^ 

wtAoous 
A an*Ai^#o*>*<a^m 

(U(VAiA<*lA^m*0*0<n^g 


uiiyujt^tyuiujtyuiyfu^iy 
r«.*^-y<>4w'V(MAaAJ7«KtK.M\ 
OtA'Mo«A«^oiAKiA*^nr' 
O>e^o  *0^O'n*t*AAj,.( 


t\t  "\t  m P j\  9 9 "n  f 


I4J  •■•  UJ  W «>l  I4J  lU  III  lal  lit  •••  til  III 

9 «Ay^«4a9ia*A^lfVfM«sl(A 
>^n>4'^.4,0«4f«.i0«9iAlt>k 
9^9 


aasaaov^aitr 


UJ  ^ 1^  to  UJ  lU  ill  III  UJ  UJ  III  III  UJ  III  III  U>  UJ  10  Ut  'iJ  UJ  ••■  III  l.l  I.I  i.i 

0>  9 WkAlA  ^^>Oif\ffi‘999‘T9  99f\i 

uNiu>».4e«K.Ajtcro«> 


<0^A«>aiuiu  ^ 9 9 49  >X>  : 


3AJ»UM(A  9 9 ^ no  lA  ■A  VI  *«i 


i^9^or>i‘^a9'990>n^9 

9\99*^'^w-t9o'^  9 O^cafT 

m«\iiMtVtA»4j)  9 ,»«^*»»Mv*rj 


*49*KilAn»49»«.t0Ut«49tU 

^•••4<M*^^l,AlA.0IU«99>« 


OfujufUJUJUJUiujuityujUiUJui 

U>AJIA^AJs0«iAlAW>V*me^ 

a 

4AtMAJAJ>o«4|0.0U»fU(0«0tMfV 


HI  ««  •«  M *0t# 


TABLE 4^  14.  FREQUENCY  AUTOCORRELA  TtON  FUNCTION  OUTPUT 


MCOUENCV  4UT0C0K«EL*TinN 


FREOU€NCT  ITU0€  ( VQlT ) PHA$t<DE6) 


• 0 

l.OOOOE^^O 

• 0 

3.90b3E«^04 

6.9289E-01 

7.732»F*0l 

7,8129C^0<» 

7.4989E-01 

I«4^62E«02 

1,1719E*08 

6.9479E-91 

-l.9d49E»02 

1•96^9E♦09 

9.6773E»01 

•S.699bP«01 

1.99)1E*09 

9.2336E-01 

-2.6ei4E«01 

2.3430E«O& 

4.6993P-01 

3.5192E«01 

2.7144»E*05 

4, 2903F-01 

9.V93«iE«0l 

3.1290€«09 

3,91’6E-0l 

1.9473E*02 

3»919<!>E«09 

3.7t>03E-Oi 

-i.4710'*02 

3.90#»3E*03 

3.0641E-gI 

-8.9006E*01 

4*2989E*09 

3.2576C-'^l 

•2.3696E*01 

4 • 6P  7sf ♦«  5 

2.9879F-01 

3. 1310^ *01 

8.*>92«*0l 

*0*> 

2,«o7ft6-01 

l.«l«»lfc*02 

S.®5>)4E*03 

2. 

-1  ,«»979C*0? 

6*290JF*''9 

2.3354e-01 

-9.  U73F*01 

2.iH39e.Ci 

•3.06»eF»Jl 

3c  *0^ 

2.17«?F-Ol 

1, 9692^*01 

7.«.21  9t  *0* 

2.l*‘llt-0l 

b.3202F*Ol 

7,»il?5r*09 

i.  v««^0t-0i 

A.  3>6<»?*''> 

0.2O3iE*O9 

?.ns32P-'>l 

-l.»>’7oe*02 

1 • e»<73P-0i 

-l.^V62P*02 

1,R046C-Cl 

-7.03A9F4w1 

9,3750f  ♦Ot> 

1 . 70bPf-wl 

1 ,b33l,F*Jl 

9.  76>oF 

7.21 1 bP*Ol 

t .Olifrc  *00 

4.295i«f  ♦02 

-4.7^''3t*02 

i.091gE^0»‘ 

1 , ^92bt-0l 

i.l32«f 

Ul7l9EM)b 

2.2^*»9«^G0 

i.?o77t-ai 

l.2900f *06 

l. 

l.2277e*^2 

♦5o 

1 , lOb6C-0l 

-1.  77i>F^‘'7 

1 , *Cft 

-1 .23'*9t*''» 

l*3672E*0*< 

l,3iPiC.A| 

-S.37b7F*^l 

1. 2c»o  U-'>* 

-1 

' . 3*^^2'-01 

1 . 1 t'2*'f-«U 

1.7<**>OF^02 

i,^POOF-0l 

1 

-4.  3<*92f  ♦02 

1 •63l6P*0© 

1 • 6<*u6 f *wO 

•A.o?30F^^;1 

l.4707F*Of 

\ .r9?ot-oi 

\.7i-iaf  *o*> 

,^,a7<»gr.«)> 

1. 0')69k*''2 

1. 7578E  *06 

4.t 

1 ,^Oo2F^n‘» 

1. 79A9E *06 

l.ir/jF-oi 

-J.37i'jF^02 

9. 3l'<*F-a2 

“7.6?b0c*0l 

l.07>Oc»O*- 

1 .03/0E-O4 

-2.l?93F*«rfl 

1.914U*^6 

7.>7^0F-02 

9.i*74F^Jl 

l.9»3U*06 

1 .02^9^-01 

9.bl09€»0l 

i.0922E*06 

«,7A94E-02 

1.*34T«^J2 

2.031 3E*0b 

7.64l9^-g2 

-l.77\3^*02 

2.n?03f *00 

“0. 72«.9e^rt1 

2. lOOAF^Oh 

». 2230E-02 

-4.4329E*0\ 

2.14a4€*06 

O.S410E-02 

2.?»32F^01 

2.t673E»06 

7.0714F-02 

7.»299E*0l 

2.2266E»0«) 

9.bl22F*02 

l.3iP6E«02 

2.2&96E706 

S.2792E-02 

-I.6069E«02 

2.3047E^06 

9.2A67F-U2 

-I.IU1F*02 

2.7438E*0fr 

9.981VE-0? 

-3.8699E*01 

2.3a26E«06 

8.9710E-U2 

9.i60SE700 

2.9219E«Uf 

7.b7SftF-02 

8.8298E»0l 

2.4609E*06 

6.9909E-02 

1. 16«2€«02 
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TABLE  4-15.  TIME  AUTOCORRELATION  FUNCTION  OUTPUT 


Tint  AUTOCTROPLAT  rUN  flJNCTJON 

TINc  AnPLITUOE (VOLT)  PHASEtOft) 


.0 

1 .n'^OOP  *00 

.3 

1.1115E-03 

». 7^67E-01 

1.  38ISOE+OZ 

2.6230E-U3 

8.6992E-01 

-8.'*75iE*3l 

3.93AAE-03 

8.0ZZ3E-01 

<..9337E*01 

3.2<.59E-03 

3.37ASP-01 

6.5S7AE-03 

Z.18Z0E-01 

-V.SZilE+Ol 

?.Hfc39E-03 

1 « b il  8t  ~0l 

8.331ZE*0T 

9. 1R03E-03 

l.nj58t-01 

-l.5035F*n? 

1 ,0‘.9ZE-02 

9.78‘'3E-0Z 

-Z.3VZ5E+01 

1.1803E-02 

Z.O13PF-02 

1.3U5E-02 

2.7763E-0Z 

-i.<.Z5<tE*nZ 

1.AA26E-02 

2.8111E-0Z 

A.'fS'.aE  + OO 

1.5739E-02 

8.<'5b?«-03 

i..08TZE*i»<: 

1.70<»9E-02 

Z.AO22F-0Z 

1. 7679E*0Z 

1.83felE-02 

T .3893F-0Z 

-Z.bbZlEAOl 

1.967ZE-02 

l.Z693t-0Z 

i.9928E+ol 

Z .098AE-0Z 

Z.8Z93E-0Z 

-1.p8ZOEaOZ 

2.ZZ95E-0Z 

1 . 1877E-0Z 

-7.Z030E+01 

Z.3607E-0Z 

1.1919E-0Z 

Z.bOA^FAOl 

Z.A918E-0Z 

l.l 376E-0Z 

-8.1A37F*nl 

Z.6Z30E-0Z 

5.b51ZE-03 

5.0590EA01 

Z.75A1E-02 

1.A7Z7E-0Z 

i.6iAAE«0Z 

z.eeszF-oz 

l.<.709e_0Z 

-<..1583E*01 

3.016AE-02 

3.933ZE-V3 

-Z.7095E-01 

3.1,‘.T5E-0Z 

1.<.Z<.1F-0Z 

-1.A381E+0Z 

3.Z797E-0Z 

1 .ZIZAF-OZ 

-Z.856ZF-01 

3.A096E-OZ 

A.7b8<tE-03 

1.9il5E*0’ 

3.»<»10E-0Z 

7.8*73F-03 

-7.135jE*01 

3.67Z1E-0Z 

l.«53<>E-03 

t).Z791E*01 

3.8033E-0Z 

3.10<.5F-03 

A.0875E+01 

3.03AAf-OZ 

b.3797£-03 

7.29<.7F*00 

7.0656F-02 

A.AA8AF.03 

-5.S80aE+0I 
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4.24.8  Spread  Calculations 


In  the  detailed  data  analysis  routine,  the  multipath  channel's  delay  spectra,  Doppler  spectra, 
freqency  autocorrelation,  and  time  autocorrelation  functions  are  processed  by  the  spread  calculation 
software.  This  routine  assumes  a density  function  of  the  form  y = f(x).  It  has  as  its  output  numerical 
calculations  of: 


a.  Vmax  corresponding  x^^^^ 

b.  The  value  of  y„,ax/“  corresponding  x-axis  intersections  (the  first  ones  to  occur  on 

both  sides  of  x^^^)  and  the  separation  of  the  x-axis  intersections  (DX) 

c.  The  same  value  as  for  item  b above  except  for  the  1/e  spread  case,  i.e., 

d.  The  first  moment  of  the  distribution 

e.  The  second  moment  of  the  distribution  as  calculated  about  the  mean  and  about  the  x^^^ 
value . 


For  each  particular  distribution  the  unit  of  the  x-parameter  output  from  the  spread  calculation 
is  to  be  interpreted  as  follows: 


Delay: 

Doppler: 

Frequency  autocorrelation: 
Time  autocorrelation: 


Tap  number 
Frequency  in  hertz 
Frequency  in  hertz 
Time  in  seconds . 


A sample  output  for  the  spread  or  moment  calculation  on  a Doppler  spectrum  is  given  in 
table  4-16. 


TABLE  4-16.  MOMENT  COMPUTATION  OUTPUT  (DOPPLER  SPECTRUM) 

XOKENT  CQXXUTtrinx  ON  nO'PLtN  SPFCTIa 


> or  Nti.  rN>i,TNtx/2.TN></c 

INTEN.  tr  .S  ■><>  lirvEI.>>l>>2>D>l 

INTia.  l./E  Nl«>  (LEVEL. XI, 12. 0«l  N.JLIOE 

INTEa.  IT  1/10  NXI,  (LtVEl.xl.x2.0XI  2.S3V1F 

fiasT  honfnt  a.iE(i2>.oo 


V.SI91F.02  1.2ta»E«02  9.JV10E.01 


♦ Oi 

6.08S1E*01 

♦01 

•6.66»7E«01 

l.OOTIf *02 

SECOND  NoacNT  asLativE  ro  n'in  <no  nix 


s.aiirF.si 


».«S20E.OI 
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4.24.9  Noise  Determination  and  Removal 

In  determining  the  characteristics  of  the  SACP  receiver  noise,  the  algorithm  calculates: 

a.  N(f)  - the  total  noise  power  spectral  density  in  a tap  region  that  precedes  the  specular 
point  return  and  is  thus  considered  to  be  multipath  free 

b.  AgMlf)  — an  estimate  of  the  thermal  noise  psd  content  of  the  total  noise  power  in  the 
multipath-free  region 

c.  Rtf)  the  difference  between  items  a and  b.  which  is  referred  to  as  the  residual  noise  con- 
tent and  is  assumed  to  be  a representative  value  for  any  given  indirect  tap 

d.  Ap(r)  - the  thermal  noise  content  for  each  tap 

e.  The  mean  value  of  the  I and  O voltage  inputs  for  each  tap 

t.  FR( n and  M(f)  the  filter  frequency  response  with  no  aliasing  and  the  filter  response 
including  the  first  four  aliased  components,  respectively. 

With  the  exception  of  the  FR(D  spectra,  all  the  above  parameters  are  provided  as  numerical 
outputs.  Samples  are  given  in  tables  4-1  7 and  4-18.  Parameters  M(f),  R(f).  and  N(f)  are  presented  is 
frequency  (hertz)  and  psd  point  pairs.  The  1 and  Q voltage  means  are  given  as  tap  number  and  mean 
voltage.  The  array  Aq(t)  is  given  as  tap  number  and  psd  point  pairs. 

4.24  10  Antenna  Pattern  FTTects  Removal 

The  primary  output  from  the  antenna  effects  removal  is  the  array.  This  array  is 

passed  into  additional  processing  algorithms,  as  shown  in  figure  .T1 , and  is  also  output  in  numerical 
and  plotted  formats.  These  formats  are  identical  to  those  previously  described  for  the  uncorrected 
Sir. cut  function.  A secondary  output,  Gtr.cu),  which  provides  a visual  indication  of  the  antenna's 
spatial  filtering  properties  as  mapped  into  the  (r.cu)  coordinate  system,  is  given  in  figure  4-27  and 
illustrates  that  for  this  representative  case  very  little  filtering  oceurs  over  the  realm  a'sponsible  for  the 
rnaionty  of  the  multipath  return  (i.e..  the  shoulders  of  the  S(r.cu)  function). 


4-8.) 


TABLE  4-17.  NOISE  PARAMETER  OUTPUTS  - M(fJ.  R(f).  N(f).  A„, 
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4,24.1  I Tap  I and  Q Amplitude  Distribution 


The  amplitude  probability  distribution  algorithm  is  applied  to  both  1 and  0 and  for  each  com- 
ponent has  as  its  primary  output  the  following  information: 

a.  Numeric  tabulation  of  each  of  the  three  histogram’s  observed  frequencies  and  expected 

frequencies 

•> 

b.  The  X"  value  and  probability  for  each  of  the  histograms 

c.  A plot  of  the  observed  and  expected  distribution  associated  with  the  histogram  con- 
structed by  subdividing  the  t-3o,  +3o)  amplitude  range  into  18  equal-length  bin  widths 

An  example  of  outputs  a and  b is  presented  in  table  4-19.  This  represents  that  portion  of  the 
output  associated  with  the  histogram  corresponding  to  bin  width  selection  based  on  equal  bin  fre- 
quency expectation  (e.g.,  5%  expected  in  each  bin;  total  population  1024). 

Figure  4-28  contains  an  example  of  output  option  c.  Plotted  on  one  page  are  the  1,  Q,  and 
phase  distributions  for  two  taps.  Bar  charts  and  continuous  curves  are  used  to  represent  the  observed 
and  expected  frequencies,  respectively.  For  the  1 and  Q distribution,  the  markers  along  the  x-axis  are 
0,  3 apart  and  span  the  range  from  -3o  to  +3o.  The  scale  mark  for  the  phase  distribution  is  given 
every  20“  over  the  range  from  -180“  to  +180°. 


4.24. 1 2 Tap  Phase  Distribution 


Numerical  and  plotted  output  formats  are  available  to  present  the  results  of  the  tap  phase  dis- 
tribution testings.  An  example  illustrating  the  plotted  format  has  been  given  in  figure  4-28.  The 
numerical  output  consists  of  a tabulation  of  the  observed  and  expected  histogram  values,  the  calcu- 
lated  x“  statistic,  and  the  probability  of  obtaining  a x"  value  less  than  the  observed.  Table  4-20 
contains  a representative  numerical  output  set. 


4.24.13  Tap  Process  Cross-Correlation 


For  complex  time-domain  processes,  such  as  represented  by  the  SACP  tap  voltages,  the  cross- 
correlation function  has  both  real  and  imaginary  components.  The  tap  process  cross-correlation 
algorithm  numerically  outputs  these  results  in  the  form  of  magnitude  and  phase  and  also  provides  a 
three-dimensional  plot  of  the  magnitude.  An  example  of  a numerical  tabulation  is  given  in  table  4-21 . 
Columns  are  lumped  in  groups  of  three  ( L AG.  AMP,  PH)  corresponding  to  time  lag  in  microseconds, 
function  normalized  amplitude,  and  phase  in  degrees.  A sample  plotted  output  is  given  in  figure  4-29. 
This  output  format  is  configured  to  display  on  one  graph  all  cross-correlations  that  are  made  with  a 
particular  tap  (referred  to  as  the  principal  tap) 
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TABLE  4-19.  TAP  ! AND  Q APD  HISTOGRAMS 
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TABLE  4-21.  CROSS-CORRELATION  NUMERICAL  OUTPUT 
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4.24  14  Tap-Guin  AutOL'orrelation  Function 


Numerical  and  three-dimensional  plot  formats  are  used  to  represent  the  derived  tap-gain  auto- 
correlation function.  Both  magnitude  and  phase  are  given  numerically,  whereas  only  the  magnitude  of 
the  function  is  plotted.  Fxamplo  numerical  results  for  the  function  magnitude  and  phase  are  given  in 
tables  4-22  and  4-2,1.  respectively.  The  output  format  in  both  cases  is  self-explanatory  (the  time-lag 
column  is  in  milliseconds).  This  particular  example  corresponds  to  the  noise-removed  option  applied 
to  horizontally  polarized  in-plane  flight  data  and  is  seen  to  be  characterized  by  a process  whose 
specular  return  falls  in  tap  25.  Corresponding  to  the  same  data  base  used  for  the  numerical  example, 
we  present  a sample  plotted  output  in  figures  4-,10  and  4-31  for  noise-present  and  noise-removed 
options,  respectively.  As  indicated,  the  receiver  noise,  which  appears  as  a spike  train  along  the  f = 0 
axis,  is  effectively  eliminated  by  the  noise-removal  algorithm. 

4.24.15  Airborne  System  Parameters  Tape 

Numerical  line  print  and  stacked  CalComp  plots  are  used  to  display  the  airborne  system 
parameters.  | 

The  printed  output  example  presented  in  table  4-24  illustrates  the  presentation  format.  Each 
page  begins  with  the  experiment  date.  Underneath  is  a columnar  listing  of  the  time  (HR/MIN/SHC); 
the  status  bits  for  the  bit  rate,  sequence  length,  antenna,  and  polarization:  followed  by  the  13  data 
channel  samples  corresponding  to  the  power  monitor  sensors  and  the  aircraft  flight  parameters. 

CalComp  plots  stacked  vertically  over  a common  time  base  are  used  to  display  the  SACP 
modulator  mode  data,  the  RF  subsystem  data,  and  any  nine  of  the  aircraft  (light  parameters  and 
power  meter  readings.  These  data  are  plotted  with  a time  scale  of  1 in/min.  In  general,  one  plot  is 
used  to  represent  a single  flight  test  leg  ( 14  to  20  min).  Figure  4-32  presents  a typical  output  plot. 
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Figure  4-30.  Tap-Gain  Autocorrelation  Function  {Noise  Present) 


TABLE  4-24.  AIRBORNE  TAPE  NUMERICAL  OUTPUT 
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Figure  4-32.  Airborne  Tape  Plotted  Results 
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5.  PHYSICAL  OPTICS  SCATTER  MODEL  PREDICTIONS 


5.1  MODEL  INTRODUCTION 

The  choice  of  an  appropriate  model  for  the  analysis  of  electromagnetic  surface  scatter  is  deter- 
mined almost  exclusively  by  the  rougliness  characteristics  of  the  retlecting  medium.  Surfaces  are 
usually  classified  as  slightly  rough,  very  rough,  or  composite  depending  on  the  magnitude  ot  the  height 
irregularities.  In  general,  different  scatter  theories  are  used  for  each  of  these  situations. 

For  the  case  of  L-band  aeronautical  satellite  oceanic  scatter,  the  surface  will  almost  always 
appear  to  he  very  rough;  this  implies  that  the  following  is  approximately  satisfied: 

> aucosd.  > 1.0  , (5-1) 

where : 

X = electromagnetic  wavelength 

0|l  = standard  deviation  of  surface  height  irregularities 

= incident  angle  of  ray  upon  the  surface  as  measured  from  the  normal . 

Analysis  of  scattering  from  very  rough  surfaces  is  usually  developed  through  the  physical 
optics  tangent-plane  method.  Commonly  called  the  Kirclihoff  approximation,  this  model  is  based  on 
the  assumption  of  a locally  plane  surface  over  the  distance  of  many  wavelengths.  This  constraint  is 
considered  to  be  satisfied  if  the  radius  of  curvatuve  of  the  surface  undulations  (5^,)  is  much  greater 
than  X ; i.e., 

X « 5^,  cos  Oj  . (5-2) 

For  this  analysis  we  employ  the  vector  formulation  of  the  physical  optics  model  and  are  thus 
able  to  properly  account  for  the  electromagnetic  polarization  dependences  of  each  particular  scatter- 
ing facet  on  the  surface.  Because  of  the  immense  complexity  of  this  model,  it  is  not  possilble  to  arrive 
at  adequate  channel  parameter  solutions  in  a closed  form.  This  is  circumvented  through  use  of  a com- 
puterized techni(iue  that  subdivides  the  sperical  scatter  surface  into  incrementally  small  areas  and  then 
determines  the  scatter  cross  section  (including  polarization  transfonnation  factors).  Doppler  shift,  and 
time  delay  associated  with  each  area.  The  complex  vector  representation  of  the  scattered  signal  is 
coupled  to  the  receiver  antenna  charaeteristics.  thereby  providing  an  estimation  of  the  received  power 
from  the  particular  surface  path  and  thus  allowing  the  channel's  delay-Doppler  scatter  function. 
S(T.tc),  to  be  constructed.  From  S(t.oj).  steps  identical  to  those  described  in  section  4. 1 2 are  fol- 
lowed to  determine  the  channel's  time-frequency  autocorrelation  function,  delay  spectrum.  Doppler 
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spectrum,  time  autocorrelation  function,  frequency  autocorrelation  function,  total  energy  content, 
and  spread  values  of  the  unidimensional  distrihutions.  Described  below  are  the  most  important  scatter 
theory  relationships  of  the  computer  routine  as  well  as  typical  sample  outputs  in  a plotted  format. 


5.:  ALGORITHM  LXLCUTION  SEQUKNCE  AND  DESCRIPTION 

A block  diagram  illustrating  the  computer  model’s  input/output  and  basic  processing  function 
is  given  in  figure  5-1 . Inputs  to  the  program  include: 

a.  System  geometry  parameters  (aircraft  altitude,  speed,  velocity  vector  direction,  grazing 
angle  at  specular  point,  and  transmitted  frequency) 

b.  Complex  polarization  vectors  of  airplane  and  satellite  antenna  (satellite  polarization 
characteristics  are  assumed  to  be  isotropic  over  the  scatter  surface;  antenna  radiation 
distribution  characteristics  may  be  accounted  for  mathematically  or  by  including  a 
azimuthal-elevation  angle  look-up  table  to  derive  complex  horizontal  and  vertical 
polarization  coefficients) 

c.  RMS  surface  slope  and  electrical  parameters 

d.  Surface  area  resolution  parameter  (detennines  the  size  of  the  incremental  surface  area 
used  in  the  surface  integration  routine  and  thus  ultimately  affects  the  fidelity  of  the 
predictions) 


Ihe  basic  processing  lunction  of  the  computer  code  may  be  broken  down  into  three  categories. 

a.  Lor  each  incremental  surface  area  (dS)  considered  in  the  integration,  calculate  the 
received  scattered  power  (dPi.  the  time  delay  (r),  and  Doppler  shift  (cu)  associated  with 
the  return. 

b.  Sort  the  power  returned  from  each  incremental  area  into  a twxvdimensional  array  (t.co) 
to  yield  an  estimate  of  the  channel's  scatter  function  S(r.cu). 

c.  Perform  Eoiirier  and  inlegnil  operations  on  S(r.cc)  to  derive  Ihe  channel's  tune-frequency 
autocorrelation  lunction.  time  aiilocorrel.ilion  function,  frequency  autocorrelation  func- 
tion. Doppler  spectrum,  delay  spectrum.  rms  scattered  energy,  and  the  spre.id  values 
of  the  uniilimensional  parameters. 


I 


For  incremental  surface 
area  determine. 

• received  power  (dP) 

• time  delay  |r) 

• Doppler  shift  (cj) 


system  geometry 

complex  polarization  vectors  of  airplane 
and  satellite  antenna 

rms  surface  slope  and  electrical  characteristics 
surface  area  resolution  parameter 


Sort  dP  into 
S (7, cu)  array 


Repeat  for  total 
significant  scatter 
region 


i: 


Integral  and  Fourier 
operations  on  S (7.  to) 


Spread 

calculations 


Output 

• delay-Doppler  function 

• time-frequency  correlation  function  R ($,  il) 

• Doppler  spectra  D (to) 

• delay  spectra  Q i7) 

• time  autocorrelation  function  R (f  , 0) 

• frequency  autocorrelation  function  R (0,  12) 

• spreads  of  D (to),  Q (7),  R (^ , 12),  R (0,  S2) 

• total  rms  scattered  energy 


Figure  5-1.  Scatter  Mode! 


Category  c utilizes  algorithmic  operations  identical  to  those  described  in  sections  4. 1 2 and 
4. 1 3.  The  second  category  ( b)  is  a relatively  straiglitforv^ard  two-dimensional  bin-sorting  algorithm 
and  needs  no  further  explanation. 

The  first  category  (a)  represents  the  surface  scatter  boundary-interlace  component  ol  the 
model  and  is,  to  a large  extent,  based  on  the  work  of  Peake  ( ref.  5- 1 ).  Draw-ing  Iron)  this  source  we 
brietly  present,  with  the  aid  of  figure  .‘^-2.  the  underlying  mathematics  associated  with  the  derivation 

ofdP,  T.  and  w.  For  ease  of  interpretation,  parameter  dP  is  calculated  relative  to  the  power  received 

over  the  direct  path  link  whereas  t and  w are  evaluated  relative  to  the  direct  path  Doppler  and  the 
specular  point  delay  values,  respectively.  The  normalized  received  power  trom  incremental  area  dS  is 
thus  expressed  as 


1 ’’d 

dP  = — 

47T 


(\fTcV 


ri“rs“  IGr 


■Gsl- 


'^xs  ‘IS, 


(5-3) 


where : 


ri,r2.rd 

^cf 


i’P 

^s 

‘'xs 

dS 


= distances  as  illustrated  in  figure  5-2 

= complex  coefficient  that  accounts  for  coupling  between  incident  polarization  vector. 

tilted  dielectric  surface  element,  and  receiver  polarization  vector 
= airplane  complex  polarization  vector  for  direct  link  path 
= satellite  complex  polarization  vector 
= surface  scatter  cross  section 
= elemental  surface  area. 


The  complex  transmission  coefficient  embodies  the  aircraft  antenna  spatial  tiltering  and  is 
calculated  as; 


T(.,-  = (DI ) I R||  (n)  ( A cos  Sj  + B sin  5j ) (C  cos  + D*in  6^) 
- R|  (a)  (B  cos  5j  - A sin  5j)  (-D  cos  + C sin  . 


154) 


5-4 


where; 


Dl 


= spherical  earth  divergence  factor 

R (a)  = parallel  Fresnel  reflection  coefficient 

II 

= (e' cos  a - E)  (e"  cos  a + F.) 

R|(a)  = perpendicular  Fresnel  retlection  coefficient 

= ( cos  a - E ) ( cos  a + E ) 

€ = complex  index  of  refraction  of  the  surface 

a = angle  between  incident  wave  and  normal  of  properly  tilted  (to  produce  retlection 

into  receiver)  surface  facet 

E = (e~  - sin-  a)  ' 

6j  = angle  between  theta  component  of  incident  wave  and  the  incident  parallel  unit 

vector  (i.e..  vertical  polarization  with  respect  to  tilted  surface) 

= angle  between  the  theta  component  of  scattered  wave  and  the  reflected  parallel  unit 
vector 

A.B.C.D  = transmitter  and  receiver  complex  antenna  polarization  vector  coefficients  as 
described  below. 


Polarization  vectors  for  the  transmit  iP-j  ) and  receive  (P|^)  antenna  system  are  given  by : 


P'l'  — A L'l'-j'  + B Lpy 


_ A A 

= C L-p^  + D L 


PR’ 


(5-5) 


where: 

A A 


Lit- Lpj 


unit  vectors  in  the  theta  and  phi  directions,  respectively,  with  respect  to  coordinates 
centered  on  transmitter 


A A 

LyR.  Lpp  = unit  vectors  in  the  theta  and  phi  directions,  respectively,  with  respect  to  coordi- 
nates centered  on  receiver 


The  physical  optics  very  rough  surface  scatter  cross  section  is  obtained  from: 


(5-(i) 


where: 


(9  .6 
n ^n 


'V»n'<-n* 


angles  as  illustrated  in  figure  5-2 
surface  slope  probability  density  function  . 
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i 


j 

For  this  particular  analysis  the  slope  probability  density  function  was  assumed  to  be  isotropic  with  the  I 

following  Gaussian  form: 


cos 

Pss<®n-'^’n>  - 

vq~ 


(- 


r?“(l+27j-)/ 


(5-7) 


where  rj*  = <tan-  = mean  square  surface  slope. 

The  delay  and  Doppler  variables  associated  with  a given  incremental  area  are  determined  from 


:7Tf, 


-np  - v + nd  -V 


n +':~'ls’''2s 


(5-8) 

(5-9) 


where  f^  is  the  transmitted  frequency,  c is  the  velocity  of  light,  and  the  other  parameters  are  as 
illustrated  in  figure  5-2. 


5.3  SAMPLE  OUTPUT 

To  facilitate  the  comparison  between  theoretical  expectation  and  experimental  results,  the 
scatter  channel  model’s  output  formats  have  been  configured  to  resemble  as  closely  as  possible  those 
associated  with  the  multipath  data  reduction  routines.  Thus,  for  the  most  part,  the  applicable  output 
format  description  given  in  section  4.23  is  pertinent  to  this  section.  To  illustrate  the  degree  of  model 
sophistication,  a set  of  sample  plotted  results  associated  with  a typical  production  nin  is  presented 
here. 

Input  parameters  for  this  sample  run  were  as  follows: 


Satellite  polarization: 
Airplane  polarization: 
Aircraft  velocity: 
Aircraft  altitude: 
Frequency 
Surface  type: 

Surface  rms  slope : 
Grazing  angle: 


Right-hand  circular 
Horiz.ontal-omni 

200  m 'sec.  great  circle  path  toward  satellite 

10  km 

lb50  Mllz 

Sea  water 

()0 

iqO 
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A cross-reference  relating  figure  number  with  the  model’s  parameter  prediction  is  given  below: 


Figure 
5-3 
54 
5-5 
5-6 

It  is  of  interest  to  note  that  the  example  run  entailed  a division  of  the  effective  surface  scatter 
area  into  more  than  roughly  80,000  elemental  patches  and  required  a total  of  380  sec  of  CDC  6600 
central  processor  time. 

SCATTCR  ruNCMON 
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Figure  5-3.  Scatter  Model  Prediction  - S(t,cj) 


fiimnun  maiimum  i«c«t«tnT 
o 100  to 

-Ito  too  *0 

e 1000  loo 


Parameter 

S(T,co) 

D(oj),  Q(t) 

iR(i.n)i 

|R(0,n)l,  |R(J.0)|. 


1 


5-8 


a 

m 


0 1.000  a.OOO  9.000  ^.OOO  9.000  i.OOO  /.OOO  0.000  t.ooo 


110'* 

fine  SCfAO»T|OH  tSCCONDSl 

Figure  5-6.  Scatter  Model  Predictions  - I R(0,m  I,  I I 


5-1  1/5- 1 : 


L 


6 MODEM/ANTENNA  DATA  PROCESSING 


The  data  gathered  for  the  digital  modem,  voice  modem,  hybrid  voice  and  data  modem,  ranging 
modem,  and  antenna  tests  on  the  KC-135  airplane  is  in  the  form  of  analog  tapes,  written  logs,  and 
onboard  strip  charts.  Processing  of  this  raw  data  is  described  below. 


6.1  SYSTEM  BLOCK  DIAGRAM 

Figure  6-1  illustrates  the  How  of  data  in  the  data  reduction  and  analysis  system.  The  diagram 
shows  the  various  machine  processes  and  human  judgment  processes  constituting  the  total  system. 

The  airborne  logs  and  strip  charts  are  used  to  determine  the  valid  voice  test  periods.  The 
selected  voice  channels  are  transcribed  to  I /4-in.  tapes  for  intelligibility  scoring  by  CBS  Laboratories. 

Tape  logs  are  used  to  determine  input  parameters  that  are  furnished  to  the  Boeing  Test  Data 
Processing  Center  iTDPCi  along  with  the  analog  source  tapes.  These  parameters  include  type  of  test 
data,  tracks  to  process,  time  segments  to  process,  analog  source  tape  speed,  and  kind  of  output 
desired  (digital  tape  or  stripout). 

Airplane  logs,  onboard  strip  charts,  and  TDPC  strip  charts  are  used  to  determine  which  time 
segments  ot  the  digital  tapes  are  to  be  analyzed.  The  time  segment  information  is  sent  with  the  digital 
tape  to  the  CDC  6600  for  processing  with  the  applicable  analysis  programs. 

Output  from  the  CDC  6600  includes  C/Nq  and  S/I  statistics,  bit-error  rates,  bit-error  statistics, 
and  ranging  data  analysis.  This  output  is  manually  merged  with  other  data  such  as  voice  intelligibility 
scores  and  aircraft  parameters  to  give  final  outputs.  Outputs  include  those  described  in  the  DRandA 
plan. 


6.2  TEST  DATA  PROCESSING  CENTER  (TDPC) 

Processing  of  the  analog  data  at  TDPC  is  subdivided  into  two  basic  systems,  with  each  system 
having  modes  dependent  on  the  data  being  processed.  One  system  processes  antenna,  voice,  and 
ranging  data;  the  other  system  processes  the  digital  communication  data.  1 he  systems  are  described 
more  fully  in  the  following  paragraphs. 
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Analog  source  tape  KC  135  Jh  ATS6 


Figure  6- 1.  Modem  Evaluation  DR  and  A System  Block  Diagram 


0.2.1  rOPT  Processing  of  .Antenna.  V'oicc*.  and  Rjnging  1 apes 


A bloL-k  diagram  ol  the  system  used  for  all  three  types  of  data  is  shown  in  figure  0-2.  A brief 
explanation  of  the  formatting  operation  for  the  antenna,  voice,  and  ranging  data  is  given  below. 

Antenna  data  on  the  source  tape  consists  ol 

IRIC-B  time  code  WBFM 

Signal  strength  (carrier  detector)  WBFM 

Oiannel  event  (A  or  B)  WBFM  . 

Tliese  signals  are  recorded  at  3-2/4  ips.  At  the  transition  of  the  event  signal,  time  is  read  and  1025 
samples  of  signal  strength  are  acquired  and  formatted  on  the  tape  as  shown  in  figure  0-3. 

Voice  data  is  recorded  at  3-3/4  ips  on  the  source  tapes  and  consists  of  the  following  signals; 

IRKFB  time  code  WBFM 

Signal  strength  (carrier  detector)  WBFM  . 

For  the  voice  data,  sampling  of  the  signal  strength  is  under  direct  control  of  the  time  code  rather  than 
by  a channel  event  signal.  Beginning  I 5 sec  after  processing  start  time  and  repeating  at  1 5-sec  intervals, 
time  is  read  and  1025  signal  strength  samples  are  taken  and  formatted  on  digital  tape,  again  as  shown 
in  figure  (i-3. 

Ranging  data  is  recorded  at  3-3/4  ips  on  the  source  tapes  and  consists  of  the  following  signals: 

IRIG-B  time  code  WBFM 

Signal  strength  WBFM 

Ranging  data  PC'M  WBFM,  Manchester  code.  305  bps  . 

Upon  detection  of  a seven-bit  synchronization  pattern  (Barker  pattern)  in  the  PCM  track,  the  next 
33  bits  of  PCM  data  are  fonnatted  with  time  and  an  identification  (ID)  word  on  digital  tape  as  shown 
in  figure  t>-4.  Fach  of  these  blocks  consists  of  ten  1 2-bit  words.  After  1 b of  these  blocks  have  been 
formatted,  1025  signal  strength  samples  are  taken  and  formatted  on  the  digital  tape  (I'lg.  b-3).  fhis 
completes  one  ranging  data  record.  .Ml  these  refonnatted  digital  tapes  are  seven-track,  8()()-bpi  tapes 
compatible  with  the  CI)C  (itiOO. 

The  input  to  the  system  diagrammed  in  figure  b-2  is  the  source  tape  recorded  on  the  airplane. 
The  tapes  are  played  back  at  30  ips.  which  is  a speed  increase  of  X8.  fhe  signals  applicable  to  the 
data  being  processed  are  demodulated  by  the  DC.S-2  demoiliilator.  I he  time  is  patched  to  the  range 
time  decoder  (R  I'D)  and  cross-coupled  to  the  digital  data  converter  (DIX  ).  The  signal  strength  and 
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111111111111 
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voice 

digital  data  bit  stream 
digital  data  amplitude  samoles 
ranging 

antenna  evaluation  (antenna  1) 


A = Event  track  status  bit  (antenna  data  only)  in  Words  6 to  1030 
1 = event  track  on 
0 = event  track  off 

A = 0 for  voice,  ranging  and  digital  data 


Figure  6-3.  Signal  Strength  Format,  800-bpi  Digital  Tape 


event  tracks  are  patched  to  the  DDC  through  the  analog  patch  bay  (APB).  The  serial  PCM  (ranging 
data)  is  decomnuitated  by  the  TMR  bit  synchronizer  and  clocked  into  the  DDC. 


The  signal  strength  data  are  sampled  and  digitized  by  the  DDC.  A 16-kHz  clock  (2-kHz  real 
time)  triggers  an  analog-to-digital  converter,  and  the  digitized  samples  are  transferred  to  the  SCC  4700 
computer  througli  the  direct  memory  access  (DMA). 

The  7-bit  PCM  synchronization  patterns  preceding  the  33-bit  ranging  data  bursts  are  detected 
by  the  DDC.  Following  detection,  the  next  33  bits  of  PCM  data  are  transferred  to  the  SCC  4700  com- 
puter in  one  9-bit  word  followed  by  two  1 2-bit  words. 

The  ATSVAR  computer  program  formats  time,  signal  strength,  and  PCM  data  into  records. 
ATSAVR  is  the  SCC  4700  computer  program  that  is  operator  controlled  through  teletype  input.  The 
operator  selects  the  mode  of  operation  according  to  the  type  of  data  to  be  processed  and  selects  the 
number  of  output  tapes  necessary  according  to  the  quantity  of  input  data. 

j Data  is  double  buffered  so  that,  while  one  record  is  being  completed  and  written  on  tape,  the 

lormatting  ot  the  next  may  begin.  This  is  necessary  in  the  ranging  mode  where  additional  synchroni- 
; zation  patterns  may  occur  before  the  signal  strength  sampling  of  the  previous  record  is  completed. 

Time  editing  is  accomplished  by  the  start-stop  switches  on  the  range  time  decoder.  Processing 
may  also  be  started  or  stopped  manually.  Each  processing  interval  constitutes  one  file  of  data  and  at 
the  conclusion  of  all  processing  intervals  the  operator  terminates  the  processing  mode.  Data  files  arc- 
separated  by  one  end-ot-tile  mark.  I he  last  data  file  on  the  tape  is  terminated  by  two  or  more  end-of- 
tile  marks. 


().2.2  Digital  Communication  Data 


The  digital  data  is  recorded  on  the  source  tapes  at  either  3-3/4  or  7-1/2  ips  and  consists  of  the 
following  signals: 


Five  tracks,  biphase  level  Pf'M 
Signal  strength  (carrier  de(ector) 
1 RKi-B  time  code 


Direct  record,  one  track  per  modem 

WBFM 

WBFM 


i 


I 
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I hc  data  rate  of  the  distital  data  is  I 200  hps.  \\hieh  results  in  tlie  following  relationship  of 
record  and  playhaek  speeds; 


Analog  souree  tape 
record  speed,  ips 


1 DIH  playhaek 
speed. ips 


During  processing  of  the  digital  data  analog  tape  tl'ig.  h-5l.  the  P(  M data  is  routed  to  the  DDC  . 
the  analog  signal  strength  to  the  DC.S-2  demodulator,  and  time  to  the  K 1 D.  The  signal  strength  is 
passed  througli  a demodulator  and  sent  to  the  DDC.  where  It  is  converted  to  a digital  value  by  an  ana- 
log-to-digital  converter  and  written  on  tape,  one  record  approximately  every  15  data  seconds. 

Five-hit  synchroni/.ers  have  been  patched  on  the  DDC  logic'  board,  one  lor  each  recorded 
modem  track.  Ihese  bit  synchronisers  obtain  bit  synchronisation  and  resolve  Manchester  phase 
ambiguity  I he  regenerated  NR/  PCM  is  written  on  tape  in  six-hit  bytes,  compatible  with  seven-track 
digital  tape 

The  frequency  synthesiser  on  the  analog  patch  bay  is  set  to  I kUs  times  the  playback  rate 
(X4)  and  patched  to  DDC  input  line  This  signal  is  routed  to  a threshold  detector.  The  leading  and 
trailing  edges  are  OR'ed  together  to  give  a 2lx  pulse/sec  sample  clock  I liis  sample  clock  is  used  for 
the  A D conversion  ol  the  signal  strength  track. 

I he  tape  formal  of  the  digital  data  is  illustrated  in  figure  b-(i. 


b.2..?  I DPC  Chart  Stripouts 

The  analog  source  tapes  contain  some  d.ita  desirable  to  examine  in  analog  visual  form 
Accordingly,  the  applicable  channels  are  stripped  out  to  yield  the  d.ita  sliown  in  table  (i-l 

■All  the  above  signals  are  recorded  on  chart  paper  for  voice,  data,  and  hy  brid  voice  and  data 
tests.  Only  the  aircralt  parameters  are  stripped  out  for  ranging  and  antenna  tests. 


()..?  ( l)(  ()(>()()  PKO(  1 SSINt, 

I he  data  tapes  lorm.itted  by  I DPt  are  processed  on  a CDC  (>ol)l)  computer  under  Ihe 
KRONOS  2.1  operating  system  f igure  b-7  gives  the  d.ita  How  of  a compiitei  run  and  is  valid  for 
all  data  processing  iligit.il.  r.mging.  voice,  .ind  .mtenna 


Figure  6-5.  TDPC  Digital  Data  Tape  Formatting  Process 


DIGITAL  DATA  FORMAT 


DATA  RECORD 


SIGNAL  STRENGTH  RECORD 


12  BITS 


12  BITS 


5 TIME  WORDS 


CHAN  1 
205  WORDS 


CHAN  2 
205  WORDS 


CHAN  3 
205  WORDS 


CHAN  4 
205  WORDS 


CHAN  5 
205  WORDS 


1030 

WORDS 


5 TIME  WORDS 


1025 

SIGNAL 

STRENGTH 

SAMPLES 


PCM  DATA  WORDS 


1 


10 


11 


12 


DATA  BIT  1 IS  THE  FIRST  BIT  TRANSMITTED 

SIGNAL  STRENGTH  IS  SAMPLED  EVERY  14  SECONDS  (APPROXIMATELY) 

FOR  1200  BPS  THERE  ARE  7 DATA  RECORDS  PER  SIGNAL  STRENGTH  RECORD 


Figure  6-6.  Digital  Data  Test  Tape  Format 


(>- 1 0 
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TABLE  6-1.  AIRCRAFT  PARAMETERS 


Information 

Format 

Aircraft  heading  (coarse) 

7.5%  proportional  to  subcarrier  FM 

Aircraft  heading  (fine) 

7.5%  proportional  to  subcarrier  FM 

Aircraft  pitch 

7.5%  proportional  to  subcarrier  FM 

Aircraft  roll 

7.5%  proportional  to  subcarrier  FM 

Hybrid  No.  1 C/N^^  (dc  voltage) 

WBFM 

Hybrid  No.  2 C/N^  (dc  voltage) 

WBFM 

The  peripheral  equipment  necessary  to  run  the  program  are  card  reader,  printer,  and  two  tape- 
drives  (one  for  the  reformatted  input  data  tape  and  one  for  the  output  plot  file).  Plotting  is  an  off-line 
process  performed  on  an  SC  4020  Stromberg-Carlson  plotter. 

Although  the  programs  are  primarily  coded  in  FORTRAN,  a few  subroutines  are  coded  in 
assembly  language  to  minimize  run  time. 

Programs  are  run  in  a batch  processing  mode.  The  program  inputs  consist  of  the  TDPC  digital 
magnetic  tapes  and  control  data  cards.  These  cards  include  start  and  stop  times  of  each  segment  and 
program  options.  The  outputs  are  a line  printer  listing  and  a plot  of  the  signal  spectrum.  A more 
detailed  description  of  program  input/output  is  given  under  each  program  section. 


Figure  6-7.  CDC  6600  Computer  Data  Flow  for  Modem  Evaluation  Data  Processing 


7.  ANALYSIS  OF  ENVELOPE  DETECTOR  OUTPUT  TO  DETERMINE  C/N^  AND  S/I 


This  section  describes  in  detail  the  inetliodology  and  underlying  theory  tor  the  determination 
of  the  carrier-to-noise  power  density  ratio.  C/N^.  and  the  direct  signal-to-niultipath  signal  ratio.  S/I. 
The  waveform  utilized  is  the  low-pass  output  of  the  carrier  detector  unit.  The  carrier  detector  unit 
processes  an  unmodulated  carrier  multiplexed  into  the  satellite-aircraft  transmission.  The  procedure 
is  involved  in  analysis  of  voice,  data,  ranging,  and  antenna  experiments. 


7. 1 BACKGROUND  THEORY 

The  carrier  detector  unit  incorporates  a bandpass  linear  envelope  detector  (full-wave  rectifier) 
to  provide  a measurement  of  channel  character.  The  nonlinear  approach  was  adopted  to  obviate  the 
necessity  for  continuous  Doppler  compensation  if  a linear  down-translation  and  analysis  were  used. 
Also,  data  sampling  rates  can  be  somewhat  smaller  with  this  approach. 

The  detector  input  signal  is  assumed  to  be  the  sum  of  ( I ) a sinusoid  of  frec|uency  f, , with 
average  power  P-/2.  i.e..  its  amplitude  is  P volts  with  a l-ohm  reference.  ( 2)  a diffuse  multipath 
component  having  center  freiiuency  f^,  and  power  Oj,,".  and  (.^ ) wideband  additive  Gaussian  noise 
with  single-sided  density  of  W Hz. 

Commensurate  with  the  commonly  accepted  scatter  model  applicable  for  an  aeronautical 
satellite  oceanic  system,  the  multipath  return  is  assumed  to  have  a complex  Gaussian  representation, 
with  spectrum  * given  by 


sr-(„)2/2B„2 


The  parameter  is  the  one-sided  e'  bandwidth  of  the  scattered  signal.  I heoretical  and  empirical 
results  indicate  that  B^.|  < 1 2.“'  Hz  lor  the  elevation  angle/airspeed  profile  llown  in  the  experiment. 

The  composite  signal  is  initially  prefiltered  with  a four-pole  Bessel  crystal  filter  to  reject  wide- 
band noise  and  pass  the  signal  components.  I he  lilter  is  modeled  as  having  a Gaussian  shape  with 

* A symmetric  spectrum  is  obtained  theoretically  for  cross-plane  llights  as  in  l ype  II  tests  but  is  not 
strictly  true  for  other  orientations. 
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bandwidth  parameter  for  analytical  simplicity.  (Bessel  and  Gaussian  characteristics  are  very 
simitar  to  the  30-dB  attenuation  points.)  Since  noise  and  multipath  are  independent,  the  total  inter- 
ference spectrum  at  the  detector  input  then  becomes 


^ «n^  -(f-fo)2/2B„2 


Measurements  have  determined  that  Bj^  = 2130  Hz  for  this  experiment,  with  the  result  that 
all  multipath  energy  is  passed  by  this  prefilter.  Figure  7-1  shows  the  total  spectrum,  including  the  cw 
signal  term. 

The  detector  is  modeled  as  an  ideal  linear,  meinoryless.  full-wave  device;  i.e.. 


y(t)  = Ar(t), 


y(t)  = output 
rtt)  = input 

A = arbitrary  gain  constant . 

2 

OfP  = multipat)>  power 

2 

Op  = noise  power 

= e '^  multipath  bandwidth,  one  sided 

;(f)  B = e '^  prefilter  bandwidth,  one  sided 


P^/2  watts  . 


\ X/57V  B„  B„  ) 


Figure  7- 1.  Detector  input  Spectrum 
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Laboratory  measurements  have  indicated  that  the  actual  detector  is  quite  linear  over  a +8  dB  range 
about  the  nominal  signal  level. 

Finally,  the  output  is  low-pass  filtered  with  a two-pole  Butterworth  filter  whose  equivalent 
one-side  rectangular  bandwidth  is  800  Hz.  This  signal  is  recorded  for  postflight  analysis.  The  objective 
of  the  analysis  is  to  determine  C/N^  and  S/I  of  the  signal. 

Following  Rice  (ref.  7-1 ),  Middleton  (ref.  7-2),  and  Davenport  and  Root  (ref  7-3).  the 
unfiltered  output  autocorrelation  is  expanded  in  a power  series  involving  the  input  autocorrelation; 


Ry(f ) + 2 *^mo^  ~jJ7 

m = 1 k = 1 

.2^ 


I cos  mcOgT. 


(7-4) 


m = I k = 1 


This  expansion  displays  the  mean-value,  signals  harmonics,  noise-cross-noise,  and  signal-cross-noise 
contributions,  respectively. 

The  coefficients  h^|^  are  a function  only  of  the  total  signal-to-input  interference  ratio. 
P"/2o",  not  the  spectral  distribution.  The  latter  may  be  determined  from 


(p2/2a^)'‘  = (p2/2a„2y  ^ (p2/2o„2y 


(7-.S) 


Now.  for  the  assumed  Gaussian  prefilter. 


(s/N„, 


So 


(7-M 


(p2/2o2)-'  = (S/I)-'  . |;s/N„) 


(7-7 1 
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llie  coet'fiL'ients  are  ( ret.  7-.1 ): 


^mk  “ ; T 1 k-1 

(m!)“  r-  ((3  - m - kt/2)  (o't 


(7-8) 


where  | F|  is  the  eontluent  hypergeometrie  function. 


By  transforming  equation  (7-4),  the  output  spectrum  is  obtained  as 

oe 

Sy(w)  = hoo^  ^ [6(w  - mw^,)  + 6(w  + m^o)] 

m = I 

OO  2 

+ ^ ^ -2^  [s^lw  - mWo>  + + mwo)]  , 


(7-4) 


k = 1 m = 0 


where  S|^(u>)  is  the  k-fold  convolution  of  the  input  interference  spectrum  given  by  equation  (7-2). 
In  this  convolution  of  spectra  lies  the  appeal  of  the  Gaussian  input  spectra. 

Although  the  spectrum  expansion  appears  rather  forniidable,  the  series  converges  rapidly  for 
the  cases  of  interest  so  that  only  a few  terms  need  to  be  considered  (m  < 2.  k < 3).  Also,  only 
products  with  (m  + kteven  yield  energy  in  the  low-frequency  zone. 


Using  the  spectrum  of  equation  (7-2),  the  following  relations  can  be  written  for  the  convolved 


spectra: 


2J2tr 


S-,(f)  = 


^ r^-(f-fo)^/2B„\^4f>fo)‘/2B„^] 

2vT:rr  B„ 

. 2fo)2/4B^‘  ^ ^ ^-(f  + 2f„)-/4B,„‘] 

4>/4»r  B_. 


_^r  4f-2fo)2/4B„2^,^-f2/4B,^^^.-(f+2tV‘/4B„-] 

4/4. 
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, --  ^ [,^f  ■ ♦ B„2|  ^ ,^-f^/2(B„2  , b„2, 

-(r*2r„|2/aB„2»B„2|] 


4j4w 


+ e 


S3,f, . [,-lf-2f„)2/6B^2  ^ ,„,2,6B„2  ^ ♦ t„,2,6B„2 


m 


8>/2>r>/3  B„ 


^^Hf+3fo)2/6B^2] 

[^-(f  - 3f«)2/6B„2  ^ - fo)‘/6B,2^  ,^)2/,b„2 

^^^f+3f„)2/6B„2] 


^ 3^-(f  + fo)2/2(2B^2  ^ b„2)^  ^ 3t„)2/2(2B„2  . b„2,] 


+ -— yy  V [e'<^  ■ ^‘'o>^/2(B^2  ^ 2B„2)  <f . r^)2/2,B,„2  + 2B„2 , 

8/2»JB„2  + 2B  2 


“'m  *'*'n 


+ ''o)^/2(Bm2  + 2B„2)^  + 3I„|2/2(B  J + 2B„2,] 


(7-10) 


Substltudnt;  ctiuation  (7-10)  In  equation  ( 7-'))  and  extracting  ordy  the  low-lrequeney  eomponenls 
gives,  lor  the  detector  oiilpiil  low-pass  spectrum  with 


“ ’"  isf.  |T;7  “B;r  ' I 

2 ■>. 


J:iL  (^/’/2B„,2/.’£^.  ''/2B,r 

^ \ B|n 


7-.S 


, l‘13  ( °m  ^ .|-/2B„- 

8vT2»  I n/5  Bn,*"  /3  B,,'"  B„ 

, -l-/4B„-l 

vT2  B„  ) 

-t-/4B„-  -I“/2B,,“| 

4/4^  1 B^  B„  B„  ) 


(7-1  I) 


From  equation  (7-1 1 ),  two  important  quantities  are  calculable.  First,  the  noise  tloor  is  taken 
to  be  that  value  of  the  spectrum  beyond  the  region  of  appreciable  multipath  energy.  In  the  above 
equation  this  consists  of  only  those  terms  with  exponentials  involving  Bj^“  since  for  f » the 
other  terms  are  negligible.  The  noise  floor  can  thus  be  calculated  as  a function  of  C/J^^  and  S/I  by 
evaluating  the  various  coefficients  and  setting  the  remaining  exponential  multipliers  to  unity.  This 
calculation  is  plotted  in  figure  7-2,  where  the  ratio  of  mean  value  squared  to  noise  floor  is  plotted. 
Fortunately,  as  discussed  below,  the  ratio  is  largely  insensitive  to  S/I  for  the  parameter  region  of 
interest. 

A second  quantity  of  interest  derived  from  equation  (7-1  I ) is  the  low-pass  variance,  mathe- 
matically defined  as  the  spectral  integral  of  (7-1 1)  multiplied  by  the  low-pass  filter  transfer  function 
squared.  The  postdetection  low-pass  filter  has  an  equivalent  one-side  bandwidth  of  800  Hz.  so  it  may 
be  assumed  that  all  terms  having  bandwidth  related  to  Bj^  are  passed  unattenuated,  while  the  terms 
involving  bandwidth  B^^  may  be  considered  as  having  a constant  spectral  density  over  the  width  ol 
the  low-pass  filter.  The  variance  of  the  low-pass  filter  output  then  becomes; 


"lp‘ 


''02‘°m'‘  . , 2 2 . ''l3^Si^  . ''22^° 

^ ■•■‘Ml  °m  ■•■ ■•■ ^ 


2„  4 
m 

4 


2B,-  (/2an\4o„,o„^hQ2--  h,,V  2/ -fill  ^ ^ \ 

Bn  lUsA’t  4NT4tt  J 2 J27T  hJ'Iv  8^2 tr /2  / 

22  \4J47r  4j4n  )] 


(7-12) 


The  mean  value  squared  of  the  low-pass  output  is  simply  The  ratio  of  mean  value 

stjuared  to  variance  can  be  evaluated  as  a function  of  C/N^.,  and  S/I.  since  they  jointly  constrain  the 
relationship  between  and  o^“  in  Fquation  (7-12).  This  ratio  is  shown  in  figure  7-3  as  a function 
of  f’/Ny  with  S/I  varying. 
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Note  that  at  high  C tlie  ratio  is  asymptotic  to  a value  detemiined  by  tlie  value  of  S/1, 
and  this  asymptote  is  an  increasing  function  of  S/I.  On  the  other  hand,  at  very  poor  C’/N^.  additive 
noise  dominates  the  multipath  and  a common  asymptote  is  approached  for  all  S/I. 

A 

•Mso  shown  in  figure  7-.^  are  performance  data  taken  in  the  laboratory  lor  a no-multipath  ease 
(S/l  =®  ).  Ilte  experimental  data  confirm  the  above  cw  theory  extremely  well. 

It  is  apparent  from  ligure  that  the  observability  of  the  S/I  parameter  is  very  poor  at  lower 
C'/Ny  values.  This  results  since  the  additive  noise  contribution  to  the  variance  dominates  that  of  the 
multipath,  and  the  curves  tend  to  converge  at  low  C/N^  and  become  independent  of  S/I.  Conse- 
quently. a modified  variance  calculation  was  employed  that  used  a 2.S0-HZ  rectangular  characteristic 
(this  liltering  was  performed  in  the  frequency  domain  by  operating  on  the  spectral  estimates).  Again 
it  may  be  assumed  that  all  multipath-related  components  are  being  passed  unattenuated,  while  the 
noise  contribution  to  the  variance  is  reduced  by  a factor  of  roughly  2.50/800.  Recalculation  of 
equation  (7-12)  using  Bj-  = 250  Hz  then  provides  a theoretical  relationship  between  the  square  of 
the  process  mean  divided  by  the  new  variance  and  the  parameters  C/N^  and  S/I.  This  relationship 
is  shown  graphically  in  figure  7-4. 


7.2  KNVFLOPK  DHTI  CTOR  OUTPUT  PROCESSING  PROCEDURE 

fhe  analog  waveform  is  sampled  at  the  Test  Data  Processing  Center  at  a 2000-Hz  rate  using  a 
lO-bit  A/D  quantizer  ( 1 024  levels).  Little  aliasing  distortion  is  thus  present.  For  every  observation 
time.  1025  such  samples  (~t).5  sec)  are  written  on  magnetic  tape  along  with  the  time  code. 

Processing  of  this  time  series  consists  of  the  following  steps,  illustrated  in  figure  7-5. 

a.  Estimate  mean  and  mean-square  value. 

b.  Remove  mean  and  linear  drift  components,  if  any,  from  the  time  samples.  This  removes 
the  dominance  ol  the  large  dc  component  in  the  neighboring  spectral  estimates. 

c.  laper  the  time  series  at  each  end  of  the  record  according  to: 

i 

1 

tr(i-I)  j 

S(i)  = S(i)  sin  , I 5i  <0.  IN 

0.2N  ( 


n(  N-i ) 

= S(i)  sin , ().‘>N<i«N 

0.2N 

= S(i)  elsewhere.  (7-I.D 


7-») 


32  36  40  44  48  52  56 

C/Nq,  dB-Hz 


Figure  7-4.  Ratio  of  Squared  Mean  to  Variance  — Low-Pass  Filter  Output  of  Envelope 
Detector,  Bf  = 250  Hi 
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Figure  7-5.  Envelope  Detector  Data  Analysis  Algorithm 
7-1  I 


This  cosine  tapering  acts  to  convolve  the  observed  random  process  with  a “filter”  having 
smoother  impulse  response;  hence  the  equivalent  spectral  window  has  lower  side  lobes. 
This  minimizes  spectral  bias  due  to  strong  components  "looking  througli”  the  side  lobe 
to  contribute  at  other  I'requencies. 

d.  Compute  the  discrete  Fourier  transform  of  the  resultant  record  and  obtain  the  magnitude 
of  the  spectrum.  With  the  data  set  given,  spectral  coefficients  are  available  roughly  every 
1 Hz  from  dc  to  1000  Hz. 

e.  Smcxjth  the  estimate  in  the  frequency  domain.  The  estimate  from  step  d typically  has 
rather  large  variance.  Estimates  of  smaller  variance  are  obtainable  by  weighting  several 
adjacent  coefficients  on  either  side  of  a desired  frequency.  In  the  procedure  given  here, 
uniform  weighting  was  used  for  10  points  on  either  side  of  the  point  in  question. 

f.  Compute  the  following  statistics: 

1 ) Noise  lloor  the  average  spectral  level  between  2.'i0  and  bOO  Hz.  a region  nominally 
devoid  of  multipath  and  within  which  the  bandpass  and  low-pass  filter  rollotf  are 
not  yet  important. 

2)  Multipath  power  - the  spectral  energy  lying  between  2 and  2.‘!0  Hz.  This  is  the 
modified  variance  described  previously,  and  which  also  includes  some  noise 
contribution. 

3)  (Mean)“/noise  fioor  and  (meanj'/multipath  power  - these  terms  are  used  to 
estimate  C/N^  and  S/I. 

g.  Using  (meanj-Znoise  lloor.  determine  C/N^  from  Figure  7-2  using  the  S/1  = KidB 
curve.  At  this  point  S/1  is  an  unknown  parameter,  but  little  error  is  incurred  over  the 
range  of  38  dB-Hz  < C/N^  < 48  dB-Hz  by  using  the  Sri  = lb  dB  curve. 

h.  Using  this  same  C/N^  and  ( mean )~/multipath  power,  interpolate  for  S/I  in  figure  7-4. 


7.3  PRfXiRAM  INPUT/OUTT’UT  l)FS(  RIP'ITON 

T he  various  programs  using  the  envelope  detector  analysis  algorithm  have  the  s;iine  signal 
strength  data  outputs  for  all  programs,  with  the  minor  exception  ol  antenna  test  cases  fhe  input 
data  to  this  analysis  algorithm  is  the  same  for  all  programs. 


Input  data  specifiod  whenever  the  signal  strength  analysis  algorithm  is  used  are  as  follows: 


Sample  rate: 

Taper  rate: 

Spectral  window  width: 

Initial  values  of  U|.  a->,  a^^: 

Low-pass  filter  bandwidth: 

Bias  value: 

Signal  sta-ngth  record  analysis  skip  factor; 
Signal  strength  plot  skip  factor: 


: kHz 
0.1 

40  Hz 

Specified  by  user 
800  Hz 

Specified  by  user 
Specified  by  user 
Specified  by  user . 


Direct  output  from  the  C'DC  bbOO  computer  is  in  the  form  of  computer  listings.  This  listing 
includes  a printout  of  the  smoothed  spectnim  as  estimated  from  the  signal  strength  time  data,  con- 
vergence data  used  in  the  least-sijuares-fit  algorithm,  and  a printout  of  the  fitted  model  spectrum.  It 
also  includes  the  following  parameters  as  a tabulation: 


Date  and  time  of  test 

Mean  value  of  the  signal 

Mean  square  value  of  the  signal 

Theoretical  spectrum  coefficients  (a  I . a->.  a^l 

Ratio  of  squared  mean/a^  (noise  fioor) 

Variance 

Standard  deviation 

Ratio  of  squared  mean/variance 

MP.  defined  as  multipath  power  plus  noise  power 

Ratio  of  squared  mean/ multipath  power 

(■'No  . 


A signal  strength  summary  plot  can  be  made  for  each  of  the  summaries  described  above.  The 
data  to  be  plotted  are  output  onto  a plot  tape,  which  is  then  run  off-line  on  the  SC  4020.  The  plots 
include  the  measured  data-smoothed  spectral  estimate  and  the  fitted  model  spectrum.  In  addition, 
several  parameters,  also  printed  on  the  listing,  are  included  on  the  bottom  of  the  plot  output. 


7.4  SAMPl  h Rl  SULTS 

Kigures  7-b  and  7-7  illustrate  the  smoothed  spectral  estimates  for  two  cases  selected  arbitrarily, 
rite  first  is  for  a Type  I test  when  no  appreciable  multipath  is  anticipated.  (Actually,  the  spectral 
analysis  approach  shows  an  absence  of  additional  energy  in  the  0-100  Hz.  region  for  Type  I tests,  indi- 
cating that  the  antenna  S I is  at  least  I ^ dB  m general. ) The  secoml  plot  (fig.  7-7)  is  for  a Type  II  run 
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with  intended  S/1  in  the  8-dB  range.  A definite  spectrum  increase  is  observed  in  the  multipath 
Doppler  zone,  and  the  curve  fitting  and  parameter  estimation  procedures  are  illustrated  below. 


[ : 


From  figure  7-7  the  squared  mean/noise  fioor  statistic  is  52,500,  and  from  figure  7-2  the 
estimated  C/N^  is46.9dB-Hz.  Although  the  squared-mean-to-modified-variance  ratio  (250-Hz 
bandwidth)  is  not  printed  on  figure  7-7,  it  is  given  in  computer  listings  and  for  this  case  was  1 1 .0. 

Use  of  figure  7-4  then  yields  an  S/1  of  6.5  dB. 

Repetitive  averaging  of  similar  cases  over  a data  run  provides  estimates  with  accuracy  ascribed 
to  be  +0.5  dB  in  C/N^  and  ±1  dB  in  S/1.  Table  7-1  is  an  example  printout,  showing  smoothed  spec- 
trum, best-fit  spectrum,  and  statistical  summary.  The  data  interval  does  not  correspond  with  the 
plots  of  figures  7-6  and  7-7. 

As  a check  on  the  agreement  between  hardware  C’/Np  measurements  (for  example  those 
obtained  with  a wave  analyzer  or  spectrum  analyzer)  and  those  generated  with  the  computer  proce- 
dure, a sample  tape  was  generated  for  hardware-determined  C/Nq  values  of  38  to  50  dB-Hz.  This 
tape  was  digitized  and  processed  in  the  usual  manner,  and  the  calibration  curve  of  figure  7-8  was 
generated.  This  curve  was  not  used  in  the  data  reduction  other  than  to  show  the  rather  good  agree- 
ment between  hardware  and  software  approaches.  The  computer  method  becomes  less  accurate  below 
38  dB-Hz  because  the  noise  bandwidth  at  the  detector  input  is  5 kHz,  implying  that  the  output 
statistics  are  quite  insensitive  to  C/N^  values  lower  than  this  point.  For  example,  figure  7-2  is  nearly 
horizontal  in  this  region. 
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Figure  7-7.  Spectra!  Analysis  Plot  of  Envelope  Detector  Output,  S/I  = 6.5  dB 
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Computed  C/N„,  dB  Hz 


Figure  7-8.  Computed  Versus  Real-Time  C/N^ 
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8.  VOICE  MODEM  DRandA  PROCESSING 


Voice  data  processing  includes  parallel  efforts  by  Boeing  and  CBS  Laboratories.  CBS  Labora- 
tories provide  intelligibility  scoring  using  listener  panels  and  supply  test  results  to  Boeing.  Boeing 
processes  the  data  to  obtain  C/N^  and  S/L  described  below. 


8.1  DETERMINATION  OF  C/N^  AND  S/I 

The  signal  strength  data  are  computer  analyzed  to  detennine  C/N^  and  S/I  for  voice  modem 
tests.  Processing  is  controlled  by  a program  that  calls  subroutines  to  position  the  data  tape  by  using 
file  numbers  and  segment  times.  The  algorithm  is  shown  in  figure  8-1.  The  determination  of  C/N^^ 
and  S/I  uses  algorithms  described  in  section  7. 

The  digital  tape  format  is  as  previously  described  in  figure  b-3.  Each  signal  strength  record 
consists  of  four  time  words,  one  identification  word,  a 1 02.5  signal  strength  sample  words.  Each 
signal  strength  word  is  I 2 bits  long.  A record  occurs  i y I .S  sec  as  controlled  by  the  IRIG-B  time 
code. 


8.2  INPUT  SPECIFICATIONS  AND  OUTPUT  FORMATS 

The  input  to  voice  program  test  segments  consists  of  the  entire  specification  of  signal  strength 
analysis  related  parameters  as  described  in  section  7.4.  Further  data  specified  for  voice  segments  that 
relate  to  tape  positioning  and  segment  identification  include: 

Date  of  test 
Tape  number 
Tape  file 

Start  and  end  time  of  segment 
Test  mode  (1  ype  1 or  II). 

Ilie  output  from  a voice  analysis  segment  is  described  in  section  7.4  and  consists  entirely  i>l  data 
resulting  from  signal  strength  analysis. 
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Figure  8- 1 . Voice  Data  Processing  Algorithm 
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9.  DIGITAL  COMMUNICATION  DATA  MODEM  DRandA  PROCESSING  ' 


Digital  data  modem  processing  is  controlled  by  a program  that  calls  subroutines  to  position  the 
tape  at  the  correct  file  number  and  time,  process  the  signal  strength  data  to  give  C/N^  and  S/1 
statistics,  and  process  the  PN  sequence  digitial  data  message. 


9.1  DIGITAL  TAPE  FORMAT 

The  digital  tape  format  for  digital  segments  is  as  previously  shown  in  figure  6-6.  This  figure 
shows  that  each  digital  data  record  is  comprised  of  four  time  words  and  one  identification  word,  fol- 
lowed by  five  blocks  of  data  bits.  Each  block  (one  per  modem)  contains  205  1 2-bit  words  (2460  bits 
per  block)  and  corresponds  to  a digital  data  channel  recorded  on  the  airplane,  thus  corresponding  to 
approximately  2 sec  of  data  at  1 200  bps.  Signal  strength  is  sampled  every  seven  digital  data  records 
(approximately  14  sec)  and  is  formatted  with  the  applicable  time  code  and  ID  shown  in  figure  6-3. 


9.2  ALGORITHM  FOR  DIGITAL  DATA  PROCESSING 

Processing  the  digital  data  involves  signal  strength  parameter  computation,  bit-error  rate 
determination,  and  error  tabulation.  Figure  9-1  illustrates  the  sequence  of  operations  performed. 

For  the  time  segment  defined  for  analysis,  data  is  input  and  subroutines  are  initialized. 

Signal  strength  records  are  separated  from  digital  data  records  and  are  processed  to  provide  C/N^  and 
S/I  parameters  as  described  in  section  7.  The  digital  data  records  are  processed  as  described  in  the 
next  section. 


9.3  DIGITAL  DATA  RECORD  PROCESSING  ALGORITHM 

The  digital  data  record  processing  algorithm  is  shown  in  figure  9-2.  Within  each  data  record, 
the  five  blocks  of  data  lone  per  recorded  modem  channel)  are  processed  in  turn  starting  with  channel 


9-1 


Figure  9-2.  Digital  Record  Processing  Algorithm 


1.  SiiK-e  not  all  five  diannels  are  always  active,  a skip  provision  is  made  lor  omitting  the  processing  lor 
selected  channels.  If  the  data  record  to  be  processed  is  the  first  record  within  a time  segment,  data 
synchronization  (reference  PN  sequence  synchronization  I is  performed.  I his  process  will  be  described 
in  detail  in  section  1.  If  the  data  record  is  not  the  l1rst  within  a lime  segment,  it  is  assumed  that 
synchronization  has  been  previously  achieved  and  bit-error  counts  are  made.  The  bit-error  count  is 
then  subjected  to  an  acceptance/rejection  test.  If  the  count  e.xceeds  a lO't  error  rate  threshold,  the 
resynchronization  mode  is  entered.  If  the  error  count  is  accepted,  the  nuitiber  ot  errors,  the  error 
locations,  and  other  data  are  tabulated.  Inter-error  spacing  and  block  error  histograms  are  also  com- 
puted and  stored  for  the  segment  summary  output. 

9.3. 1 Synchronization  of  Input  Data  to  Reference  Data 

To  process  the  digital  data,  two  items  associated  with  synchronization  must  be  accomplished: 

a.  I he  error-free  reference  data  sequence  must  be  generated. 

b.  The  reference  sequence  must  be  synchronized  in  phase  with  the  received  data  lor 
bit-by-bit  comparison. 

The  transmitted  PN  data  sequence  is  of  length  2047.  described  by  the  generator  polynomial 
x'  * -t-  X-  + 1.  An  identical  reference  PN  sequence  is  generated  in  the  computer.  Starting  with  each 
bit  of  the  reference  sequence.  2047  bO-bit  words  are  formed  and  stored  in  memory.  Because  of  the 
properties  of  a PN  sequence,  each  of  these  2047  60-bit  words  has  a unique  I I -bit  pattern  in  the  I 1 
most  significant  bit  positions.  The  remaining  49  bit  positions  are  filled  with  the  normal  continuation 
of  the  PN  sequence.  A cross-reference  table  is  created  so  that,  given  a particular  1 1-bit  pattern,  the 
60-bit  word  with  the  same  I 1-bit  pattern  in  the  most  significant  bit  positions  can  be  found. 

The  algorithm  for  synchronization  is  illustrated  in  figure  9-3.  Basically  the  method  consists  ol 
extracting  the  first  I 1 bits  from  the  input  string  and  using  the  cross-reference  table  to  determine  which 
of  the  2047  60-bit  reference  words  has  the  same  I 1-bit  pattern  in  the  first  I 1 bits.  1 he  reference 
words  are  then  called  in  the  proper  order  to  form  the  reference  se(|uence.  Bit-by-bit  comparisons  are 
then  made  with  the  input  data  string  for  the  entire  block  of  data  bits  m the  record  lor  that  channel. 

If  the  two  strings  are  not  synchronized,  the  error  rale  will  be  approximately  .‘'04  and  Ihe  sync  posi- 
tion will  be  rejected  by  the  error  rate  threshold  test  ( lO'i  errors)  th.il  follows.  If  rejection  occurs.  Ihe 
synchronization  algorithm  is  repeated,  this  time  using  the  second  group  of  I I bits  from  Ihe  input  data 
stream.  I his  process  is  continual  until  synchronization  is  achieved  or  the  end  of  Ihe  channel  block  is 
reached.  If  synchronization  is  not  achieved  (all  trials  yield  error  rales  greater  Ih.m  0.  lOl,  the  24(i()  bits 
Ifom  that  particular  recoril  and  channel  block  are  rejected.  When  Ihe  eiul  ol  the  d.ila  block  lor  this 
channel  is  reached.  Ihe  synchronization  .ilgorilhm  is  applied  to  Ihe  d.it.i  block  lor  Ihe  next  ch.inncl. 
Synchronization  is  performed  independently  for  each  ol  Ihe  loe  ch.innels  since  bits  Irom  the  various 
demodulators  are  not  necessarily  identically  limed. 


START 


SET  POINTER 
TO  FIRST  11  BITS 
OF  INPUT  STRING, 
SET"SYNCHRO'' 

TO  NO 


TEST  FOP 
END  OF 
•SEARCH. 


RETURN 


NER  < MAX 


SET  SYNCHRO" 
TO  YES 
SET  END  OF 
SEARCH 


MOVE  POINTER  TO 
NEXT  11  BITS  IN  THE 
INPUT  STRING 
SET  END  OF  SEARCH 
IF  END  OF  STRING 


Figure  9-3.  Digital  Message  Synchronization  Algorithm 
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It  is  readily  apparent  that  an  error-free  group  of  I I bits  in  the  input  data  string  is  required  for 
synchronization  to  be  achieved  with  the  above  method.  At  the  bit-error  rates  encountered  in  this  test 
program,  this  condition  has  a high  probability  of  being  met  with  the  first  group  of  1 1 bits  read.  Even 
if  several  attempts  are  required  to  achieve  synchronization,  all  2460  bits  in  a channel  block  are 
examined  when  bit-error  counts  are  made.  This  is  achieved  simply  by  using  the  reference  words  to 
construct  the  reference  PN  sequence  for  the  entire  record,  both  forward  and  backward  from  the 
synchronization  point. 


To  minimize  computer  time  requirements,  resynchronization  is  not  performed  unless  required. 
Once  synchronization  has  been  achieved,  the  position  of  the  last  bit  in  each  channel  is  saved  from  one 
record  to  the  next,  thus  making  it  possible  to  resume  comparing  the  reference  string  with. the  input 
string  without  resynchronizing.  The  reference  sequence  and  the  2047  60-bit  reference  words  stored 
in  memory  are  generated  only  once  for  each  computer  processing  run. 


9.3.2  Block  Error  Statistics 

Block  error  statistics  and  inter-error  spacings  provide  a measure  of  the  channel’s  tendency  for 
error  bursting.  To  provide  block  error  statistics,  the  input  data  of  one  channel  is  considered  as  a single 
string,  from  start  to  end  of  an  entire  time  .segment.  The  string  is  divided  into  adjacent  blocks  of  n bits 
(n  = 24).  The  number  of  errors  in  each  block  is  counted  and  tabulated  to  give  the  frequency  of  occur- 
rence for  a given  number  of  errors  per  24-bit  block.  This  tabulation  is  output  as  part  of  the  segment 
summary. 


9.3.3  Inter-Error  Spacing 

The  histogram  ot  inter-error  spacing  gives  the  frequency  of  a given  separation  in  bits  between 
successive  errors.  As  for  the  block  error  statistics,  the  input  data  is  considered  to  be  a single  string 
for  a given  channel  and  time  segment. 


Figure  9-4  illustrates  the  received  string,  with  some  error  positions  (x)  and  the  inter-error 
spacing  that  is  given  by  the  difference  in  error  positions.  The  tabulated  output  gives  the  frequency  of 
occurrence  for  a given  inter-error  spacing  and  is  included  as  part  of  the  segment  summary. 
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Figure  9-4.  Inter-Error  Spacing  Definition 


9.4  PROGRAM  INPUT  AND  OUTPUT  DESCRIPTION 

Inputs  to  the  digital  data  processing  program  include  all  the  signal  strength  processing  param- 
eters explained  in  section  7.4.  In  addition,  the  following  parameters  are  specified  to  position  the 
tape,  identify  the  test,  and  specify  the  channels  to  be  analyzed: 

Data  of  test 
Tape  number 
Tape  file 

Start  time  of  segment 
End  time  of  segment 
Test  mode  (Type  I or  II) 

Channels  to  be  analyzed. 

Outputs  provided  by  the  digital  data  analysis  are  as  follows. 

a.  Error  locations  are  tabulated  for  each  record  processed  ( table  9- 1 ).  The  time  of  each 
record  is  printed  and,  if  errors  are  present,  the  channel  number  containing  errors,  number 
of  errors,  and  bit  locations  of  the  errors  in  the  data  block  (locations  run  from  bit  1 to 
2460)  are  also  given. 

b.  At  the  end  of  each  time  segment,  a segment  summary  is  provided  (table  9-2).  This 
summary  includes  the  date,  start  and  end  time  of  segment,  experiment  mode  identifica- 
tion, tabulation  of  number  of  bits  analyzed  for  each  of  five  channels,  number  of  errors, 
and  bit-error  rate.  Also  included  is  the  block  error  histogram  for  each  channel  for  the 
segment  and  the  inter-error  spacing  histogram  (table  9-.^).  At  the  bottom  of  the  inter- 
error histogram,  (he  number  of  resynchronizations  required  during  the  segment  is  given 
for  each  channel. 
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TABLE  9-2.  DIGITAL  SEGMENT  SUMMARY  AND  BLOCK  ERROR  HISTOGRAM 


TABLE  9-3.  INTER-ERROR  SPACING  HISTOGRAM 


10.  RANGING  MODEM  DRandA  PROCESSING 


Analysis  of  TSC  and  NASA  ranging  data  is  provided  by  the  same  ranging  software.  The  analy- 
sis, performed  on  time  segments  up  to  5 min  in  length,  includes  both  relative  ranging  error  and  signal 
strength  analysis.  The  time  segments  analyzed  are  chosen  for  their  relatively  constant  aircraft  param- 
eters and  for  the  constancy  of  C/N^  and  S/I. 


10.1  DIGITAL  TAPE  FORMAT  (TSC  RANGING) 

Source  tapes  are  converted  to  digital  tape  for  the  TSC  ranging  modem  at  TDPC.  NASA  infor- 
mation is  received  in  digital  form.  The  digital  tape  format  for  the  TSC  digital  ranging  modem  data  was 
shown  in  figure  6-4.  Table  10-1  describes  the  meaning  of  each  of  the  status  bits  and  code  information 
bits.  Each  bit  in  the  range  information  corresponds  to  25  nsec  in  range  time.  For  each  range  value 
there  is  a corresponding  time  and  for  every  16  range  readings  there  is  a signal  strength  record. 


TABLE  10-1.  TSC  DIGITAL  RANGING  TAPE  LEGEND 


BIT 

1 THRU  27 

0 filler 

BIT 

28 

70  MHz  PLL  LOCKED  + 

BIT 

29 

70  MHz  PLL  LOCKED  - 

BIT 

30 

CLOCK  PLL  LOCKED 

BIT 

31 

CORRECT  CORRELATION 

bit 

32 

... 

BIT 

33 

CODE  FORMAT  1 

BIT 

34 

CODE  FORMAT  3 

BIT 

35 

CODE  FORMAT  5 

BIT 

36 

CODE  FORMAT  7 

BIT 

37 

CODE  FORMAT  2 

BIT 

38 

CODE  FORMAT  4 

BIT 

39 

CODE  FORMAT  6 

BIT 

40 

CODE  FORMAT  8 

BIT 

41 

MSB 

RANCE  DATA 
1 BIT  = 25 

BIT 

60 

LSB 

lO-l 


10.:  RANGING  DATA  ANALYSIS  ALGORITHMS 


The  overall  ranging  data  processing  block  diagram  is  given  in  figure  10-1.  The  signal  strength  1 

records  are  separated  from  the  ranging  data  records  and  are  analyzed  as  described  in  section  7 to  yield 
estimates  of  C/N^  and  S/1.  All  processes  described  in  that  section  apply  to  the  ranging  data  signal 
strength  samples. 

The  major  blocks  of  the  ranging  analysis  algorithm  are  shown  in  figure  10-2.  These  blocks  are 
described  below. 

10.2.1  Purge  of  Unreasonable  Range  Values 

Before  the  range  values  are  analyzed,  the  entire  array  of  range  values  is  examined  to  determine 
if  the  values  are  reasonable.  Specifically,  the  second  range  value  is  initially  assumed  to  be  correct  and 
any  change  from  this  range  to  the  range  at  the  next  data  time  will  be  less  than  some  specified  threshold 
value.  This  test  is  repeated  on  the  entire  array  with  the  differential  in  range  being  calculated  for  all 
adjacent  ranging  values.  For  acceptance  , 

lr(tj+,)-r(ti)  |<  Rp  (10-1) 

where: 

r(tj)  = range  at  time  tj 

r(tj^.j ) = range  at  the  next  range  reading  output  time 

Rp  = differential  range  threshold  that  may  not  be  exceeded  (Rp  usually  10,000  meters). 

The  values  of  range  that  do  not  meet  this  criterion  are  purged  from  the  ranging  data  array  and 
the  rejected  data  is  printed  out.  When  a range  value  is  rejected,  the  last  valid  value  of  range  is  saved 
and  used  as  r(tj)  until  a point  is  found  that  meets  the  above  criterion. 

10.2.2  Least-Squares  Fit  and  Error  Array 

Values  of  time  and  relative  range  are  read  from  the  digital  tape  into  an  array  in  the  computer. 

Following  the  initial  purge  of  unreasonable  values,  a least-squares  fit  to  the  empirical  data  is  made 
using  a second-order  curve; 

r(t)  = at-  + bt  + c.  ( KV2) 

ITie  method  of  the  least-squares  fit  is  explained  in  appenilix  A. 


II 

j 


10-2 


When  the  least-squares-t'it  curve  has  been  determined,  an  error  array  is  formed  using  the  least- 
squares  curve  as  the  assumed  actual  path.  The  array  is  formed  using  the  following  algorithm. 


ArU;)  = ri(tj)-r(tjl.  (10-3) 

where: 

rj(tj)  = real-time-measured  range  at  time  tj 

r(tj)  = assumed  actual  range  at  time  tj  as  estimated  by  the  least-squares-tit  curve  trotn 
equation  ( 1 0-2). 


I 

\ 

) 

I 
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10.2.3  Krror  Statistics 

Error  analysis  is  performed  on  the  error  array  formed  above.  The  mean  and  mean  squared 
values  of  the  error  array  are  calculated  as; 


n 


(10-4) 


( 10-5) 


where  n is  the  number  of  range  values  in  the  array. 

The  variance  and  standard  deviation  are  then  calculated  from 

= Ar  - (Ar)  (10-(iI 

A second  purge  is  now  performed  on  the  data.  Any  points  found  outside  the  5o  limits  are 
deleted  since  they  are  assumed  to  be  caused  by  equipment  malfunction  or  incorrect  ambiguity  resolu- 
tion and  would  corrupt  the  processed  outputs.  I he  purge  is  accomplished  by  eliminating  the  ranging 
values  from  the  data  corresponding  to  points  outside  the  5a  limits.  Purged  events  are  printed  as  part 
of  the  output.  If  this  second  purge  operation  occuis,  the  least-squares  fit  is  repeated  to  form  a new 
error  array  and  new  values  are  computed  for  mean,  mean  square,  and  variance. 
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10.2.4  Error  Distribution  Histogram 


The  error  distribution  histogram  is  formed  as  part  of  the  analysis  that  leads  to  a chi-squared 
goodness-of-fit  test.  The  lower  limit  for  the  frequency  distribution  is  chosen  to  be  - -o  and 

the  upper  limit  is  Xy  = + 2a.  Typically,  eight  bins  are  used  in  computing  the  lustogram.  The 

midpoint  of  each  bin  is  found  by 


AX 

Xl-— 


AX 


AX  AX 
. Xii  + 

-l  b -1 


where: 

A.X  ( X - .X  ^ t ; ( K - 2 ) 

K = number  of  intervals  or  bins . 

In  forming  the  actual  error  histogram,  the  two  tail  intervals  cotitain  all  the  points  outside 
the  upper  and  lower  frequency  distribution  limits.  The  remaining  error  points  are  contained  within 
their  proper  intenal.  between  the  upper  and  lower  limits  of  (he  liistogram  as  determined  above. 


10.2.5  Chi-Squared  Cioodness-of-Fit  Test 


.A  chi-squared  goodness-of-fit  test  is  made  on  the  range  error  data.  The  mean  and  standard 
deviation  of  the  observed  error  array  are  used  as  estimation  parameters  to  detennine  a theoretical 
Gaussian  distribution  with  the  same  total  number  of  points  as  the  observed  sample  of  range  error 
values.  This  theoretical  distribution  is  used  to  construct  an  expected  observations  histogram  having 
the  same  number  of  cells  and  bin  size  as  the  error  histogram  constructed.  From  these  two  histograms, 
the  following  test  statistic  is  computed: 


X 


2 


(0:  ■ 


(10- ■7) 


where: 

K = number  of  bins 

Oj  = number  of  observed  ranging  errors  in  sample  bin  i 

Fj  = theoretically  expected  number  ol  error-,  within  sample  bin  i. 
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The  test  statistie  is  used  in  combination  with  the  number  of  degrees  of  freedom  to  deter- 
mine the  probability  that  a sample  of  the  same  size,  randomly  chosen  from  a known  normally  distri- 
buted population  of  large  size  and  having  the  same  mean  and  standard  deviation,  will  have  a 
statistic  greater  than  or  equal  to  that  of  the  observed  sample. 

The  number  of  degrees  of  freedom  is  known  to  be  equal  to  K - k,  where  k = I + number  of 
parameters  estimated  from  the  observed  data.  Since  mean  and  standard  deviation  are  both  estimated, 
^ ~ ^ +2-3  and  the  degrees  of  freedom  are  found  to  be  K - 3.  The  probability  values  are  deter- 
mined by  a look-up  table  contained  in  a subroutine  known  as  the  chi-square  distribution  percentage 
point  table. 


10.3  INPUT  SPECIFICATION  AND  OUTPUT  FORMAT 

The  input  to  ranging  program  test  segments  consists  of  the  entire  specification  of  signal 
strength  analysis  related  parameters  (described  in  sec.  7.4).  Further  data  specified  for  ranging 
segments  relate  to  tape  positioning  and  segment  identification: 

Date  of  test 
Tape  number 
Tape  file 

Start  time  of  segment 
End  time  of  segment 
Test  mode  (Type  I or  II). 


The  output  Irom  the  CDC  6600  consists  of  ( I ) signal  strength  summaries,  plots,  and  time 
segment  information  discussed  in  section  7.4  and  (2)  range  data  information.  The  range  data  informa- 
tion outputs  include  the  following: 


Input  data  read  from  the  digital  ranging  test  tape,  as  shown  in  the  example  of  table  10-2. 
The  items  listed  are: 


I)  Point: 


Number  n.  identilying  n^^^  measurement  of  segment  being 
analyzed 


2)  Status  bit: 


1 70  MHz  PEL  LOCKED  + 

2 70  MHz  PLL  LOCKED - 

3 CLOCK  PLL  LOCKED 

4 CORRECT  CORRELATION 

5 — 
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3)  Time:  GMT  at  which  range  measurement  was  made. 

4)  Range:  Range  in  meters 

5)  Range  in  nsec:  Range  in  nanoseconds 

b.  Purged  range  errors,  censored  range  error  values,  and  range  error  array  (output  not 
shown). 

c.  The  ranging  segment  summary,  as  shown  in  table  10-3.  The  mode  type  refers  to  Type  1 or 

: 11  and  the  code  status  bits  refer  to  table  10-1.  Other  items  are  self-explanatory. 

1 

i 

I 

10.4  PLACE  RANGING  DATA  ANALYSIS 

PLACE  ranging  tapes  recorded  and  furnished  by  NASA  for  analysis  are  reformatted  from  nine- 

■ to  seven-track  digital  tape  compatible  with  the  Boeing  CDC  6600  computer.  The  slot  and  row  data 

■ pertaining  to  the  KC-135  are  used  to  determine  which  records  are  to  be  analyzed.  Time  and  range  are 

: read  from  the  chosen  records  and  put  into  storage  in  memory.  At  the  end  of  the  segment  the  range 

values  that  have  been  stored  in  memory  are  analyzed  using  the  algorithms  described. 

The  ranging  data  analysis  outputs  are  essentially  the  same  in  format  as  those  described  in 
section  10.3  except  that  printouts  specific  to  the  TSC  digital  ranging  modem  and  signal  strength 
analysis  are  not  included.  An  example  output  is  given  in  table  10-4. 
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TABLE  10-4.  RANGING  SEGMENT  SUMMARY,  PLACE  MODEM 
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1 1 ANTENNA  TEST  DRandA  PROCESSING 


The  antenna  test  data  is  processed  in  the  same  manner  as  the  voice  data  with  one  minor 
deviation;  it  is  identified  as  antenna  data  in  the  identification  word  as  shown  in  figure  (>-3.  The 
signal  strength  records  described  below  are  examined  for  the  I or  0 antenna  bit  to  determine  to 
which  antenna  the  signal  strength  record  applies.  The  algorithm  (low  is  essentially  the  same  as  for 
the  voice  case  (fig.  8-1 ).  The  signal  strength  data  is  processed  to  obtain  C'/Nq  and  S/I  using  the 
algorithm  described  in  section  7. 


I l.l  DIGITAL  TAPE  FORMAT 

The  digital  tape  format  for  antenna  tests  was  shown  in  figure  b-3.  Each  signal  strength  record 
consists  of  four  time  words  (hours,  minutes,  seconds,  and  milliseconds),  one  identification  word,  and 
1025  signal  strength  sample  words.  The  most  significant  bit  of  each  signal  strength  sample  word  is 
either  I or  0 depending  on  the  receiver  to  which  the  carrier  detector  unit  is  connected.  The  number 
and  spacing  of  signal  strength  samples  are  under  the  control  of  the  receiver  multiplex  signal,  which  also 
determines  whether  the  signal  strength  bit  is  I or  0. 


11.2  INPUT  SPECIFICATION  AND  OUTPUT  FORMAT 

The  input  to  antenna  program  test  segments  includes  the  entire  specification  of  signal  strength 
analysis  related  parameters  as  required  (described  in  sec.  7.4).  Further  data  specified  for  antenna  test 
segments  relate  to  tape  position  and  segment  identification: 

Date  of  test 
Tape  number 
lape  file 

Start  time  of  segment 
End  time  of  segment. 

The  output  from  an  antenna  analysis  segment  is  as  described  in  section  7.4.  consisting  entirely 
ol  data  resulting  trom  envelope  detector  data  analysis.  The  antenna  bit  is  printed  on  the  envelope 
detector  analysis  output  to  indicate  to  which  antenna  the  data  is  applicable. 
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APPENDIX  A 

LEAST  SODARES-FIT  ALGORITHM 


The  least-si)uares-lit  method  is  used  several  times  in  the  processing  of  the  modem  evaluation 
data.  Althougli  this  method  is  well  known  and  documented,  some  explanation  is  necessary  to  describe 
the  software  used  in  one  of  the  ATS-b  processing  cases  (signal  strength  processing).  Solving  a least- 
squares  fit  consists  mainly  in  finding  the  minimum  of  a function  of  one  or  more  variables.  \n  the  case 
of  signal  strength,  a method  was  found  that  is  an  efficient  tradeoff,  insofar  as  implementation  is 
concerned,  between  a numerical  and  an  analytical  approach.  The  least-squares  method  used  in  modem 
evaluation  programs  is  presenteil  in  general  terms  first,  then  each  particular  case  is  explained. 


A.  I GENERAL  CASE 

Suppose  n measurements  (Xj,  Yj)  are  given.  The  least-squares  fit  is  a method  for  finding 

m unknown  parameter  values  of  an  analytical  model  Y = F(X.a|.a-, aj^,).  that  give  the  best  fit 

to  these  points  in  the  squared-error  sense.  The  fit  is  considered  best  when  an  “error”  El.  function  of 
Uj  parameters  and  input  points  Xj.  Y|  is  minimum.  E is  defined  as: 

n . 

""  2 (FlXj.aj.as a^)-Y|)T.  (A-l) 

i=l 

In  the  general  case  the  partial  derivative  of  F.  with  respect  to  aj  .aT a,.,.,  is  taken  and  set 

equal  to  zero.  The  resulting  m equations  are  then  solved  simultaneously  for  a|.as a^^. 

Specific  cases  of  least-squares  fit  used  in  the  analysis  of  ATS-b  data  are  discussed  individually 
in  the  following  sections. 


A.:  FITTINC;  Wini  a linear  FUN(  lion  EOR  linear  DRIFI  REMOVAl, 

The  function  F(X.aj)  chosen  is: 

F = a I + a 1 X . 
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Ilie  same  previously  described  steps  yield  linear  equations  with  three  unknowns,  which  are  solved  tor 
a I . a a 
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A.4  SPECTRUM  MODEL 

The  spectrum  model  used  is  not  linear  relative  to  the  aj  parameters,  a fact  that  complicates 
the  derivative  expression.  Short  of  an  analytical  solution,  we  use  a iterative  numerical  solution. 
Convergence  by  iteration  involves  computing  the  function  with  an  initial  set  of  values  for  aj  and 

varying  the  parameters  in  search  of  the  minimum  E(a|.a. a^^).  The  function  derivatives 

3E/3aj  provide  useful  information  to  increment  the  parameters.  Note  that  most  of  the  operations 
needed  to  compute  E are  common  to  the  computation  of  each  3E/3aj.  The  function  evaluation, 
derivative  computation,  and  the  convergence  algorithm  are  the  main  parts  of  the  least-squares  fit  in 
this  case. 


A.4. 1 Function  and  Derivative  Computation 
The  function  used  is; 
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This  may  be  rewritten  as 


F(X.aj)  = a,  K,  (ai.X)  + a3  Ft  (X) 


By  substitution  in  equation  (A-IF. 


E = 


£ (a,  F,  (aT.Xi)  + a3FT(Xj)-Yi). 


The  following  notation  is  used  for  simplilication: 


F,i  for  F,  (aT.Xj).  Ft;  for  Ft  (Xj), 


FPj  for 


9F  j(aT,Xj) 
9aT 


The  derivatives  of  E are; 
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In  this  relation  F||.  Fij.  FPj  arc  an  array  ot  values.  Fij  is  indepeiuienl  of  the  aj  parameters.  This 
array  is  computed  only  once.  F|j  and  FPj  are  functions  of  a i and  should  be  recomputed  at  each 
iteration.  The  array  FPj  is  computed  from  Fjj  values: 


FP: 


3 E‘  j (a .\j I 
3a  -i 


:X|‘  F,  ta^.Xj) 


A.4.2  Convergence 


The  convergence  algorithm  used  to  find  the  minimum  of  the  function  F is  given  in  figure  .\-l . 
The  algorithm  computes  the  function  and  its  derivatives  with  a given  set  of  parameter  values.  The 
parameters  are  then  increased  or  decreased  depending  on  the  sign  of  the  corresponding  derivative,  Fhe 
increment  step  is  initially  defined  to  be  10  i of  the  parameter  value.  The  following  rule  is  used  to  vary 
the  step  size:  il  in  two  consecutive  steps  the  corresponding  derivative  does  not  change  sign,  the  step 
size  is  doubled  (up  to  a maximum  of  20';  of  the  parameter  value);  if  the  derivative  does  change  sign, 
the  step  size  is  halved.  Convergence  is  stopped  after  a given  number  of  iterations. 

In  the  algorithm  of  figure  .X-l . the  previous  step  parameters  are  denoted  .\j.  the  correspond- 
ing function  and  derivatives  are  denoted  FA  and  3F.\  3A^.  and  B is  used  for  next  iteration  value. 

Aj  is  the  increment  size. 


APPENDIX  B 


SEA-STATE  DATA  REDUCTION 


B.  1 MEASURED  PARAMETERS 

The  sea-state  measurements  available  from  the  buoy  can  be  listed  as  follows: 

Compass  bearing 
Roll  angle 
Pilch  angle 

Acceleration  along  three  axes 
Wave  staff  heiglit  measurements. 

ITiese  variables  are  discussed  below  in  terms  of  their  code  and  calibration  properties. 


B.  I . I Magnetic  Compass  Bearing 


The  compass  bearing  is  obtained  as  a parallel  digital  output.  The  code  is  a seven-bit  Gray  code 
shown  in  table  B-l,  with  one  division  equivalent  to  approximately  3*^  (i.e.,  360/128).  The  time  con- 
stant of  the  compass  is  roughly  3 sec.  The  reference  line  on  the  buoy  is  in  the  pitch  plane,  as  shown 
In  figure  B-l.  Clockwise  rotation  results  in  an  increasing  angle. 


B.  1.2  Roll  Angle 


Roll  angle  is  available  as  a 1 2-bit  output  from  a synchro/digital  converter.  The  roll  angle  is 
offset  so  that  in  a calm  sea  the  digital  output  is  2048.  i.e.,  an  offset  binary  code  with  one  least  signifi- 
cant bit  corresponding  to 


A = 360  X 2"  * “ degrees 
= 5.4  minutes. 

Depression  of  the  buoy  port  side  (relative  to  the  reference  line)  results  in  an  increasing  roll  angle 
indication. 

The  roll  (and  pitch)  angles  are  derived  from  displacement  gyros  that  are  gimballed  to  indicate 
angular  displacement  from  horizontal  or,  more  precisely,  angular  displacement  from  the  long-term 
(5-min)  average.  Figure  B-l  indicates  the  orientation  of  pitch  and  roll  angles. 
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1 

1 

0 

0 

0 

0 

0 

64 

96 

180 

1 

1 

0 

0 

0 

0 

1 

65 

97 

183 

1 

1 

0 

0 

0 

1 

1 

66 

99 

186 

1 

L 

0 

0 

0 

1 

0 

67 

98 

188 

1 

1 

0 

0 

1 

1 

0 

68 

102 

191 

1 

1 

0 

0 

1 

1 

1 

69 

103 
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64 

32 

Weights 
16  8 4 

2 

1 

Order 

Weighted 

Decimal 

Degrees 

Code 

1 

1 

0 

0 

1 

0 

1 

70 

101 

197 

1 

1 

0 

0 

1 

0 

0 

71 

100 

200 

1 

1 

0 

1 

1 

0 

0 

72 

108 

202 

1 

1 

0 

1 

1 

0 

1 

73 

109 

205 

1 

1 

0 

1 

1 

1 

1 

74 

111 

208 

1 

1 

0 

1 

1 

1 

0 

75 

110 

211 

1 

1 

0 

1 

0 

1 

0 

76 

106 

214 

1 

1 

0 

1 

0 

1 

1 

77 

107 

217 

1 

1 

0 

1 

0 

0 

1 

78 

105 

219 

1 

1 

0 

1 

0 

0 

0 

79 

104 

222 

1 

1 

1 

1 

0 

0 

0 

80 

120 

225 

1 

1 

1 

1 

0 

0 

1 

81 

121 

228 

1 

1 

1 

1 

0 

1 

1 

82 

123 

231 

1 

1 

1 

1 

0 

1 

0 

83 

122 

233 

1 

1 

1 

1 

1 

1 

0 

84 

126 

236 

1 

1 

1 

1 

1 

1 

1 

85 

127 

239 

1 

1 

1 

1 

1 

0 

1 

86 

125 

242 

1 

1 

1 

1 

1 

0 

0 

87 

124 

245 

1 

1 

1 

0 

1 

0 

0 

88 

116 

247 

1 

1 

1 

0 

1 

0 

1 

89 

117 

250 

1 

1 

1 

0 

1 

1 

1 

90 

119 

253 

1 

1 

1 

0 

1 

1 

0 

91 

118 

256 

1 

1 

1 

0 

0 

1 

0 

92 

114 

259 

1 

1 

1 

0 

0 

1 

1 

93 

115 

262 

1 

1 

1 

0 

0 

0 

1 

94 

113 

264 

1 

1 

1 

0 

0 

0 

0 

95 

112 

267 

1 

0 

1 

0 

0 

0 

0 

96 

80 

270 

1 

0 

1 

0 

0 

0 

1 

97 

81 

273 

1 

0 

1 

0 

0 

1 

1 

98 

83 

276 

1 

0 

1 

0 

0 

1 

0 

99 

82 

278 

1 

0 

1 

0 

1 

1 

0 

100 

86 

281 

1 

0 

1 

0 

1 

1 

1 

101 

87 

284 

1 

0 

1 

0 

1 

0 

1 

102 

85 

287 

1 

0 

1 

0 

1 

0 

0 

103 

84 

290 

1 

0 

1 

1 

1 

0 

0 

104 

92 

292 
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64 

32 

Weights 
16  8 4 

2 

1 

Order 

Weighted 

Decimal 

Degrees 

Code 

1 

0 

1 

1 

1 

0 

1 

105 

93 

295 

1 

0 

1 

} 

1 

1 

1 

106 

95 

293 

1 

0 

1 

1 

1 

1 

0 

107 

94 

301 

1 

0 

1 

1 

0 

1 

0 

108 

90 

304 

1 

0 

1 

1 

0 

1 

1 

109 

91 

306 

1 

0 

1 

1 

0 

0 

1 

110 

89 

1 

0 

1 

1 

0 

0 

0 

111 

88 

312 

1 

0 

0 

1 

0 

0 

0 

112 

72 

315 

1 

0 

0 

i 

0 

0 

1 

113 

73 

316 

1 

0 

0 

1 

0 

1 

1 

114 

75 

321 

1 

0 

0 

1 

0 

1 

0 

115 

74 

323 

1 

0 

0 

1 

1 

1 

0 

116 

78 

326 

1 

0 

0 

1 

1 

1 

1 

117 

79 

329 

1 

0 

0 

1 

1 

0 

1 

118 

77 

332 

1 

0 

0 

1 

1 

0 

0 

119 

76 

335 

1 

0 

0 

0 

1 

0 

0 

120 

68 

338 

L 

0 

0 

0 

1 

0 

1 

121 

69 

1 

0 

0 

0 

1 

1 

1 

122 

71 

343 

1 

0 

0 

0 

1 

1 

0 

123 

70 

346 

1 

0 

0 

0 

0 

1 

0 

124 

66 

349 

1 

0 

0 

0 

0 

1 

1 

1'’5 

67 

352 

1 

0 

0 

0 

0 

0 

1 

126 

65 

354 

1 

0 

0 

0 

0 

0 

0 

127 

64 

357 

A 


z 


COMPASS 

REFERENCE 

ROLL  AXIS 


Figure  B- 1.  Buoy  Roll  and  Pitch  Coordinate  System 

B.  I.3  Pitch  Angle 

This  measurenieiit  is  similar  to  that  for  roll  angle.  The  angular  deviation  from  horizontal  is 
output  as  an  offset  binary  number  ( 1 2 bits).  The  resolution  is  5.4  min.  As  shown  in  ligure  B-l . pitch 
and  roll  are  defined  relative  to  horizontal,  so  that  depression  of  the  pitch  line  causes  a positive  increase 
in  pitch  angle. 

Note  that  satisfies  the  usual  definition  of  roll  axis,  i.e..  the  axis  about  which  rotation  takes 
place.  Similarly,  y satisfies  the  usual  definition  of  pitch  axis,  with  x-z  referred  to  as  the  pitch  plane. 


V 


B.  1.4  Three-.Axis  Acceleration 


Strap-down  accelerometers  provide  the  three  components  of  acceleration  of  the  buoy.  Fhe 
output  in  each  case  is  a I 2-bil  word.  olTset  so  that  zero  acceleration  indicates  204S.  I he  sensitivity  ol 
each  accelerometer  can  be  summarized  in  the  following  way: 


Vertical  acceleration  (zl  +()..?g 


Resolution 


0 (ig 

40<)(, 


B-t< 


Horizontal  acceleration  (x  and  y)  +0.5  g 


Resolution  = ■ 

4096 

It  is  assumed  that  the  x-  and  y-axes  lie  in  the  pitch  and  roll  planes,  respectively. 


I 

B.  1.5  Wave  Staffs 


Wave  height  relative  to  the  buoy  is  sensed  by  means  of  a lO-elemeni  array.  The  array  elements 
are  numbered  and  spaced  as  shown  in  figure  B-2.  with  the  reference  axis  (numbers  1.3.  2)  lying  in  the 
pitch  plane  of  the  buoy.  i.e..  orthogonal  to  the  y-axis  shown. 

All  the  variables  are  sampled  at  a 25-Hz  sampling  rate  and.  with  the  exception  of  the  compass, 
all  sensors  are  capable  of  25-Hz  response. 

The  wave  staffs  are  variable-resistance  devices,  with  resistance  decreasing  as  a function  of 
immersion  depth.  .Ideally,  the  resistance  value  should  be  sensed  electrically  by  driving  the  wave  staff 
with  a constant-current  source  (ac)  and  measuring  voltage.  However,  a simpler  approach  is  to  employ 
a constant-voltage  source  with  series  resistance  (see  fig.  B-3).  When  the  series  resistance  is  con- 
siderably larger  than  the  maximum  wave  staff  resistance,  the  net  effect  is  equivalent  to  a constant- 
current  generator.  For  smaller  values  of  R^.  the  wave  staff  voltage  can  be  expressed  as 


V 


s 


where: 

Vj  = measured  staff  voltage 
Rj  = wave  staff  resistance 
Rq  = series  resistance 
V = voltage  source. 

In  processing  the  sea-state  data,  a relationship  between  the  measured  voltage  and  wave  height 
is  needed.  Thus. 


Ro  V-V, 


B-7 


Wave  Staff  Diameter  = 0.128  in. 


Figure  B-2.  Wave  Staff  Geometry 
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Figure  B-3.  Wave  Staff  Circuitry 

Assuming  a linear  relationship  between  height  h and  resistance  ratio  R 

s o 


h = 


we  find  that 


h 


v-v/ 


The  resistive  values  of  interest  are  kncwn  to  be; 


R = 4i,i,n 

O 1 

Rs 


= S(>4.:5n  . 


This  maximum  valiio  OLX'urs  with  the  wave  stall  out  ol  the  water.  With  the  stall  partially  immersed,  it 
is  assumed  that  that  portion  of  the  winding  below  the  surface  is  totally  short-circuited.  With  lull 
immersion,  the  value  of  should  be  zero.  The  length  of  the  wave  staff  is  3 I in.  (78.8  cm),  so  the 
constant  a is  computed  to  be: 


0.788  X 4d9 
8b4.25 


meters 


= 0.45, 


and.  hence. 


h = 0.45 


,V-V, 


It  should  be  noted  that  three  analog  calibration  signals  are  supplied  from  the  buoy.  These 
signals  are  generated  by  the  same  means  as  the  wave  staff  voltages  except  that  a dummy  resistor  is  used 
in  place  of  a physical  wave  staff.  These  resistors  were  chosen  to  simulate  the  resistance  oi'a  wave  staff 
when  it  is  extended  by  roughly  lO'T,  50'7,  and  90'7  of  full  scale.  In  practice,  the  resistance  values 
corresponded  to  wave  staff  extensions  of  1 1.57C7. 45.36'^.  and  86.78'7.  The  corresponding  voltages 
will,  of  course,  not  indicate  this  same  percentage  of  full  scale  because  of  the  nonlinear  relationship 
between  resistance  and  measured  voltage. 

By  observing  the  calibration  signals  and  maximum  staff  resistance  signals  while  the  buoy  is 
on  the  Coast  Guard  cutter  deck,  it  is  possible  to  provide  a check  on  the  unknown  constants  in  the 
expression  for  R.^.  Four  points  are  plotted  in  figure  B-4  corresponding  to  the  voltage  V^,  that  results 
for  1 \7r.  469f,  86''.7.  and  100'7  of  the  maximum  R^.  The  voltage  scale  is  expressed  as  an  integer  where 
the  maximum  value  of  4096  corresponds  to  the  A/I)  converter  reference  voltage. 

With  the  aid  of  these  plotted  points,  it  is  found  that 


3910  (same  units  as  V^) 


and 


= 0..59 


B-IO 


riius. 


with  expressed  as  an  integer  in  the  range  0 to  40‘)5. 

The  enLiivalent  value  of  R,,/R^  obtained  from  knowledge  of  the  resistor  values  is 
" ^inax 


= 0.578. 

R 

^inax 


With  the  measured  full-scale  value  of  2440,  V is  computed  to  be 


V = 3860. 


B.2  AN.XLOG  TAPt  FORMAT 

The  basic  data  frame  consists  of  the  compass,  roll,  pitch,  acceleration,  and  wave  staff  samples 
to  which  are  appended  three  calibration  samples  indicating  voltages  representative  of  107?,  50'’?,  and 
90'7  of  the  full-scale  analog  voltages  in  the  system  under  normal  operation.  A 33-bit  sync  train  is  also 
added  at  the  head  of  each  frame.  Fhe  composite  trame  is  indicated  in  table  B-2.  Kach  frame  is 
extracted  from  the  analog  tape  in  serial  form,  with  most  significant  bits  of  each  sample  appearing 
first.  In  the  analog-to-digital  tape  conversion,  the  software  is  set  up  to  handle  I 2-bit  bytes  of  data, 
rite  total  data  segment  of  223  bits  is  therefore  extended  to  cover  five  bits  of  the  sync  code  to  give  a 
multiple  of  I 2;  i.e.,  the  equipment  is  set  up  to  search  for  the  abridgeil  sync  code  (octal); 


00  3 6 0 7 0 6 2 2 (28  bits). 


i 

I 

I 

! 


TABLE  B-2.  DATA  FRAME  FORMAT 


01703431L25 


L).-vr.\ 

BUS 

DA  FA 

BUS 

Compass 

7 

^1 

12 

Roll 

12 

^2 

12 

Pitch 

12 

A 

12 

^10 

12 

z 

A 

12 

Calj^ 

12 

X 

A 

12 

Cal2 

12 

y 

Cal^ 

12 

I'he  rciiiainini>  Hvl’  hits  are  then  treated  as  data  and  appear  in  tlie  five  most  significant  bit 
locations  of  the  1 3-bit  compass  word.  The  actual  compass  bearing  must,  (herelore.  be  extracted  by 
first  eliminating  these  bits  and  then  decoding  the  seven-bit  (iray  code. 


B..^  CONVHRSION  OF  ANALCX'.  l APFS  TO  DIOITAL 

The  FMR  tape-stripping  facility  at  TSC  was  used  to  convert  the  analog  tapes  to  digital  form. 

.Analog  tape  speed  was  3-.^/4  ips  and  the  data  rate  was  (i.4  Kbps.  I'he  FMR  bit  and  frame  synchroni/er 
were  loaded  with  the  following  contents  (octal): 

Bit  Synchroni/er 

Register 

1 

s 

4 

.s 

B-I.I 

1 

I 


Contents 

I 4 4 
0 I 
: : I 
0 4 3 
0 3 0. 


1 


Frame  Synchronizer 


Register 

I 

1 

} 

4 

5 

6 

7 

8 
9 

10 


13 

14 

15 

16 
17 


Contents 

0  0 0 
I 7 0 

1 6 I 

0 1 7 

2 0 7 
2 2 () 

1 5 5 
0 04 
000 
0 0 0 
0 0 0 
0 7 0 
0 2 3 
0 4 I 
0 0 0 
14  1. 


The  program  DFS.LOAI  I 2. 1 I was  loaded  into  the  PDP-1  I by  means  of  the  absolute  loader  and 
the  desired  parameters  were  entered.  A concise  summary  of  the  parameters  required  is  given  in  section 
3.2,  along  with  a listing  of  the  digital  tape  file  structure  for  each  analog  tape  stripped.  The  data  seg- 
ments are  referenced  by  experiment  time  in  the  form:  hours:  minutes:  seconds. 


B.4  PROChSSlNC.  OF  THFl  JOINT  WAVF  SLOPF  PROBABILITY  DFNSITY  FllNC  lTON 

In  this  section  three  topics  relative  to  the  measurements  leading  to  the  joint  probability 
density  functions  for  the  wave  slopes  are  discussed: 

a.  Bias  removal 

b.  Measurement  geometry 

c.  Coordinate  rotation . 

In  item  c,  transformation  to  a fixed  coordinate  system  is  required  because  of  the  nonnegligible 
time  variability  of  the  magnetic  compass  bearing.  I'lie  fixed  reference  was  chosen  to  be  in  the  direc- 
tion of  the  satellite  so  that  the  results  could  be  easily  related  to  the  multipath  scattering  measurements 
for  in-plane  and  cross-plane  IJights. 
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B.4. 1 Bias  Removal 


During  the  programming,  a ciuestion  arose  relative  to  the  removal  ol'  bias:  At  what  point  in 
the  processing  should  the  bias  be  removed?  We  chose  to  remove  the  bias  Irom  the  raw  roll  and  pitch 
data;  this  choice  guarantees  that  an  unbiased  angle  0 with  value  at  or  near  0“  is  mapped  into  a slope 
value  via  that  part  of  the  tan  0 curve  near  the  origin. 


B.4,2  Measurement  Geometry 

The  conventions  used  to  describe  roll  and  pitch  are  illustrated  in  Hgure  B-5a.  The  x-axis  is 
the  reference  direction  and  points  to  the  “bow.”  The  negative  y-axis  points  in  the  starboard  direction. 
When  the  +x  axis  is  depressed  in  the  x-z  plane,  this  is  considered  to  be  an  increase  from  0.  i.e.,  a 
positive-valued  pitch  angle  6.  Positive  values  of  0.  the  roll  angle,  result  from  depression  of  the  +y  axis 
(port  arm)  in  the  y-z  plane.  * 

The  sense  of  angular  dependence  is  appropriately  illustrated  in  figure  B-5a.  We  point  out  that 
all  the  tape  data  on  sea  slope  are  referred  to  the  buoy  coordinate  system  according  to  the  convention 
illustrated  in  figure  B-5a. 


B.4.3  Coordinate  Rotation 

The  fact  that  the  buoy  coordinate  system  itself  rotates  is  one  major  reason  for  the  need  to 
rotate  each  data  point  to  a single  fixed  coordinate  system.  Because  of  convenience,  the  aircraft 
heading  is  chosen  as  the  reference  direction.  The  transformations  required  are  illustrated  in  figure 
B-5b.  To  summarize  briefly,  note  that  all  sea-slope  measurements  arc  referenced  to  the  buoy  coordi- 
nate system,  which  has  a nonconstant  bearing  a.  It  is  desired  to  reference  these  quantities  to  the 
constant  aircraft  heading  at  13.  To  do  this,  a coordinate  transformation  through  the  variable  angle 
£ = ^ - a is  required.  Using  primed  variables  to  indicate  the  location  of  a point  in  the  “new”  coordi- 
nate system,  the  transformation  is  governed  by: 


X = x'  cos  £ + y'  sin  £ 
y = y'  cos  £ - x'  sin  £ 

z = z' . ( B- 1 ) 

'Actually  the  pitch  depression  is  an  increase  above  180'*  ( I 72.7®  was  measured  for  an  8®  lift ).  where- 
as a depression  of  the  starboard  arm  causes  a decrease  in  0 relative  to  180®  (a  I 2®  depression  re- 
sulted m a measured  value  of  lb7..S®).  We  have  subtracted  the  180®  from  the  data,  however. 


a l>  « > H > iu> 


(b)  Transformation  to  Fixed  Reference  System 

ROLL  ANGLE 
PITCH  ANGLE 
AIRCRAFT  HEADING 
MAGNETIC  COMPASS  READING 


Figure  B-S.  Sea-State  Geometries 


B-16 


In  the  "old"  cooidiiiatf  sy^tcnl.  the  oqualiun  lor  a plane  passing  through  the  origin  was 
governed  by: 


Ax  + By  + ( / = 0. 


Using  the  tact  that  a veetor  with  eomponenis  A.  B.  and  ('  is  nomial  to  the  plane,  it  is  an  easy  matter 
to  show  that 


Additionally,  one  can  show  by  using  the  translormation  properties  in  equation  (B-l ) that  the  equation 
lor  the  plane  in  the  "new"  coordinate  system 


A'x'  + B'y'  = 0, 


has  the  coelTicients 


A'  = A cos  { - B sin  { 


B'  = A sin  { + B eos  J 


(•'  = ( 


Relations  ( B-.U  and  ( B-4)  hold  lor  the  primed  eoordin.ite  system  as  well  as  the  unprimed.  By 
substituting  equation  B-(>  into  ( B-.U  and  (B-4),  it  is  easily  shown  that  the  roll  and  pitch  angles  in  the 
new  coordinate  system  are  governeil  by 


Ian  0'  = Ian  0 sin  ^ + Ian  0 eos  { , 


B-l  7 


.Iiul 


Ian  o'  = tan  0 cos  J - tan  0 sin  ^ . (B-S) 

-\s  an  important  example,  we  look  at  the  ease 

t]  = :7o‘’. 


For  titis  case,  sin  { = - cos  a and  cos  { = - sin  a.  so  that  finally  we  have 


tan  0'  = - tan  0 cos  q - tan  0 sin  a , 


* 

and 


tan  O'  = tan  0 sin  a + tan  0 cos  a . 


(B-IO»  » 


the  sign  assignment  in  equations  (B-‘>)  and  ( B-IO)  can  be  checked  quickly  by  drawing  a few  sketches 
and  testing  for  a = 0^’  and  Q = 

Summari/ing.  we  note  from  figure  B-5  that  the  slopes  in  the  x-  and  y-directions  are  given  by 
r)^i  = tan  O'  and  rj^.,  = tan  0'.  Thus,  in  the  final  reference  system  (positive  x-axis  aligned  with  aircraft 
heading),  the  wave  slopes  are  given  by: 


= tan  0 cos  J - tan  0 sin  J 


\ 


rj^.i  = tan  0 sin  J + tan  0 cos  { 


{ = ti-a 
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